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ON THE IMPOSSIBILITY OF ALL OVERLAPPING TRIPLET CODES 
IN INFORMATION TRANSFER FROM NUCLEIC 
ACID TO PROTEINS 


By 8S. BRENNER 


MEDICAL RESEARCH COUNCIL UNIT FOR THE STUDY OF THE MOLECULAR STRUCTURE OF 
BIOLOGICAL SYSTEMS, CAVENDISH LABORATORY, CAMBRIDGE, ENGLAND 


Communicated by G. Gamow, June 10, 1957 


It is a generally accepted view that nucleic acids control the synthesis of pro- 
teins, and it has been proposed more specifically that the sequence of amino acids 
in a polypeptide chain is determined by the order of nucleotides in ribo- or deoxy- 
ribonucleic acid. The problem of how this determination is effected has come 
to be known as the ‘‘coding” problem. The formal aspects of this problem can be 
investigated theoretically, and most of the work done in this field has recently 
been reviewed by Gamow, Rich, and Yéas.! 

Since there are only four different nucleotides in RNA or DNA to determine 
twenty different amino acids, it is clear that more than one nucleotide must be used 
to code for each amino acid. Most codes have been constructed on the basis that 
each amino acid is determined by a set of three nucleotides. Such triplet codes, 
however, have an excess of information, since there are sixty-four different triplets 
for the twenty amino acids. In Gamow’s original diamond code,’ several triplets, 
chosen in a particular way, coded for any given amino acid; the code was therefore 
“degenerate.”” This code was also of the overlapping type—that is, the number 
of nucleotides in the nucleic acid was equal to the number of amino acids in the 
polypeptide chain. Gamow’s diamond code does not, in fact, code for known se- 
quences, and the same is true for the major-minor code, another overlapping triplet 
code, invented by L. Orgel.! These are, however, only two examples of a large 
number of possible codes of this type which can be obtained by choosing different 
ways of degenerating the triplets. To test all of these systematically is clearly 
impossible, and hence it is necessary to have some general theorem about such 
codes. 

The general overlapping triplet code has the following properties. 

(i) The coding triplets are chosen from four nucleotides, A, B, C, and D, giving 
sixty-four different triplets. 

(ii) Coding is overlapping, each triplet sharing two nucleotides with the sueceed- 
ing triplet in a sequence. Thus the sequence ABCDA codes for three amino acids: 
ABC for the first, BCD for the second, and CDA for the third. 

(iii) An amino acid may be represented by more than one triplet; that is, the 
sixty-four triplets are degenerated into twenty sets. 

687 
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Since any dipeptide sequence is represented by a sequence of four nucleotides, 
there cannot be more than 256 different dipeptides. On the other hand, if all di- 
peptide sequences were possible, 400 would be expected. Thus overlapping codes 
introduce restrictions in amino acid sequences. The number of dipeptide se- 
quences known is less than 256, and, although statistical studies have suggested 
that all dipeptides are likely to be found, the significance of this’ result has been 
difficult to assess.!_ The sample of proteins studied is highly selected, a large number 
of sequences are fragmentary, and the methods used to study sequences further 
bias the data. 

However, sufficient sequences are known to prove that it is impossible to code 
them with overlapping triplets. The proof is simple and does not depend on any 
special way of degenerating the triplets. It consists in the demonstration that 
sixty-four triplets are insufficient to code the known sequences. 

Proof: Since successive triplets share two nucleotides in common, any given 
triplet can be preceded by only four different triplets and succeeded by only four 
different triplets. In an amino acid sequence j.k.l., we call 7 an N-neighbor, and / 
a C-neighbor, of k. For every four different N-neighbors (or C-neighbors) or part 
thereof, k must have one triplet assigned to it. Thus the minimum number of 
triplet representations for each amino acid can be counted from a table of neighbors. 

The available sequences are given in the Appendix. From these sequences a 
grid is constructed and the different neighbors counted for each amino acid. The 
number of triplets assigned to each amino acid is based on the larger number of its 
neighbors. These data are given in Table 1, from which it can be seen that seventy 

TABLE 1 
Minimum Minimum 
No. of No. of 


Amino Triplets Amino Triplets 
Acid C-Neighbors N-Neighbors Required Aci C-Neighbors N-Neighbors Required 


Lys 17 5 Pro 12 4 
Ser 7 13 5 Tyr 10 3 
Gly 15 : Glu 11 3 
Leu 15 . Glun 9 3 

3 


Cys i . Asp i 
Arg . j é Asn ¢ 10 3 


Ala 4 f d lleu ¢ 3) 3 


Val é : His j 9 3 
Thr : Met 5 7 2 
Phe : Try ‘ 3 1 

Total 70 


triplets would be required to code the sequences. We conclude, then, that all 
overlapping triplet codes are impossible. 

This result has one important physical implication. The original formulation of 
overlapping codes was based on the similarity of the internucleotide distance in 
DNA to the spacing between amino acid residues in an extended polypeptide 
chain. It was supposed that each amino acid was spatially related in a one-to-one 
way with each nucleotide on a nucleic acid template. The present result shows 
that this cannot be so and that each amino acid is stereochemically related to at 
least two, if not three, nucleotides, depending on whether coding is partially over- 
lapping or nonoverlapping. The difficulties raised by this can easily be overcome 
by assuming that the polypeptide sequence is in contact with the nucleic acid tem- 
plate only at the growing point, and detailed schemes can be readily proposed. 
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As far as the coding problem is concerned, it now appears that all amino se- 
quences are likely to be found and that it will not be possible to effect a “decoding” 
by discovering restrictions in sequences. The nonoverlapping of triplets implies 
that there must be some way of determining which triplets in a sequence are coding 
triplets and which are not, and a very interesting code has recently been proposed 
by Crick, Griffith, and Orgel,’ in which this problem is dealt with in a novel manner. 


APPENDIX 
Amino Acip SEQUENCES 


In writing the sequences, the same conventions used by Gamow et al.' have been followed. 
Wherever doubt exists as to whether glutamic acid is present as such (glu) or as the amide 
(glun), it has been assigned as ‘“‘glux,” and the same rule has been followed for aspartic acid 
and asparagine. All the longer lysozyme sequences suggested by Thompson (Biochem. J., 
60, 507; 61, 253, 1955) have been omitted, since some of these appear to be incorrect when 
compared with those established by the French workers. Sequences established by carboxy- 
peptidase digestions alone are given at the end of the list but are omitted from the grid. 
The same applies to the pepsin sequence of Williamson and Passmann (J. Biol. Chem., 222, 
151, 1956), as there are conflicting reports about the N-terminal group (Van Vunakis and 
Herriott, Biochim. et Biophys. Acta, 23, 60, 1957). 

The grid (Table 2) shows the number of times dipeptide sequences are found. Identical 
sequences from the closely related proteins vasopressin and oxytocin and corticotrophin and 
melanophore-stimulating hormone are only recorded once. Dipeptide sequences from 
lysozyme are not recorded if the same sequence is found in a longer peptide. When both 
glu and glun are absent, glux is counted as a neighbor, and the same rule is followed for asp, 
asn, and ax. 

SEQUENCES USED IN THE GRID* 
-thr. gly. ileu’- 
-thr. ser. ileu® ¢- 
-ala. gly. val?- 
Insulin A: Gly. ileu. val. glu. glun. cys. cys. ala. ser. val*. cys. ser. leu. tyr. glun. leu. glu. 
asn. tyr. cys. asn.*~® 


Insulin B: Phe. val. asn. glun. his. leu. cys. gly. ser. his. leu. val. glu. ala. leu. tyr. leu. val. 
cys. gly. glu. arg. gly. phe. phe. tyr. thr. pro. lys. ala.* ® 


Oxytocin: Cys. tyr. ileu. glun. asn. cys. pro. leu. gly. NH2.’ 


-arg*- 


Vasopressin: Cys. tyr. phe. glun. asn. cys. pro. lys?. gly. NH2.’ 


Corticotrophin: Ser. tyr. ser. met. glu. his. phe. arg. try. gly. lys. pro. val. gly. lys. lys. arg- 
arg. pro. val. lys. val. tyr. pro. asp. gly. ala. glu.® asp. glun. leu. ala®. glu. ala. phe. pro. leu- 
-glux. ala. ser‘- 
glu. phe.~'8 


Glucagon: His. ser. glun. gly. thr. phe. thr. ser. asp. thr. ser. lys. tyr. leu. asp. ser. arg. 
ala. glun. asp. phe. val. glun. try. leu. met. asn. thr." 
-ser’- 
: . b . 
Melanophore-stimulating hormone: Asp. glu’. gly. pro. tyr. lys. met. glu. his. phe. 
gly. ser. pro. pro. Lys. asp.!—!" 


** Cattle. ° Pig. “Sheep. “ Whale. * Horse. ‘Salmon. % Chicken. * Man. 
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-ileu’- 
Hypertensin: Asp. arg. val. tyr. val*. his. pro. phe. his. leu.'*: 3 
-ser*- 


f 


. glun. cys. his. thr. val. glu. lys-’ 


Cytochrome ec: -val. glun. lys. eys. ala® ® © 


99 


Trypsinogen: Val. asp. asp. asp. asp. lys. ileu. val. gly-?* 
Ribonuclease: Lys. glu. thr. ala. ala. ala. lys. phe. glun. arg. glu-*4~* 
-tyr. cys. asn. glun. met. met. lys. ser. arg. asn. leu. thr. lys. asp. arg. cys.*4~”? 


-lys. asn. val. ala. cys. lys. asn. thr-** 27 


-cys. asn. arg. glu. ser. thr. ser. gly. lys. tyr. pro. asn. ala. cys. tyr. lys. thr. thr. asn 


ala. Lys. his-*6 27 

-tyr. glun. ser. tyr-*4 

-phe. asp. ala. ser. val.?4, °8. 29 

Lysozyme: -lys. asx- 
-arg. his. lys-*! 
Lys. val. phe. gly. arg- 
-ala. lys. phe. glux- 
-asx. tyr. arg. gly- 
-arg. gly. tyr. ileu. leu- 
-asn. ala. tyr. gly. ser. leu. asn- 
-thr. pro- 
-leu. pro-** 
-asn. arg- 
-ileu. arg- 
-thr. pro. gly. ser. arg- 
-val. ala. try. arg-** 
-gly. cys. arg. leu. 
-phe. glu. ser. phe. asp. glu. ala. thr. asp. arg- 
-cys. glu. ala. leu. ala. ala. met. lys. arg-** 
-ala. ala- -asx. ileu- -ileu. arg- 
-ala. leu- -asx. leu- -ileu. asx- 
-ala. lys- -cys. ala- -ileu. val- 
-ala. met- -cys. arg- -leu. ala- 
-arg. aSX- -CyS. asx- -leu. cys- 
-arg. Cys- -cys. glux- -leu. leu- 
-arg. gly- -cys. ileu- -lys. gly- 
-arg. leu- -cys. lys- -met. asx- 
-asx. ala- -gly. leu- -met. lys- 
-asXx. arg- -gly. met- -phe. asx- 
“ASX. ASX- -glux. ala- -phe. glux 
-asx. gly- -glux. leu- 


Ovalbumin: -ala. gly. val. asx. ala. ala-* 
-asx. ser. glux. ileu. ala- 
-glux. ser. ala-*’ 
-cys. ala- -thr. cys- -cys 
-cys. gly- -val. cys- -gly. 
-cys. val- -cys. glux- -asx. 
-ser. Cys- -phe. cys- 

-val. ser, pro.” 
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Papain: eu. pro. glux-* 
ser. &sx- -CYS. ASX- -val. cys-“ 
-ASX. CYS- -cys. gly. 
Hemoglobin: Val. glun. leu- 
Val. leu-*! 
-lys. arg- leu- -val. lys- 
-lys. leu- lys- -arg. lys- 
-ser, arg- phe- -tyr. arg- 
-ala. arg- arg- -phe. arg-*” 
Myoglobin: Gly. leu-** “4 
y-Globulin: Ala. leu. val. asx. glux-* 
B-Lactoglobulin: -val. glux- -thr. lys- -ala. lys- 
-val. leu- -lys. gly- -leu. lys- 
-asx. lys- -pro. lys- -phe. lys-” 
-glux. lys- -lys. pro- 
Carboxypeptidase: Asn. ser- 
-ser. thr-* 
-thr’- 


Serum albumin: Asx. ala®-" 
Lactogenic hormone: Thr. pro. val. thr. pro-” 


Tobacco mosaic virus: -thr. ser. gly. pro. ala. thr.*! 
Pro. ileu. glux-** 
Casein: Lys. leu. val. ala. glux. asx-** 
-leu. gly- -ser. pro- _—_-ser. leu-” 
Chymotrypsinogen: -gly. leu. ser. arg. ileu. val-** % 
-tyr. thr. asn. ala-* 
-gly. asp. ser. gly-*” 


SEQUENCES NOT USED IN THE GRID 


Serum albumin: -gly. val. ala. leu’. 
-ser. val. thr. leu. ala. ala*.® 


Actin: -his. ileu. phe.* 
Tropomyosin: -ala. ileu. met. thr. ser. ileu.® 
Pepsin: Leu. gly. asp. asp. his. glu- (ef.5) 
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INFECTION OF PROTOPLASTS BY DISRUPTED T2 VIRUS 
By JouHNn SpizizEn* 


DEPARTMENT OF MICROBIOLOGY, WESTERN RESERVE UNIVERSITY SCHOOL OF MEDICINE, 
CLEVELAND, OHIO 


Communicated by H. G. Wood, May 29, 1957 
INTRODUCTION 


The functions of the components of virus particles have been the subject of much 
investigation in two widely different virus systems. In the case of tobacco mosaic 
virus (TMV) it has recently been demonstrated that the ribonucleic acid isolated 
from the virus retains infectivity® ? and possesses the property of genetic deter- 
mination.* In contrast with the low infectivity of TMV, bacteriophage T2 can 
be studied with very high precision. The protein coat of T2 isolated after osmotic 
rupture has been found to adsorb to the receptor site on the cell wall of the host 
and to stop multiplication of the host cell.‘ After injection through the tail of T2 
virus into the host cell, the deoxyribonucleic acid (DNA) component is found to be 
associated with the replicating activity of the virus.» However, no direct biological 
activity of the isolated DNA has been demonstrated. It would be of much interest 
to determine whether the biological activity of the protein moiety subsequent to 
attachment*: ° is a necessary prerequisite for the function of the DNA portion or, 
indeed, whether DNA can be active in the complete absence of protein. 

Since cell wall components and the protein tip of T2 virus interact prior to injec- 
tion of DNA, it was considered possible that the removal of the cell wall might per- 
mit the replicating moiety of the virus to enter the cytoplasm directly. Removal of 
cell walls can be achieved by lysozyme treatment of Escherichia coli B in hyper- 
tonie sucrose to produce forms known as “‘protoplasts.’’” * Intact T2 virus will not 
infect these protoplasts. We have attempted to infect protoplasts directly with DNA 
isolated from T2, without success. However, when osmotically disrupted prepara- 
tions of T2 virus* were mixed with protoplasts of EF. coli B, significant increases in 
the numbers of infective particles were observed. In fact, similar increases in infec- 
tive units were obtained with protoplasts of a number of bacteria whose cells are nor- 
mally resistant to infection. Although the osmotically shocked preparations em- 
ployed in these studies contained all the components of the original virus, we have 
obtained evidence to indicate that the activity observed may be due to DNA units 
protected by a protein shell. This communication describes some of the data ob- 
tained on the interaction of protoplasts and disrupted T2 virus, as they relate to the 
functional requirements of virus-cell interaction. 
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METHODS 


Preparation of Protoplasts.—Protoplasts of £. coli B and other bacterial strains 
were prepared by a modification of a procedure suggested by Roberts,’ similar to 
that of Zinder and Arndt.’ The bacteria were grown in Difeo nutrient broth for 
3 hours at 37° C. on a reciprocal shaker, starting with 2 X 10’ cells per milliliter. 
The cells were then harvested by centrifugation and resuspended in the same volume 
of 0.5 M sucrose and 0.01 M tris(hydroxymethyl)aminomethane at pH 9, contain- 
ing 12.5 wg. lysozyme per milliliter. The mixture was incubated for 5 minutes at 
37° C. and then centrifuged, and the sediment was resuspended in the same volume 
of 0.25 M sucrose containing nutrient broth at one-half concentration. The proto- 
plast spheres so produced were enumerated in a Petroff-Hauser chamber, using the 
phase-contrast microscope. With E. coli B under the conditions described, an 
average of 4 X 10° spheres per milliliter were produced from 2 X 10° viable bac- 
teria per milliliter. 

Preparation of T2 Shockates.—Phage strain T2r+ was propagated on F. colt B 
in a glucose-ammonium salts medium’? and concentrated by differential centrifuga- 
tion in an SS-1 Sorvall centrifuge, first at 3,000 * g to remove cellular debris and 
then at 20,000 X g for 45 minutes. In the preparation of osmotically shocked 
suspensions, phage concentrates were diluted to 2 X 10'! plaque-forming particles 
per milliliter in 317 NaCl for 5 minutes at room temperature and rapidly diluted 
40-fold in 1 per cent gelatin solution, resulting in 99 per cent or better reduction in 
infectivity. T2 shockates heated to 71° C. for 15 minutes or T2 phage suspen- 
sions heated to 71° C. for 15 minutes followed by osmotic shocking were also em- 
ployed in these studies. The latter preparations will be referred to as “heated 
shockates.”’ 

“Infection” of Protoplasts with Heated Osmotic Shockates of T2.—It was found that 
heated preparations of T2 shockate which contain <10? infectious particles per milli- 
liter as determined by plaque counting with whole cells produced infective T2 when 
incubated with protoplasts of 2. coli B or EF. coli B/2. T2 phage increases as high 
as 10° plaque-forming particles per milliliter were observed in 60 minutes at 37° C. 
Infective phage formation was not observed when heated shockates were mixed 
with protoplasts which were suspended in broth or water without sucrose.'! Nor 
were any viable phage formed when suspensions of whole cells with or without 
sucrose were used. The possibility of some type of “multiplicity reactivation’’!* 
with intact cells was ruled out by using diluted suspensions of cells; in no instance 
was any viable T2 formed. 

It was’ of some interest. to investigate in detail the production of virus in proto- 
plasts from a resistant strain, since reinfection with any intact virus of cells re- 
maining after lysozyme treatment cannot occur. L. coli B/2 was used for these 
studies, and it was found that the general characteristics of viral formation were 
similar to those with Z. coli B. Figure 1 illustrates the kinetics of phage formation 
in E. coli B/2. Aliquots of the mixture of protoplasts and heated T2 shockate 
were removed at intervals and plated with sensitive 2. coli B. Although it is not 
certain that all “intraprotoplast” infective particles were thus being assayed, 
some idea of the formation of infective particles was obtained. Thus, within the 
first 5 minutes, we observe the appearance of 2 X 10* T2 phage per milliliter, which 
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may be due to disaggregation and release of intact T2, since in a number of prep- 
arations made in identical fashion we could show by dilution and incubation in the 
absence of protoplasts that approximately this quantity of intact virus was present. 
Nevertheless, it represents only a very small fraction of the phage subsequently 
formed. This is followed by a lag period up to 11 minutes, a period of rapid multi- 
plication lasting a little over 20 minutes, and a period of diminished multiplication 
after 30 minutes, yielding 7 X 10’ T2 units per milliliter. 
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Fig. 1.—Formation of infectious T2 phage from heat-treated and 
osmotically shocked T2 (noninfectious preparation) in protoplasts of 
E. coli B/2 at 37° C. 





No quantitative estimates of the rates of adsorption and the yields of virus per 
protoplast have as yet been made. Experiments were performed in which mixtures 
of protoplasts and heated shockate were centrifuged after incubation for 5 minutes 
at 37° C. and the protoplasts resuspended in 0.25 M sucrose and broth. The in- 
fective virus T2 produced after 60 minutes at 37° C. was then determined and found 
to be nearly the same as that formed in the uncentrifuged mixtures. This indicates 
that rapid adsorption of the disrupted virus to the protoplast occurs under these 
conditions. 

Certain conditions have been found essential for infective T2 virus synthesis. 
Thus versene (0.01 M) completely abolished the activity, and Mg ions (0.001 M) 
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gave some stimulation. Sodium chloride at concentrations of 0.1 M or greater were 
found to inhibit. Virus formation was completely prevented in the absence of 
casein hydrolyzate or nutrient broth. Finally, certain metabolic poisons have been 
found to inhibit, namely, potassium cyanide, potassium arsenite, and potassium 
arsenate at 0.001 M concentration or higher. 

Characteristics of Heat-treated T2 Shockate-—When heat-treated T2 shockate was 
incubated with crystalline trypsin (Nutritional Biochemicals; 50 ug/ml) for 30 
minutes and soybean trypsin inhibitor (Worthington Biochemical Company; 1 mg/ 
ml) added prior to mixing with protoplast suspensions, only 0.1 per cent as much 
of the infective T2 was produced as compared with that formed in the absence of 
trypsin treatment (Table 1). Crystalline desoxyribonuclease (DNAase) (Worth- 


TABLE 1 
Inrectiviry oF HEATED T2 SHockaTEes ON PROTOPLASTS 
or E. coli B 
——-+ E. coli B PrororpLtastst———— 
Initial Infectivity Titer/Ml, 
TREATMENT OF Heatep SHocxkate* Titer/MI 60 Min. 
None <10? 8 X 10’ 
50 ug. trypsin, 30 min., 37°, + 1 mg. antitrypsin <10? 2X 104 
50 wg. DNAase + 0.005 M Mg, 30 min., 37°, + 0.01 M 
citrate <10? 2 xX 108 
50 ug. DN Aase + 0.01 M citrate, 30 min., 37° <10? 1 X 10° 
Versene, 0.01 M <10? <10' 
Homogenized in Waring Blendor or Mickle agitator, 2 min. <10? <108 
Ultraviolet, 1-mm. layer, 30 inches from GE 15-watt lamp, 
5 see. <10? <108 


* Heated shockate: T2 virus at 10? particles per milliliter, heated to 71° C., 15 minutes; dissolved in 3 M 
NaCl, 5 minutes, and rapidly diluted 1:40 in 1 per cent gelatin. Infectivity titer, <10? plaque-formers per milli- 


jiter. 
t Protoplasts: 3-hour broth growth of EZ. coli B, treated with 12.5 ug/ml lysozyme in 0.5 M sucrose, centrifuged 
and resuspended in 0.25 M sucrose-broth. 


ington Biochemical Company) and a highly purified pancreatic DNAase (kindly 
supplied by Dr. Werner Baumgarten, Merck and Company) destroyed less of the 
activity than did trypsin; however, a comparable destruction occurred in the pres- 
ence of citrate, which would inhibit the enzymatic action on free desoxyribonucleate. 
The effect of DNAase on heated shockate may therefore be due to contaminating 
enzymes. Both trypsin and DN Aase had little or no inhibitory action when added 
10 minutes after protoplasts and heated shockate were mixed. 

The activity of heated shockate was completely inhibited by homogenization 
for brief intervals in a chilled Waring Blendor or by agitation in the Mickle disin- 
tegrator (Table 1). In addition, complete inactivation was effected by a 5-second 
exposure in a 1-mm. layer to a GE 15-watt germicidal ultraviolet lamp at 30 inches, 
under conditions where little T2 phage inactivation could be observed. 

These data indicate that “heated shockate” contains disrupted viral particles 
capable of inducing virus synthesis in the presence of protoplast suspensions. The 
sensitivity of these disrupted viral forms to the action of trypsin contrasts with the 
resistance of normal viral particies. That the nucleic acid of the subparticles which 
participate in viral synthesis is probably in a more exposed form than in intact 
virus T2 is indicated by the extreme lability to ultraviolet and homogenization. 
It is tempting to visualize that in the process of osmotic disruption and heat treat- 
ment of T2, packets of DNA are produced which are kept intact by the protein 
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layer of the head membrane, which serves to keep the DNA from forming co- 
polymers. ‘This is further suggested by the fact that in experiments with DNA (at 
| mg/ml) isolated from T2 virus (using the phenol extraction procedure of Gierer 
and Schramm! and the acetate-duponol method of Mayers and Spizizen"*), it was 
found that no phage synthesis occurred on mixing with protoplasts. Discrete 
packages of DNA may then be required for infection of protoplasts, since DNA 
has been clearly shown to be required for T2 phage synthesis.° 

Protoplast “Infection’”’ with Unheated Osmotic Shockates of T2.—Phage increases 
have also been observed when an unheated osmotic shockate of T2 was mixed with 
protoplasts derived from E. coli B. Table 2 shows data from an experiment in 
which, after incubation at 37° C. for 60 minutes, an increase of 7 X 10" infective 
units per milliliter was observed. When citrate was added to DN Aase-treated 
T2 shockate or soybean trypsin inhibitor to trypsin-treated shockate prior to mixing 
with protoplasts, no inhibition of virus formation occurred, unlike the effects pro- 
duced with heated shockate. The addition of crystalline desoxyribonuclease, 
crystalline ribonuclease, or crystalline trypsin to the mixture of shockate and 
protoplasts completely inhibited the increase in infective particles. Thus the 
effects of the enzymes are probably on protoplasts (cf. Spiegelman) and in this 
respect. differ from intact-cell experiments in which virus growth is unaffected by 
these enzymes. ‘T2 virus which was not shocked but was diluted to the same in- 
fectivity titer present in shockate did not multiply in the presence of protoplasts. 

Centrifugation experiments with T2 shockate using the Spinco ultracentrifuge 
(Table 2) showed that the supernatant fluid, after centrifugation at 50,000 x g 


TABLE 2 
INFECTION OF E. coli B Prorop.asts witn T2 Osmotic 
SHockaTes (UNHEATED) 
-—— + E. coli B Prororptastst———— 
Initial Increase in 
Infectivity Infectivity/M1l, 
TREATMENT OF Osmotic SHOCKATE* Titer/Ml 60 Min. 
None 5 10° 7 X10? 
DNAse, 100 ug/ml, + 0.005 M MgS0,, 15 min., 37° 5 10° 0 
DNAase, 100 ug/ml, + 0.01 M citrate 5 108 7X 10° 
RN Aase, 100 ug/ml, 15 min., 37° 5 10° 0 
Trypsin, 100 ug/ml, 15 min., 37° : 108 - 
Trypsin, 100 ug/ml, 15 min., 37°, + 1 mg. antitrypsin 10 7x10’ 
Supernatant of Spinco centrifugation, 50,000 X g, 60 min. 105 x 10° 
Phage control, diluted in 1 per cent gelatin 10° 0 
* Osmotic shockate: T2 virus at 2 X 10!! particles per milliliter in 3 M NaCl for 5 minutes at room temper- 
ature; rapidly diluted 40-fold in 1 per cent gelatin. 


+ Protoplasts: 3-hour broth growth of Z. coli B, treated with 12.5 ug/ml lysozyme in 0.5 M sucrose, centrifuged 
and resuspended in 0.25 M sucrose-broth. 


for 60 minutes, retained only a fraction of its original activity. Thus the supernate 
contained 4 per cent of the infectious virus, and the increase in T2 on incubation 
with protoplasts was only 5.8 per cent of that obtained with the uncentrifuged 
shockate. The bulk of the activity could be found in the sediment. Centrifuga- 
tion of shockate at 20,000  g for 30 minutes, however, did not sediment the bulk 
of the activity observed with protoplasts, although more than 90 per cent of the 
intact virus was thus sedimented. Thus the activity with protoplasts is probably 
due to disrupted particles with similar but not identical sedimentation behavior to 
intact phage. 
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Virus T2 Formation in Protoplasts of ‘Resistant’? Bacteria.—Of much interest 
was the observation that disrupted T2 viral preparations could produce infectious 
virus in protoplasts from bacterial species which are resistant to viral infection. 
Table 3 summarizes the results."° The resistance of those cells in the intact state 


TABLE 3 
PHaGe T2 PropuctTion FrroM INACTIVATED SHOCKATE IN PROTOPLASTS 
oF DiFFERENT BACTERIAL SPECIES 
Sensitivity of Intact Phage T2 Formed/M1 
Cells to Infective T2 in 60 Min. from Inactivated 
Organism Phage (104/M1) T2 Shockate* in Protoplastst 
E. coli B + (>) 108 
E. coli B/2 0 (>) 108 
). coli K12 + 5 x& 103 
E. coli Crookes + 1.5 x 105 
Anaerogenic E. coli (WR-2) 0 104 
Atypical Klebsiella (WR-3) 0 5 xX 108 
Providence group paracolon (WR-6) 0 1.3 X 10' 
A. aerogenes 417 0 2 X 106 
Proteus vulgaris (WR-1) 0 0 
B. megaterium 899 0 0 
Pseudomonas saccharophila 0 0 


* T2 heated to 71°, 15 minutes, and shocked; infectivity for Z. coli B intact cells, < 10* plaque-forming particles 
per milliliter. 
t Protoplasts prepared from 3-hour growths in nutrient broth, as deseribed in text. 


whose protoplasts were found to produce virus is apparently due to lack of proper 
receptor sites in the cell walls. We have also found a number of resistant organisms 
whose protoplasts will not support virus T2 synthesis from heated shockate, namely, 
the Proteus, Bacillus megateritum, and Pseudomonas saccharophila strains tested. 
The nature of the resistance of the latter organisms remains to be determined. 


DISCUSSION 


Multiplication of bacteriophage was observed by other investigators in protoplasts 
of B. megaterium,": ® B. subtilis,® and EF. coli,’ provided that infection with virus 
was performed prior to removal of surface components containing the specific 
receptors for viral adsorption. Direct infection of protoplasts with intact virus 


has not been observed. 

The present experiments show that disrupted particles of T2 virus can produce 
infectious virus in protoplasts from sensitive and certain resistant strains of bacteria. 
By the removal of cell wall components, it is possible to determine whether the 
protoplasm is or is not able to support viral propagation. Thus Aerobacter aerogenes 
417 is highly susceptible in its protoplast state to coliphage T2 shockates, whereas 
the intact cells are completely resistant to infection with T2 virus. The specific 
cell receptor—viral tail combination is thus replaced by a less exacting protoplast- 
disrupted viral interaction. 

Two types of disrupted T2 preparations have been used in these studies of pro- 
toplast infection. Unheated osmotic shockate of T2 is unaffected by trypsin or 
desoxyribonuclease. The heated osmotic shockate apparently contains more 
degraded forms of T2 virus, since trypsin almost completely destroys its ability to 
infeet protoplasts. This indicates the requirement of some protein for the pro- 
duction of T2 phage in protoplasts. This latter protein may be part of an envelope, 
protecting the DNA from forming large polymers. The formation of large molecu- 
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lar aggregates of DNA, as occurs during its isolation free of protein, may prevent 
its entrance into protoplasts. On the other hand, the protien may be involved 
in attachment to the protoplast membrane or may actually participate in the 
formation of viral particles. 

Intact protoplastic structure appears to be essential for phage production, since 
“shocking”’ of protoplasts by removal of hypertonic sucrose prevents the formation 
of infective virus. The requirement for amino acids and the inhibitory action of 
metabolic poisons suggest that complex biochemical events are probably involved 
in viral formation from shockate in protoplasts. 

The data presented do not establish definitely that new viral particles were 
being synthesized by some form of replication. In most experiments, in fact, 
less virus was recovered than was originally employed for the preparation of dis- 
rupted viral particles. In certain experiments, however, it has been possible to 
recover more infective T2 after incubation of heated T2 shockate with protoplasts. 
This was also found to be the case when protoplasts of resistant strains were em- 
ployed, where infection of whole cells with T2 virus could not account for any of 
the increases observed. Nevertheless, the possibility must be considered that 
some form of reactivation might have occurred, involving replacement or synthesis 
of some portion of the virus, which then makes it possible to infect intact sensitive 
cells. It is also possible that some type of multiplicity reactivation may have 
taken place in the protoplasts, involving the participation of a number of damaged 
particles. Since phage multiplication is known to occur in protoplasts, our results 
may be accounted for by viral reactivation within protoplasts with subsequent 
replication prior to liberation. Studies are in progress to elucidate this problem. 

Thus the use of protoplasts provides us with a technique for the study of the 
“minimal unit” capable of viral formation. 


SUMMARY 


Disrupted T2 virus preparations possessing no measurable infectivity have been 
found to produce infectious T2 virus after incubation with protoplasts from LE. coli 
B or a number of other bacteria normally resistant to T2 infection. Virus pro- 
duction occurs only in intact protoplasts, that is, in hypertonic sucrose, and requires 
amino acids. Although protein-free DNA isolated from T2 virus cannot produce 
infective virus, evidence is presented that DNA units protected by a protein shell 
may be involved in the production of T2 virus in protoplasts. 
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INCORPORATION OF RADIOACTIVE CARBON INTO A PTERIDINE 
OF A BLUE-GREEN ALGA 


By C. Van Baaten, H. 8S. Forrest, AnD J. Myrrs* 


LABORATORY OF ALGAL PHYSIOLOGY AND GENETICS LABORATORY, DEPARTMENT OF ZOOLOGY 
UNIVERSITY OF TEXAS, AUSTIN, TEXAS 


Communicated by Theophilus S. Painter, June 10, 1957 


In a preliminary paper! we have reported a striking relationship between the loss 
of photosynthetic activity in a blue-green alga, Anacystis nidulans, and the release 
from the algal cells of a group of closely related pteridines. Two of these com- 
pounds make up the greater bulk of the fluorescent materials released; one is a 
yellow crystalline compound which was described briefly in our previous paper.! 
The second is a blue fluorescent compound (hereinafter called “Compound C’’) 
which appears to be the principal material obtained from the cells when they are 
killed with acid. 

Compound C has been obtained chromatographically pure in small quantities, 
and some observations have been made on its chemistry, although these have 
necessarily been limited by the quantities so far available. In the first place, it 
is readily oxidized by alkaline permanganate to 2-amino-4-hydroxy-6-carboxypteri- 
dine (J in the accompanying diagram), as previously reported, and it has ultra- 
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violet and infrared absorption spectra (Figs. 1 and 2) very similar to that of 2- 
amino-4-hydroxy-6-(1',2’-dihydroxypropy!l) pteridine (biopterin; /7 in the dia- 
gram).2* It is readily distinguished from the latter compound, however, by its 
Rf values in various chromatographic solvents (Table 1). Compound C is quite 
stable to alkali, being recovered unchanged after 2 hours of heating with 2 NV 
sodium hydroxide at 100° C. On the other hand, when treated with 2 N hydro- 
chloric acid under the same conditions, it is degraded quantitatively to a different 
fluorescent substance which proved to be identical, paper-chromatographically 
and spectrophotometrically (ultraviolet and infrared), with biopterin. These 
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facts, then, suggest that Compound C contains an aldehyde linked in an acetal 
grouping across the vicinal hydroxy! groups of biopterin. So far, we have not been 





/\ 


ee a 


wcugent sod x 


. 


‘, 
a 
iY 
‘ 











sa 
iL 
240 © —285 0 0 400 
WAVE LENGTH (my) 
Fic. 1.—Ultraviolet absorption spectrum of Compound C: 
-X —~ X = in 0.1 N sodium hydroxide: X --- X --- X 
= in 0.1 N hydrochloric acid. 


As T 





TT T T 


1 we 


T T T T T T T 
i 4 1 4 1 4 4 4 
Ld 2 








1 
“ so 
WavE LENGTH © MECRONS 
Fic. 2.—Infrared spectrum of Compound C. Potassium bromide wafer containing approxi- 
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able to identify the fragment which is released in this hydrolysis. It appears to be 


a low-molecular-weight compound which is not acidic (since Compound C is vir- 
tually neutral, as shown by paper electrophoresis). 

Although the structure of Compound C is not yet fully defined, we wish to report 
a series of experiments in which it has been shown that, over a short period of time, 


TABLE 1 
Solvent Biopterin Compound C 
0.31 0.10 4 per cent sodium 0.60 
citrate 
0.37 sec-Butanol, formic 0.64 0.47 
acid, water 


(8:2:5) 


Solvent Biopterin Compound C 

n-Butanol, acetic 0.73 
acid, water (4:1:1) 

n-Propanol, 1 per 0.46 
cent ammonia 


(2:1) 
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radioactive C'4 as C'4O, is incorporated by A. nidulans into the acid hydrolyzable 
fraction of Compound C but not into the 1’,2’-dihydroxypropy! side chain. 

Experimental.—Freshly harvested cells grown in a continuous culture chamber* 
were washed and suspended in phosphate buffer (pH 7.5; 0.02 M). Eighteen milli- 
liters of cell suspension containing approximately 3 mm*/ml of cells were placed in 
a special, flat Warburg vessel (total volume, 60 ml.; area, 57 cm.*), fitted with a 
Siamese side arm. Approximately 2 mg. of BaC'Os; (0.36 mc/mg) were placed in 
one of the side arms and 0.5 ml. of 10 per cent hydrochloric acid in the other. The 
vessel was uniformly illuminated by a neon grid of 4,100 wwatts/em* (<7,000 A). 
The cell suspension was shaken for 5 minutes in the light, and the C'*O, was then 
released in the side arm. After a further 5 minutes of shaking in the light, the 
cells were killed by adding the remaining acid from the side arms. To insure kill- 
ing and release of the pteridines, the flask was heated for 2 minutes at 100° C., after 
which time the contents were filtered through Celite, to give a colorless, acidic 
extract. 

The fluorescent materials in this extract, previously adjusted to pH 4, were ad- 
sorbed on a Magnesol’ column (1 X 15 cm.), which was then washed very thoroughly 
with water, followed by 5 per cent aqueous acetone. At this point fluorescent 
materials were beginning to appear in the eluate, and elution of these was then 
completed with 30 per cent aqueous acetone. Combined fractions containing 
Compound C were evaporated to small bulk, and the concentrated solution was 
applied to a single spot on a sheet (46 X 56 cm.) of Whatman No. 3 mm. filter 
paper. The chromatogram was developed in the first direction with sec-butanol, 
formic acid, water (8:2:5) and in the second direction with n-propanol, 1 per cent 
ammonia (2:1). The fluorescent material corresponding to Compound C was 
eluted from the paper with water and rechromatographed against authentic Com- 
pound C in three more solvents, its radioactivity and quantity being measured 
after each chromatogram (see Table 2). The quantitative retention of radio- 
activity with fluorescence after each step is taken as convincing evidence of the 
purity of the radioactive Compound C. 


TABLE 2 


MEASUREMENT OF CoMPOUND C AND OF Its Rapioactivity* 
Recovered 
Total Compound C 
Counts/Min. Measured 
Recovered Fluorometri- 
Corrected for cally 
Sample Solvent Background (ug.)T 
sec-Butanol, formic acid, water (8:2:5) 161 5.10 
i Sir ry 1 per cent ammonia (2:1) 156 5.10 
n-Butanol, acetic acid, water (4:1:1) 112 3 56 
Biopterin recovered n-Butanol-acetic acid, water (4: 1:1) 3 2.00 
after hydrolysis 
Dark sec-Butanol, formie acid, water (8:2:5) 0 3.08 
and n-Propanol, | per cent ammonia 
(2:1) 


* All chromatcgrams were ascending, on Whatman No. 3 mm. Fluorescence measurements were done with a 
Farrand Fluorometer Model A. Radioactivity was determined with a Tracerlab windowless flow counter, SC-16, 


and a Tracerlab Utility Scaler. 
+ Estimated against a weighed standard of Compound C. 


Light 


After the final chromatogram, Compound C was hydrolyzed in 2 N hydrochloric 
acid (0.2 ml.) for 30 minutes at 100°, and, after removal of the acid over potassium 
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hydroxide in a desiccator, the residue was chromatographed against authentic 
biopterin. As shown in Table 2, the recovered biopterin was not radioactive. 

This experiment has been carried out three times in all, using various different 
procedures for the purification of the Compound C (e.g., Filtrol* columns, instead of 
Magnesol, used before or after paper chromatography; different solvent combina- 
tions for paper chromatography, including the use of dilute acetic acid as an alter- 
nate solvent); the results in each experiment were the same as those recorded 
above. 

When an exactly similar experiment was carried out with cells incubated in the 
dark, no radioactivity was found in the Compound C recovered (Table 2). 

Discussion.—From the work of our previous paper! and that reported herein, the 
following facts may be marshaled. From A. nidulans and several other blue-green 
algae a number of pteridine compounds can be extracted in quantities of about 0.1 
per cent of the cell dry weight. Variation in distribution of the various pteridines 
with the method of extraction suggests that they result from one or two parent 
compounds in the cells. Concomitant loss of pteridines and photosynthetic ac- 
tivity during incubation of Anacystis at 4° C. suggests a role in photosynthetic 
metabolism. 

Nathan, Hutner, and Levin’ recently have shown for the parasitic trypanosomid 
Crithidia fasiculata a growth-factor requirement which can be supplied by biopterin 
independently of folic acid. It was postulated that “Crithidia factor,” an uncon- 
jugated pteridine, identified with biopterin, provides for a metabolic function other 
than that performed by folic acid. Our data raise the question whether biopterin 
itself is the metabolically active compound in the cells. No significant quantity 
of biopterin is extractable from Anacystis, although at least one of the extractable 
pteridines is readily converted to biopterin by acid hydrolysis. On the other hand, 
our data provide direct evidence for the postulate that compounds of the biopterin 
type are active metabolically in a role different from those of the folic acid group. 

Of the several pteridines extracted from Anacystis, one, Compound C, is con- 
verted by acid hydrolysis to biopterin and a fragment which has not yet been 
identified. Since the ultraviolet and infrared spectra of biopterin and Compound 
C are almost identical, the unidentified fragment must be attached to the side 
chain of biopterin and must be relatively inert spectrophotometrically. During 
photosynthetic metabolism of Anacystis, C'4 from COs is rapidly incorporated into 
the unidentified fragment of Compound C but not into the biopterin portion (cf. 
Table 2). 

Since the identity of the C'*-labeled fragment of Compound C is unknown, there 
remains the possibility that the recovered radioactivity is due to an impurity 
which carries through all the chromatographic separations and finally separates 
chromatographically only after acid hydrolysis. We consider this possibility 
highly unlikely in view of the strict correspondence between the radioactivity and 
measured fluorescence of Compound C in successive chromatograms. 

The metabolic locus of action of Compound C (or its cellular parent compound) 
is not specified by the evidence at hand. Labeled carbon is incorporated into the 
unidentified fragment only in the light. However, Anacystis is an obligate photo- 
autotroph with a very low rate of endogenous (dark) metabolism. The question 
whether the locus of action is associated with the early stages of carbon dioxide 
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fixation and reduction is not answered by the evidence presently available. A 
cautious hypothesis must be limited to this statement: A pteridine related to 
biopterin acts as a metabolic carrier for a carbon fragment. 

Work is in progress on the complete identification of Compound C and its rela- 
tion to the other pteridines previously described in Anacystis. Such information 
should allow a more precise statement of the metabolic role of these compounds. 

Summary.—A labeled pteridine has been obtained from A. nidulans during 
photosynthesis in the presence of C'4OQ.. After acid hydrolysis, this pteridine 
yielded unlabeled biopterin and a labeled and as yet unidentified fragment. The 
possible metabolic significance of this and related compounds is discussed. 

* This work was supported by the Rockefeller Foundation and by the Robert A. Welch Founda- 
tion, Houston, Texas. 
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IONIC SEMIPERMEABILITY AS A BULK PROPERTY* 
By Spencer L. Barro, Jr., GEoRGE KARREMAN, HELMUTH MUELLER, AND ALBERT 
SZENT-GYORGYI 
INSTITUTE FOR MUSCLE RESEARCH, MARINE BIOLOGICAL LABORATORY, WOODS HOLE, MASSACHUSETTS 


Communicated May 23, 1957 


One of the most fundamental facts of biology is the unequal distribution of K+ 
and Na* in and outside the cell, K+ being found in a relatively high concentration 
inside and Nat outside. This unequal distribution is attributed to the function of 
the cell membrane, to its selective permeability and “pumping action.” No gen- 
erally accepted picture has yet been developed for the mechanism of this hypotheti- 
cal activity of the membrane, and so it seems justifiable to consider alternate pos- 
sibilities. 

The water molecule has such a distribution of charges, including lone pairs of 
electrons, as to lend it a tendency to link molecule to molecule in a tetrahedral ar- 
rangement. In ice this regular association extends throughout the whole structure. 
In liquid water the structure is broken up by heat agitation into small fragments, 
and there is no long-range regularity. The fragments of the crystalline structure 
are very small and have a temporary existence, associations being constantly formed 
and broken. 

There is strong evidence that around surfaces the small crystallites may join to 
more extensive lattices, establishing a long-range regularity reaching many hundreds 
of molecules deep into the fluid.'| So it is not impossible that most of the water 
within densely packed biological structure is in a more or less ordered state. 
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While in K+ the nuclear charge is screened by three electron shells, in Na* it is 
screened only. by two. Moreover, the dimensions of K+ are very similar to those of 
the water molecule, which is not the case with Na*. Consequently we can expect 
that K+ does not disturb the regularity of the water structure, while Na* does and, 
by ordering the water dipoles around itself, creates disorder in the water structure; 
K+, on the other hand, may readily be substituted for water molecules. This dif- 
ference is borne out by the temperature dependence of the conductivity of NaCl 
and KCl as well as by the “B factor” of both cations.2,. When Na* is introduced 
into lattice-ordered water structure, bonds have to be broken which hold the water 
molecules together, and work will have to be done. Correspondingly, we can also 
expect the system to tend to eject Na* already introduced and to balance its negative 
charges with preference by K+. This allows us to construct a theory which picture 
the ionic semipermeability of cells as a bulk property, while earlier theories ascribed 
semipermeability to the ‘“‘“membrane”’ only. 

Such a theory suggests experiments and calculations along different lines. First, 
we can attempt to demonstrate the deranging activity of Na* and the inactivity of 
K+. Second, we can try to demonstrate actually the preference of biological struc- 
tures for K+ in the absence of a membrane. Third, we can calculate whether the 
dipole forces of water are sufficient to account for the known distribution of K+ and 
Na* in biological material. 

It has been suggested by one of us! that the phosphorescent behavior of various 
fluorescent substances in a frozen watery solution depends on the orderly structure 
of water and its interaction with the excited molecule. It has also been shown that 
Na* causes a strong disturbance in the light emission of Rhodamin B, while K* 
does not doso. Equally inactive is NH{, the physical properties and dimensions of 
which are very similar to those of K+, while Li+ was found to create even graver 
disturbance than Nat. The reaction with Rhodamin was too complex to allow a 
simple interpretation. The situation is simpler with tryptophan, adenosine tri- 
phosphate, and nucleic acid. All these substances are excited to triplets’ and show 
in the presence of 0.25 per cent glucose a long-lived phosphorescence if their watery 
solution is frozen in dry ice and illuminated by ultraviolet light. If their triplet 
excitations are stabilized by the water structures, then the afterglow must be 
quenched by any substance which interferes with the regularity of these structures, 
and we have to expect Nat to cause quenching, while K* is inactive. The experi- 
ment showed that KCl, in physiological concentrations such as 0.1 M, had no effect 
on the lifetime of the triplet excitation of the above systems and did not materially 
reduce the duration of the afterglow, while NaCl, in the same concentration, had a 
strong quenching action. 

For the demonstration of the preference of cellular structures for K+ in the ab- 
sence of an active membrane, we used glycerinated muscle fibers. According to our 
experience, glycerination destroys the “membrane,” as indicated by the complete 
absence of electric excitability and conduction. Moreover, in the presence of the 
high K* concentrations employed, the membrane, even if present, could be expected 
to be inoperative. So we tried to decide in the experiment whether such a glycer- 
inated material is capable of distinguishing between K+ and Nat. The difficulty of 
the experiment lies in the fact that we can expect undisturbed water structures in 
resting muscle only. and the charge distribution in this material must be greatly in- 
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fluenced by the tetravalent ATP anions present. However, it is very difficult to 
produce a glycerinated muscle which is completely relaxed in the presence of ATP. 
To do this, the so-called ‘‘relaxing factors’? have to be present and the glycerol 
treatment has to be of short duration. The objection may be raised that in such a 
carefully treated material the “membrane”? may not completely be destroyed. 

Our experiments were performed in the following manner: 2-3 mm. thick bundles 
of the psoas muscle of freshly killed rabbits were tied to applicator sticks and im 
mersed for 3-6 hours in 50 per cent glycerol at 0° C. Then the muscle fascicles were 
transferred into a muscle extract prepared according to Bendall.4 The muscle 
fascicles remained in this extract for another 2 hours, in order to provide them with 
a sufficient amount of relaxing factor to prevent their contraction in the next bath, 
which consisted of the same extract diluted with an equal volume of salt solution. 
Here the final concentration of KC] and NaCl was 0.08 M, that of MgCl, was 0.004 
M, and that of ATP was4mM. After 1-hour incubation at room temperature the 
muscle was decomposed into 0.1—0.2-mm. thin fiber bundles which were left in this 
liquid for another 20-40 minutes. After this time a number of fiber bundles (20-40) 
were liberated from adhering fluid by placing them on filter paper, then dissolved in 
concentrated nitric acid and analyzed for K and Na on the flame photometer. The 
results are summarized in Table 1. As the table shows, no significant difference 
was found between the relative concentrations of K and Na in the bath and in the 
fibers. 


TABLE 1 
Experiment K/Na Ratio K/Na Ratio 


No. in Bath in Muscle 


I 0.96 (iid 
1.14 
II 1.09 1.12 


III 1.12 1.10 


As is well known, Ca plays a dominant role in muscle function and ionie distribu- 
tion. In order to maintain its normal ionic K/Na balance, the muscle must be 
provided with Ca?*. For this reason, the above experiments were completed by 
adding CaCl, in 0.5 mM final concentration to the bath. In the first experiment 
the intracellular K/Na ratio rose within 5 minutes from 1.07 to 1.31. In the second 
experiment the quotient rose within 3 minutes from 1.07 to 1.33. After 4 minutes 
it went up to 1.8. In the third experiment the quotient rose within 4 minutes from 
1.10 to 1.52 and within 9 minutes to 1.77. In a fourth experiment the fibers were 
immediately placed in a Ca-containing bath; the K/Na ratio was found within 4 
minutes to be 1.30, in 12 minutes, 1.35, while in the bath it was 1.05. 

On more prolonged incubation the quotient dropped again, which may have been 
due to the disorganizing effect of the contracture induced by the addition of Ca. 
It seems that two effects were counteracting one another: the Ca tended to increase 
the quotient of K/Na, while the developing contraction tended to equalize the 
ionic distribution. It is unfortunate that no relaxed state could be maintained in 
the presence of Ca. 

The experiment is a rather imperfect one, and the quotients obtained are still far 
from the physiological quotient 10; nevertheless, the experiment definitely shows an 
unequal distribution, a trend toward eliminating Na* and retaining K*, 
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The results discussed suggested an attempt to calculate the free energy of the 
ordering of water molecules. As is well known, the membrane potential of cells is 
about 0.1 volt. Therefore, the question arises whether the free energy of the order- 
ing of water molecules is equal to or somewhat greater than 0.1 electron volt, which 
free energy would be necessary, according to the proposed mechanism, to produce 
the unequal ionic distribution found in cells. 

Quantitative treatment of the ordering aspects of molecules in condensed phases 
is not easy. Whereas the distances between the molecules are largely determined 
by attractive and repulsive central forces, orientational forces also play an important 
role for the ordering of the molecules. To obtain an estimate of the free energy due 
to the orientational forces, we will consider only the interaction of nearest neighbors, 
for which we take only dipole-dipole interaction. Recently, Pople® has studied, on 
the basis of statistical mechanics, systems with noncentral force fields and has been 
able to derive an expression for the orientational free energy. His expression re- 
duces, for the case of point-dipole interaction, to 


Fo” vn l B («) 
ny 6 kT \a?/ ’ 


where F'°” represents the orientational free energy per molecule, 2 the number of 
nearest neighbors of a molecule, u the dipole moment of a molecule, a the nearest- 
neighbor distance, 7’ the absolute temperature, and k the Boltzmann constant. 

Since the water molecules tend to arrange the other water molecules around 
themselves in a tetrahedron, we will take z = 4. The value of a in water is 2.9 A 
(Morgan and Warren®). Substituting these values together with » = 1.85 & 107'8 
esu (Sanger and Steiger’) and the value of k, we find, at 7 = 300° F°” = —0.193 
electron volt (which, of course, is of the order of magnitude of the energy of the H- 
bond). It therefore seems that there is enough energy already in the locally 
ordered structures of molecules in liquid water, even taking into account only 
orientational forces and dipole-dipole interactions, to sustain electrical potentials of 
the order of 0.1 volt. If the order were more extensive than in liquid water, the 
energy would certainly be enough to account for the order of magnitude of bio- 
electrical potentials. 

* This research was supported by a grant from the Commonwealth Fund, the National Heart 
Institute (H-2042 R), the Muscular Dystrophy Associations of America, Inc., the National 
Science Foundation, the American Heart Association, the Association for the Aid of Crippled 
Children, and the United Cerebral Palsy Foundation. 
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THE EFFECT OF GIBBERELLIN UPON FLOWER FORMATION* 


By ANTON LANG 
DEPARTMENT OF BOTANY, UNIVERSITY OF CALIFORNIA, LOS ANGELES 
Communicated by W. H. Chandler, May 6, 1957 


As reported in a preliminary note,' gibberellin was found to induce flower forma- 
tion in a cold-requiring plant, a biennial variety of Hyoscyamus niger (henbane). 
This result was of interest in two respects: Gibberellin had been known only to 
promote shoot growth by cell elongation but not to have any effect on the initiation 
of flowers or other plant organs; chemical induction of flower formation has been 
largely an unsolved problem in plant physiology. 

The studies on flower induction by gibberellin have since been extended to a 
number of other plants. The purpose of this paper is to summarize the results. 
In the meantime, this line of investigation has been taken up by other authors, too, 
and a definite picture is beginning to emerge about the effect of gibberellins in 
flower formation. 


MATERIAL AND METHODS 


Plants were chosen for this study in which flower formation is controlled either 
by low temperature or by photoperiod and which therefore can be maintained in 
the vegetative condition at will. Following is a list of the species and varieties 
that were used. 

a) Biennials.—Plants with an obligatory cold requirement in flower formation: 
Hyoscyamus niger L. (henbane), two biennial strains; Daucus carota L. (carrot), 
varieties Early French Forcing and Danver’s Half Long; Petrosilenum crispum 
(Mill.) Mansf. (parsley), varieties Hamburg and Plain; and Brassica napus L. 
(turnip), variety Purple Top Globe. 

b) Winter Annuals.—Flower formation hastened by cold treatment but may 
ultimately take place without it: Secale cereale L. (Petkus winter rye). 

c) Long-Day Plants.—Hyoscyamus niger L. (henbane), an annual variety; 
Crepis tectorum L.; Samolus parviflorus Raf. (water pimpernel); Silene armeria 
(garden or Sweet William catchfly) ; Secale cereale L. (Petkus spring rye). Hyoscya- 
mus, Crepis, Samolus, and Silene, when grown on short day, remain vegetative 
for indefinite periods of time; spring rye initiates flowers on long and short day, 
but in the latter case with great delay. 

d) Short-Day Plants.—Glycine max (L.) Merr. (soybean), variety Biloxi; Xan- 
thium pensylvanicum Wallr. (cocklebur). 

The plants were grown and treated in the greenhouse under strictly nonflowering 
(noninductive) conditions. Cold-requiring plants were kept at temperatures 
above 17°C., long-day plants were kept on a 9-hour day (light, 8:00 a.m. to 5:00 P.m.), 
and short-day plants were kept on photoperiods of at least 19 hours (supplementary 
light of approximately 30 foot-candles from 5:00 p.m. to 2:00 A.M.). 

Prior to the start of the gibberellin treatment, the plants were allowed to reach 
an age at which they were fully capable of responding to thermo- or photoinduction. 

Two types of gibberellin preparations were used. One, supplied by Dr. F. H. 
Stodola (USDA Northern Utilization Research Branch, Peoria, Ill.), consisted of 


709 





710 BOTANY: A. LANG Proc. N. A. §, 
a mixture of gibberellin A; (= gibberellin A) and gibberellin A; (= gibberellin X = 
gibberellic acid). The other, received from Dr. C. Leben (Lilly Research Labo- 
ratories, Indiznapolis, Ind.), contained only gibberellin As and a small admixture of 
an unknown material which apparently has no gibberellin activity. The prep- 
arations were applied in aqueous solution, to which a wetting agent, Polyglycol-31 
(Dow Chemical Company), had been added at a rate of approximately 0.5 ml. per 
1-liter solution. Controls were treated with water plus wetting agent. Rye was 
treated either by administering the solution (0.1 ml.) into the angle between the 
youngest fully emerged leaf blade and the following leaf sheath or the stem, or by 
cutting the tip of one of the younger leaves and inserting the stump into a small 
vial with gibberellin solution. The second method gave a much greater response 
than the first. In the other species, 0.1 or 0.2 ml. of solution was applied once 
daily to the growing region of the shoot. The treatment was continued either for a 
specified period of time or until flower buds became visible. Plants that showed 


no visible flowering response by the end of the treatment were dissected (usually sev- 
eral weeks later) and examined for presence or absence of microscopic floral 
structures. 


EXPERIMENTAL 


Results with Cold-requiring Plants —The results of gibberellin treatment of bi- 
ennial plants are given in Table 1.2 The experiments were conducted under long- 


TABLE 1 
Errect oF GIBBERELLIN TREATMENT ON FLOWER FoRMATION IN BIENNIAL PLANTS 


Hyoscyamus 
niger 
(R«ENNIAL) 
Purple- Yellow- 


Brassica 
Daucus carota napus 
Early Danver's PURPLE 
(MIcroGRams/ flowered flowered French Half Top sativum 
Piant/ Day) Strain Strain Forcing Long GLOBE Hamburg Plain 


No. plants treated 16 15 14 10 14 5 5 


GIBBERELLIN Petrosilenum 


20 


50 


LOO 


No. plants flowering 
Weeks to anthesis 
No. plants treated 
No. plants flowering 
Weeks to anthesis 
No. plants treated 
No. plants flowering 
Weeks to anthesis 
No. plants treated 
No. plants flowering 
Weeks to anthesis 
No. plants treated 
No. plants flowering 
Weeks to anthesis 
No. plants treated 
No. plants flowering 
Weeks to anthesis 
No. plants treated 
No. plants flowering 
Weeks to anthesis 
No. plants treated 
No. plants flowering 
Weeks to anthesis 


0 0 


28 
25 
8 to 12 
27 
26 


8 to 12 to 8 


_ Sees @ 
60 
7 to9 


0 0 0 0 ( 


9 
0 


ae 

3* 

>27 

14 

14* 

18 to >27 
14 

14 
9 to 25 

10 

10* 
20 to >27 
4 
3” 
21 to >27 
4 


+r 


21 to >27 


10 
0 


10 
0 


ie: 

7* 

>27 

10 

10* 
23 to >27 
9 

g* 

>27 

5 

4* 
26 to >27 
5 


4* 
18 to >27 


* Plants, or part of plants, failed to reach anthesis before end of experiment. 
t Tip injured in plants that failed to respond. 


of 
0 
14 
0 
14 
4* 
>27 
14 


day conditions. Since biennial plants, as a rule, require not only a cold treatment 
in order to form flowers, but also long photoperiods (which must. be given following 
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the cold treatment), only one of the two blocks to flower formation which may be 
created in these plants was present in the experiments. 

All four biennial species responded to gibberellin application with flower forma- 
tion; the controls remained strictly vegetative.. Optimal response was reached 
(except, perhaps, in parsley) with a daily dose of 10 ug. of gibberellin. Higher 
doses seem to be somewhat less effective; this point, however, needs confirmation. 

The degree of response varied considerably, depending on the species. Hyoscya- 
mus responded promptly and uniformly, carrot somewhat more slowly and less 
uniformly, the onset of anthesis extending over many weeks. In both species the 
gibberellin-induced flowering response was quite similar to the cold-induced one 
(Fig. 1). The flowers were normally developed and, at least in Hyoscyamus, fully 


Fie. 1.—Carrot, Early French Forcing. Left: control; right: 8 weeks’ 
cold treatment; center: 10 ug. of gibberellin daily. The flowering plants 
are 1 meter tall. 


fertile. (Seed set in carrot was generally poor in the greenhouse.) In turnip, 
only flower buds had been produced at the end of 6 months of continuous treat- 
ment, while in parsley, only a few plants had formed microscopic inflorescences 


over the same period of time. 
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Petkus winter rye was treated for 3-4 weeks with gibberellin solutions ranging 
from 0.1 to 50 mg/liter. Stem elongation was markedly enhanced, particularly 
by the higher concentrations, the stem becoming thin and the leaves small. How- 
ever, flower formation, as measured by the stage of ear development and by the 
number of leaves produced prior to ear initiation, was not significantly affected. 

Results with Long-Day Plants.—The results with Samolus, Crepis, Silene, and 


D) 


annual Hyoscyamus are presented in Table 2.2. Samolus and Crepis responded 


TABLE 2 
Errect or GIBBERELLIN TREATMENT ON FLOWER FoRMATION IN LONG-Day PLANTS 


Gibber- 
ellin Hyoscyamus 
(Micrograms Samolus Crepis Silene niger 
Plant /Day) parvifloris tectorum armeria (Annual) 
0 No. plants treated 32 9 12 25 
No. plants flowering [* 0 0 
Weeks to anthesis 8 cca 
No. plants treated 29 8 
No. plants flowering 29 e Ss 
Weeks to anthesis 4to 10 4to7 
No. plants treated 29 6 
No. plants flowering 29 
Weeks to anthesis 4to7 
No. plants treated 2 
No. plants flowering 
Weeks to anthesis 
No. plants treated 
No. plants flowering lt 
Weeks to anthesis are 
20 No. plants treated 6 
No. plants flowering 3t 


Weeks to anthesis 9 to >20 
6 


50 No. plants treated reg 
No. plants flowering ; oo 6 
Weeks to anthesis ree ae 12 to 15 
100 No. plants treated 6 BED 6 
No. plants flowering 6 ee 6 
Weeks to anthesis 4 ays 14 to 18 


* One plant flowered spontaneously for unknown reasons. : 
+ Plants, or part of plants, failed to reach anthesis before end of experiment. 


rapidly with flower formation, and the number of flowers or inflorescences and their 
development were similar to those in plants transferred to continuous long day 
(see Fig. 2). However, the seed set in Samolus was definitely below that of long- 
day-induced flowers. (The seed set of Crepis, like that of carrot, is usually poor 
under greenhouse conditions. However, gibberellin-induced inflorescences did 
form some viable achenes.) The optimum of the flower-inducing effect is reached 
at about 10 ug. of gibberellin daily. Higher doses result only in a slightly earlier 
response. 

Silene took a longer period of time and larger doses of gibberellin to initiate 
flowers, the plants first making a considerable amount of stem growth (Fig. 3). 
However, those receiving 50 and 100 ug. of gibberellin daily all came ultimately 
to flowering. The inflorescences were much less compact than in long-day-induced 
plants, but the flowers were well developed and produced at least some viable 
seed. In annual Hyoscyamus, only part of the plants initiated flowers; the flower 
buds were abortive. Daily treatment with 10 ug. of gibberellin or more sometimes 
‘aused injury and death of the tip of the elongating stem, accounting in part for the 
incomplete response. However, it is evident that, as in the biennial plants, the 





Fig. 2.—Samolus parviflorus. Left to right: controls; plants treated with 1, 2, 5, 
10, and 20 yg. respectively, of gibberellin daily. The tallest plants are 40 cm. in 
height. 


Fic. 3.—Silene armeria. Left to right: controls; plants treated with 2, 5, 10, 
20, and 50 yug., respectively, of gibberellin daily. The tallest plant on the right 
is 1.5 meters high. 
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degree of flowering response which can be obtained by gibberellin treatment of 
long-day plants varies considerably, depending on the species. 

Short-day-grown spring rye responded to gibberellin application similarly to 
non-cold-treated winter rye, that is, stem elongation was promoted but flower 
formation was neither promoted nor inhibited. 

Experiments with Short-Day Plants—Soybean and cocklebur were treated with 
doses of up to 10 ug. of gibberellin daily. The periods of treatment were 4 and 6 
weeks, respectively. Stem elongation was very pronounced, the soybean plants 
assuming a twining habit of growth. However, the plants remained vegetative. 

Specificity of the Gibberellin Effect—In order to determine the specificity of the 
effect of gibberellin on flower formation in biennial and long-day plants, Samolus, 
Silene, and annual and biennial Hyoscyamus were treated with auxin (indoleacetic 
acid and naphthaleneacetic acid) and with kinetin, separately and in mixture. 
The maximal amounts (20 ug. of auxin, 2 wg. of kinetin) and the periods of applica- 
tion (4 weeks or more) were such that a similar treatment with gibberellin would 
have caused at least a pronounced bolting response. However, the plants con- 
tinued to grow in the habit of rosettes and remained vegetative. The only visible 
effect was some epinasty of the petioles of the younger leaves in plants receiving 
relatively high auxin dosages, 


DISCUSSION 


Flower formation was induced in all four biennial plants and in four out of five 
long-day plants that were treated with gibberellin. No flowering response was 


obtained in the two short-day species that were tested, and flower formation was 
not affected in winter and spring rye. Positive results have also been observed in 
cabbage, rutabaga, and a number of other biennial plants;* in endive, a winter- 
annual plant;* and in Crepis leontodoides,>’ Lapsana communis,’ dill, spinach, radish, 
and several other long-day plants.* On the other hand, no induction of flower 
formation was found in the short-day plant Kalanchoé blossfeldiana.?’ Of particular 
interest is the finding of Biinsow and Harder® in Bryophyllum crenatum and B. 
daigremontianum. These two species are long-short-day plants,’ and gibberellin 
treatment enabled the plants to bypass the long-day part of photoinduction but 
apparently not its short-day part. It thus appears that application of gibberellin 
allows numerous plants to overcome cold and long-day requirements in flower 
formation but that it does not substitute for any short-day requirement. 

These findings amplify our understanding of the role of gibberellin in growth and 
development of plants, and they provide us with a tool for the biochemical control 
of flower formation. 

Gibberellin was discovered as a promoter of shoot elongation.'® In the course 
of the last year, however, quite a variety of other gibberellin actions have been 
found: the effect in flower formation and effects on cell division,'! seed germina 
tion,'® fruit set,'® inhibition and dormancy phenomena,'* and others. On the 
other hand, the existence of gibberellin-like substances in flowering plants has 
been established by Phinney, West, Ritzel, and Neely" and others.’6 There seems 
little doubt that we have before us a new group of native plant regulators with a 
spectrum of functions which, quantitatively, is not less wide than that of auxin 
but is quite distinct from the latter qualitatively. The flower-inducing action of 
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gibberellin was not duplicated by treatment with auxin, kinetin, or mixtures of the 
two and seems therefore to be rather specific. 

Induction or promotion of flower formation by treatment with defined chemical 
compounds has so far been accomplished only under highly special conditions or in 
isolated instances. Photoinduction in long-day plants can be promoted by treating 
the plants with sugars” or auxin,'* but the effect is present only under photoperiodic 
threshold conditions and not in strictly noninductive daylengths. Auxin also 
promotes flower formation in the pineapple,'® the litchi,” and certain varieties of 
sweet potato,?! while maleic hydrazide was found to induce flower formation in 
young celery plants.2*? However, in numerous other plants the same compounds are 
either inactive with regard to flower formation or have an inhibitory effect on this 
process. The gibberellins are the first substances that have been found capable of 
inducing flower formation under unequivocally nonflowering conditions and in at 
least numerous representatives of large physiological groups of plants, the cold- 
requiring and the long-day plants. In view of the presence of gibberellin-like 
materials in higher plants," ‘© it seems reasonable to assume that such materials 
play an important part in the processes that lead to flower formation. 

The precise function of the gibberellins in flower formation will have to be the 
subject matter of further study. However, several comments appear in order at 
the present moment. While the majority of cold-requiring and long-day plants 
that have been treated with gibberellin have responded with flower formation, some 
so far have not. This is true for rye, and Lona*: ** lists several other species (all 
long-day plants). In those species that did respond, the degree of response was 
variable, ranging from one similar to the effect of optima] thermo- or photoinduction 
to very late and weak flowering. This variation is not associated with differences 
at least in the photoperiodic sensitivity of the various species, for Samolus, Crepis, 
Silene, and annual Hyoscyamus do not differ appreciably in this respect. Two 
possible explanations come to mind. (a) The work of Phinney et al. shows that 
the gibberellin-like materials from higher plants are not chemically identical with 
the known gibberellins, which are products of certain strains of the fungus Fusarium 
moniliforme. The evidence suggests the existence of a variety of different gibberel- 
lin-like substances in higher plants, and it is conceivable that these different gib- 
berellins have some degree of species-specificity. In that case, a given plant may 
not respond optimally to a “foreign” gibberellin. (b) Gibberellin may not be 
the only factor limiting flower formation in cold-requiring and long-day plants 
but may be acting in conjunction with some other factor, for example, an inhibitory 
material. In such an event, the amount or activity of this second factor will de- 
termine the effectiveness of added gibberellin in flower formation. One preliminary 
result is in agreement with this assumption. In a small experiment with carrot 
and turnip it was found that the flowering response was more rapid and uniform 
when the temperature was relatively low, although not low enough to effect direct 
thermoinduction. A similar experience is reported by Wittwer and Bukovace.* 
It has been suggested that the effect of low temperature on flower formation in 
cold-requiring plants consists in changing the balance between promotive and 
inhibitory processes, low temperature shifting this balance in favor of the former.*4 
If gibberellin is part of the “‘promotive processes,” it should promote flower forma- 
tion most actively at comparatively low temperatures. The reason for the failure 
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of some plants to form flowers upon gibberellin treatment may have been that 
these plants were treated under conditions in which the balance of the promotive 
and inhibitory processes of flower formation was too far on the side of inhibition. 

Another question that will require careful study is whether the effect of gibberellin 
in flower formation is direct or secondary. Both thermo- and photoinduction re- 
sult in the appearance of a transmissible floral stimulus in the plants. Attempts 
at isolating this material have failed, but results of grafting experiments suggest 
that the floral stimuli of biennial and of photoperiodically sensitive plants are 
different, while those of long- and short-day plants are interchangeable and thus, 
presumably, identical.*4 If this conclusion is valid, gibberellin, which causes flower 
formation in long- but not in short-day plants, is not the floral stimulus itself. It 
may, however, be a factor involved in the formation of this stimulus in long-day 
(and cold-requiring) plants. This interpretation is in agreement with a difference 
in the flower-inducing action of low temperature and long days, on the one hand, 
and of gibberellin, on the other. Most cold- and long-day-requiring plants, when 
grown under noninductive conditions, have a rosette habit of growth. All bien- 
nial and long-day species used in the present work belong to this type. Induction 
in these plants causes not only flower formation but also stem growth (bolting), 
and the two processes occur usually almost simultaneously. If the same plants 
are treated with gibberellin, stem elongation in most cases precedes flower initiation 
in a conspicuous manner. It appears from this that the primary effect of gibberel- 
lin is on stem growth and that flower formation is induced indirectly, the bolting 
plant becoming capable of forming floral stimulus. 


SUMMARY 


Treatment with gibberellin induced flower formation in four non-cold-treated 
biennial plants (biennial Hyoscyamus niger, carrot, turnip, parsley) and four long- 
day plants (Samolus parviflorus, Crepis tectorum, Silene armeria, annual Hyoscyamus 
niger) kept on short day. The degree of the flowering response varied from one 
similar to the effect of optimal thermo- or photoinduction to a late and incomplete 
response. Gibberellin treatment did not affect flower formation in unvernalized 
winter rye and short-day-grown spring rye and did not induce flower formation in 
two short-day plants (Biloxi soybean and cocklebur) maintained on long day. It 
is concluded that gibberellin-like materials play a part in the processes of flower 
formation but that they probably function in conjunction with other factors and 
that their effect on flower formation may be a secondary one. 


The author is greatly indebted to Mrs. Patricia F. Finn and Miss Joan Roach 
for their unfailing and conscientious help in this work and to Drs. F. H. Stodola 
and C. Leben for the generous supply of gibberellin. 


* Research supported in parts by grants from the National Institutes of Health, United States 
Public Health Service (RG-3939), the Lilly Research Laboratories, and the National Science 
Foundation (G-3388). 
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CHANGES IN PROTEIN TOPOGRAPHY UPON OXYGENATION 
By Irvine M. Kiorz anp Ricuarp E. HEINEY 
DEPARTMENT OF CHEMISTRY, NORTHWESTERN UNIVERSITY, EVANSTON, ILLINOIS 
Communicated by Theodore Shedlovsky, June 11, 1957 


The site of attachment of oxygen in the oxygen-carrying pigments has long been 
specifically identified as the metal in these proteins, the metal being attached to 
porphyrin groups in hemoglobin and chlorocruorin or directly to side chains of 
certain amino acid residues in hemerythrin and hemocyanin. It seemed of interest 
to inquire further whether the bulk of the macromolecule participates in the process 
of oxygen uptake. Recent work! has indicated that changes in optical rotation 
of proteins reflect modifications in over-all configuration. Measurements of spe- 
cific rotation would thus seem in order. 

A few measurements for hemoglobin are available in the older literature.” § 
The significance of the values observed [ares mu Cv +15°] is complicated by the 
presence of the porphyrin prosthetic group and of co-operative interactions between 
sites. Consequently it seemed appropriate to consider first the behavior of a pig- 
ment in which these complications do not occur. An examination was therefore 
made of the variations in optical rotation upon deoxygenation and oxygenation 
of hemocyanin, in which oxygen is held by copper attached directly to the protein 
fabric. In samples of this protein obtained from Busycon, the active sites act 
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independently of each other in isoionic solutions; the oxygenation curve fits a Hill 
equation® ' withn = 1. 

Since solutions of hemocyanin are deeply colored, measurements of optical rota- 
tion were made with a Bellingham and Stanley polarimeter to which a photomulti- 
plier detector had been attached. Furthermore, in view of the strong absorption 
of yellow light by this protein, a mercury-cadmium lamp, rather than a sodium lamp, 
was used, with a strong line at 4358 A, near the minimum in the absorption curve 
of the protein. 

Some observations are summarized in Table 1. Oxyhemocyanin from Busycon 
canaliculatum showed some small variations in specific rotation (—56° to —62°) 
from sample to sample. A sample with a = +61.8° was treated with NaHSOs, in 
an arhount approximately five times the copper content of the protein solution. 
Within 15 minutes a dropped to —34.5°. Concurrently the blue color, and the 
associated absorption peak at 580 mu, disappeared. Likewise, the much more 
intense absorption band of the oxygenated protein at 340 mu vanished completely. 
When the solution was then allowed to stand open to air, the characteristic visible 
blue and the near-ultraviolet peak at 340 my of the oxygenated protein returned, 
though not to its full intensity, and a@ rose again to —59.2°. | 

Similar results were obtained with less NaHSOs;, in an amount about twice the 
copper content (Table 1). 

TABLE 1 
Sreciric RorTatTions OF DEOXYGENATED HEMOCYANIN 
AND OXYHEMOCYANIN AT 25°C, 
Solution Specific Rotation, a 
Oxyhemocyanin; dialyzed sera —59.0°, —56.2°, —57.7°, —61.8° 
(approx. 1 per cent protein) 
Deoxygenated hemocyanin; Solution 1 + NaHSO;; —34.5° 


NaHSO,;/Cu & 5 
Solution 2, exposed to air for 40 hours —59.2° 
Oxyhemocyanin : —62.3° 
Deoxygenated hemocyanin; Solution 4 + NaHSO;; —37.7° 
NaHSO;/Cu = 2 


Oxyhemocyanin : —59.5° 
Deoxygenated hemocyanin; Solution 6 + glucose oxi- —21.8° 
dase-catalase (0.003 per cent) and glucose (0.018 per cent) 


As a further check, it seemed desirable to avoid a method of abstracting oxygen 
from the protein by reduction and to use a method of removal which depleted the 
dissolved oxygen in solution in equilibrium with the hemocyanin. Despite in- 
tensive efforts, complete displacement of dissolved oxygen by means of bubbled 
high-purity nitrogen could not be attained; although the weak 580 my peak of 
hemocyanin disappeared, it did not prove possible to reduce the intense 340 mu 
absorption to the background value. An enzymatic method'! was therefore used. 
A commercial preparation of glucose oxidase and catalase at a concentration of 
0.003 per cent in the presence of 0.02 per cent glucose completely removed the oxy- 
gen from solution, as judged by the disappearance of the 340 mu band of he- 
mocyanin. As in the bisulfite treatments, this deoxygenated hemocyanin 
showed a pronounced drop in optical rotation; as compared to oxyhemocyanin 
(Table 1). 
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It is apparent, therefore, that in the oxygenation of hemocyanin, in which an 
O2 molecule is bound at a specific localized site, the protein fabric is affected in a 
very extensive fashion. Such a possibility has also been suggested for hemoglobin, 
by Haurowitz'? to account for crystallographic differences between ferro- and oxy- 
hemoglobin, and by Wyman and Allen'* and by St. George and Pauling" as an 
explanation of the co-operative interactions in the uptake of successive O2 mole- 
cules. From the present experiments it is clear that substantial changes in topo- 
graphical or structural configurations occur even in the absence of porphyrin 
prosthetic groups and even when the binding sites do not show any co-operative 
interactions. The very act of taking up an oxygen molecule makes itself felt over 
wide regions of the macromolecule. 

The nature of these configurational changes is still obscure. The general trend 
has been to attribute them to steric displacements or rearrangements in polypep- 
tide chains. It is now clear thet these steric effects, if they exist, do not depend 
on the presence of a large porphyrin moiety. It is perhaps pertinent to point out 
also that the direction and magnitude of the change of optical rotation upon oxy- 
genation are similar to those observed on denaturation of proteins and hence would 
imply a substantial increase in disorganization of the macromolecule. Neverthe- 
less, oxygenation-deoxygenation processes can be carried through many cycles 
without appreciable change in the properties of the protein. 

We should like to consider, therefore, an alternative model which attributes 
changes during oxygenation to modifications in the “frozen water’’®~'* of the hy- 
dration lattice of the oxygen-carrying protein. From this point of view, the 
binding site of the deoxygenated protein may be visualized as having an aqueous 
covering sheath whose lattice structure is determined by the nature of the residues 
(and metal) at the active site and by the amino acid side chains in the neighborhood 
of the site. Under appropriate conditions—for example, hemocyanin above its 
isoelectric point, or hemoglobin at any pH—these crystalline water islands may 
be sufficiently extensive that they merge into an interwoven lattice. The attach- 
ment of an oxygen molecule to the binding site would produce an incongruity in 
this lattice, because of the larger size'* of O. as compared to H.O and because of 
the differences in hydrogen-bonding properties. As a result, the hydration “‘ice- 
berg” might be disturbed for an appreciable distance around the active site. This 
disturbance could affect the hydration structure of neighboring sites and hence 
change their affinity toward oxygen molecules. Simultaneously, it would manifest 
itself by changes in optical rotation. Changes in degree of interaction between 
sites of hemocyanin due to changes in pH may be attributed to changes in hydration 
as the charge of the protein is varied. 

In terms of changes in the structure of the hydration lattice, one could also 
understand why so many different substances affect the affinity, and hence the 
shape and position of oxygenation curves, of the oxygen-carrying pigments. For 
hemoglobin these substances may vary from presumably innocuous salts?® which 
affect only the ionic strength, to specific reagents such as mercurials*! which com- 
bine with mercaptan groups, to generally reactive substances such as formalde- 
hyde.*?. All these substances can increase the oxygen affinity of hemoglobin, and 
they affect the degree of interaction between sites; all these substances would also 
modify the lattice structure of the aqueous hydration sheath. 
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From this point of view the observed dielectric behavior of hemoglobin also 
seems reasonable. It has been reported by Takashima”* that the high-frequency 
dielectric increment of hemoglobin changes with degree of oxygenation. As has 
been emphasized by Jacobson,‘ dielectric behavior reflects the character of the 
hydration sheath of a macromolecule. It is only reasonable, then, that binding of 
oxygen should change the dielectric properties if it affects the water lattice of 
the protein molecule. 

It would seem, therefore, that the character of the “frozen’”’ water on the surface 
of the protein may play an important role in determining the behavior of the oxygen- 
carrying pigments. 


This investigation was supported in part by a research grant (RG-4134) from 
the National Institutes of Health, United States Public Health Service. 
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VECTOR MODEL FOR INDIRECT PROTON HYPERFINE 
INTERACTIONS IN xr-ELECTRON RADICALS* 


By Harpen M. McConne.i 


GATES AND CRELLIN LABORATORIES OF CHEMISTRY, f 
CALIFORNIA INSTITUTE OF TECHNOLOGY, PASADENA 4, CALIFORNIA 


Communicated by Richard M. Badger, June 24, 1957 


Significant information on unpaired electron distributions in molecules, liquids, 
and solids can sometimes be obtained from magnetic resonance studies of electron- 
nucleus hyperfine interactions. A particularly important problem in this field 
involves the proton hyperfine interactions that have been observed in the electron 
magnetic resonance spectra of a large number of aromatic radicals in solution.'~* 
The problem concerns the theoretical basis for a proposal that the observed hyper- 
fine splitting ay due to aromatic proton N in an aromatic radical can be used to 
estimate the unpaired spin density py at carbon atom N to which proton N is at- 
tached, according to the equation” *~* 


ay = Q py. (1) 


A semiempirical value of Q derived from observed splittings ay and theoretically 
calculated spin densities py is —30 + 5 gauss, or —85 + 15 Mc.? Equation (1) 
has been derived recently using molecular orbital theory with arbitrary xz con- 
figuration interaction. We show here how the simple proportionality in equa- 
tion (1) may be derived using the Dirac vector model. 

Let S,, S,, S, be the spin vector operators for electrons A, 6, x. Electron hisa 
ls hydrogen atomic orbital centered on aromatic proton N. Electron é¢ is in the 
sp? hybrid orbital centered on carbon atom N and directed toward aromatic proton 
N. Electron z is in a 2p,7 orbital centered on carbon atom N. The nodal plane of 
the 2p, orbital passes through proton N. Let |n°> be a basic set of spin eigen- 
functions including both o and 7 electrons. Assume the Hamiltonian matrix 
H,,,, to be diagonal with respect to the set |n°> except for nondiagonal terms due 
to é—z exchange coupling. Treat this é~x interaction as a first-order perturbation 
on the ground-state spin function, |0°>: 


xc = const. — '/.(1 + 48,-S,)J ez. (2) 


In equation (2) J,, is the exchange integral between electrons 6 and x. To first 


order, the exact ground-state spin function |0> is 


<n° |S,+S, |0°> 
u AE’ - in°>, (3) 


l0> = |0°> + 2/., 


where AE, = E, — Eo. We use |0> in equation (3) to calculate the hyperfine 
splitting due to proton N, which is 


> <0|S;- |0> 
an = ay g ‘ 


(4) 


z 


721 
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Here, in standard notation,® 6 


Ds ee 8g I | Mu 
“4 a Te 


1,420 Me 


S, = <0|> S,,|0>. (6) 
k 


The sum in equation (6) is over all spins in the molecule. From equations (3) and 
(4) we obtain, to first order, 
cc Ad oe 17 _ <0°|S,-S, In°> <n° Snz|0°> 


a as AE, 


ay (7) 
Explicit calculations show that only a very limited number of o—7 excited spin 
states can contribute to the summation in equation (7), and if we assume that all 
of these o—z excited states have roughly the same excitation energy, AF, then 
equation (7) can be simplified to give 

4] an 

ay = —“— <0°|S,-S,S,,|0°>. (8) 

AE § 2 
To evaluate the above diagonal matrix element, we note that if \1I°> is the nor- 
malized pure z-electron spin component of |0°>, then the « component contains 
the factor (a(é) B(h) — B(6) a(h))/+/2, and when the integrations over the spin 
variables of A and o (and the other o electrons) are carried out, we obtain 


_ fda\ , « <tt*|8,,[i°> 
all AE] ” S, ; 


(9) 
Since py = <II° [Ses 11°>/S,,° we obtain equation (1) from equation (9) whenQ = 
—J,, 4y°/AE. Previous calculations indicate J,,/AE to be of the order of mag- 
nitude of 0.1, so that this theoretical value of Q is in order-of-magnitude agreement 
with the semiempirical value.® 

The critical step in the present derivation of equation (1) is the assumption of 
an approximately constant AE ~ AE, for all the interacting states in°> in equa- 
tion (7). A sufficient condition for the validity of this assumption is that all 
mn electron exchange integrals with magnitude Fave are very small compared 
with the CH bond integral \Jon I. In this case AE, is determined by J,, alone. 
It is significant that the conditions we find here for the validity of equation (1) 
using the vector model are compatible with the molecular orbital derivation of 
equation (1), with arbitrary 2-2 configuration interaction. In both cases the 
single proportionality between spin density and proton hyperfine splitting in equa- 
tion (1) becomes exact as i Jes | — 0 and ii Jen| — (0. Detailed cal- 
culations, including nonorthogonality effects, are of course required to test the 
adequacy of these limiting approximations. 

* Supported in part by the Office of Ordnance Research, U.S. Army. 

+ Coniribution No. 2214. 
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CRITIQUE OF THE THEORY OF OPTICAL ACTIVITY 
OF HELICAL POLYMERS 
By Wiiui1am Morrirr 
DEPARTMENT OF CHEMISTRY, HARVARD UNIVERSITY 
AND 
Donap D. Firrs* anp Joun G. Ktrkwoop 


STERLING CHEMISTRY LABORATORY, f YALE UNIVERSITY 


Communicated June 10, 1957 


1. Introduction.—In a series of recent articles,'~* we have severally attempted 
to use coupled-oscillator theory in order to understand the optical activity of poly- 
peptides. This approach, the development of which stems from Born® and Kuhn, 
was first used in quantum mechanics by Kirkwood’ as the basis for his polarizability 
theory of rotatory power. In recent months we have been clarifying the relation- 
ship between our respective treatments, leaving out of account differences which 


stem from the specific means by which we have attained tractability. As a result, 
we have been able to show that a general formula, derived by Kirkwood’ under 
special (nondegenerate) conditions which do not apply to polymers, has a wider 
range of validity and may also be used in the present context. But, paradoxically 
enough, we have thereby found an error in Moffitt’s work,* one of the purposes of 
which was to acknowledge explicitly the degenerate nature of this problem. 

In the following section, Kirkwood’s formula is rederived under conditions apper- 
taining to a sequence of identical residues. And in a final section this prescription 
is more specifically applied to helices. Some discussion is given of the way in which 
the result may be used. It appears that a particularly important factor (A,,) is 
composed of two terms, one of which (A,,‘) was correctly assessed by Moffitt, 
whereas the other (A,,‘”) was absent from his earlier derivation. The reason for 
this omission is analyzed and shown to originate, somewhat unexpectedly, in the 
failure of the Born-von Karman boundary conditions which he had used in the 
course of his exciton treatment. 

We are forced to conclude that the close agreement between the data presented 
by Moffitt and Yang,* and Moffitt’s calculations for the a-helix,‘ is illusory and 
does not enable a definite screw sense to be assigned to these polypeptides. Moffitt’s 
calculations were based on certain A,,‘"’ terms alone, and our estimates indicate 
that the corresponding A,,” contributions are of the same order but of opposite 
sign. However, since the qualitative form of the dispersion shown by the various 
partial rotations was correctly described in his work, this quantitative failure in 
no way vitiates the empirical correlations put forward by Moffitt and Yang; any 
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reservations we may have on this score have already been enumerated by these 
authors. 

2. Degenerate Coupled-Oscillator Theory.—In this section we shall rederive Kirk- 
wood’s formula in such a way as to include, explicitly, the case of identical residues. 
It is convenient to use a modified form of the notation employed by Moffitt;’ thus 
we adopt Section 2 of that paper, with only a change of phases (see below), and 
those parts of Section 3 which precede the paragraph containing equation (24). 
We shall avoid specifying the particular geometry of the polymer at this stage, so 
that our results are equally valid for random coils and for helices. 

Our polymer therefore consists of an array of N identical residues, each of which 
is identified by a serial index m. The ground and excited states of the mth residue 
are represented by the real functions y,» eal Vom, respectively.’ It is supposed 
that the interaction between different residues is sufficiently weak—or that the 
symmetry of the array is sufficiently high—that the respective excitation energies 
Eom = €, are the same for all residues. 

The ground state of the polymer as a whole is represented, in zeroth-order, by the 
simple product IT yYo,,(7,,), which we call '0); tm Symbolizes the co-ordinates of the 

m 
electrons associated with the mth residue. At first, its excited states are also repre- 
sented by simple products; if only one residue is excited (‘‘fundamentals’’), 


| em) _ 0) ( (Wom) ~ ous (1) 


if two or more residues are excited (‘‘overtones’” and “‘combinations’’), 


| om7n) sa 0) (WomVon) ~' Vom¥rn (2) 
= | rnom), (n # m); 


and so on. Owing to the possibility of synchronizing the instantaneous charge 
distributions in different residues, these functions must be combined linearly to 
give the correct zeroth-order functions, namely, 
|oK) = > |om) (om| cK), 
m (3) 
lorK) = Y% DY |omm) (omrn| orl), , 
m n#m 
where the concomitant transformations are all unitary and serve to diagonalize 
the Hamiltonian locally: typically, 
( qian, "= KE KoKK’: (4) 
The stationary states of the polymer, correct to first-order terms, are now repre- 
sented by 
0) = |0) + 222) |erV) (prU'\0) 
\oK) = | cK) -_ > |e) (rLi\oK) + Ld | ere) (orpW '|oK), 
THE 
with similar expressions for the upper states of combinations bands. Here 
(prU||0) = —(prU!3 0)/(e, + «€ a 
(rL||oK) = —(rL|X| cK) /(e, -— €, 
(orpW||\ cK) = —(orpW \5¢| cK) )/(€p = ie), 
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and it has been assumed that the numerators are first-order small quantities with 
respect to the denominators. 

We now compute, correct to first order, the matrix elements of a typical ‘“one- 
electron’’ operator @, such as the linear or angular momentum, corresponding to 
the excitation of afundamental. Setting 


V or | 


mn = (rm | om), 
Og = (0| | om) = (om| | 0)*, 


om 


(7) 
we find, after a little rearrangement, that 


(0) @| ok) a a {Pom ad Zz Fue Say” 2€. 


m n=zm 


_ x py Vinn’ (Prn*/(€, + €¢) + Prn/(€, — e,) |} (om| cK). (8) 


THe n#~m 


With our phase choice, we may enumerate two particular forms which this expres~ 

sion may take, depending on whether @,,, is real or imaginary. If @ is to represent 

the total linear momentum P = >> #,, (or, again, the total angular momentum 
m 


L), Pom = Pom is imaginary and therefore 
{0| P| ok) rT > \ Pom = ‘2 at aa 2€, 
m nxm 
2 > th V an Dente/(€,7 — €.”)} (om/cK). (9) 


THO nem 


On the other hand, taking 0 = Q = >> gq», where q,, is the sum of the position 


m 


vectors of all the electrons in residue m, we see that qg,,, is real, so that 


(0| Q| cK) = do dem — DLS GenV mn”’/2€e 


m nzm 


-2> > Vue Geiiel i = e,2)} (om/cK). (10) 


TH? nz=m 
Since Wom, Vem are the eigenfunctions of a self-consistent field operator 3,,, 
Pen _ (m, th) €oQom; (11) 


where m, is the electronic mass. Using this relation and the properties derived 
above, it is now easily verified that, correct to the first order, 


(0| P| cK) = (m,/ih)E,x (0| Q| cK), (12) 


where, for simplicity, we have taken as origin for the energy that of the ground state 
in zeroth-order. Thus the internal consistency of our model allows us to use with 
impunity either the dipole moment. or the dipole velocity expressions relating to the 
optical properties of the polymer. 

It is a straightforward matter at this stage to write down, correct to first-order, 
the rotational strength associated with any transition in the polymer. Since this 
quantity is the scalar product of two vector matrix elements, namely, of Q and of 
L, it is clear that, to our approximation, at least one of these elements must be non- 
vanishing in zeroth-order. Accordingly, we may ignore the contributions from 
overtones and combinations, confining ourselves to the fundamentals whose rota- 
tional strengths are typified by 





CHEMISTRY: MOFFITT ET AL. Proc. N. A.S. 


R,x = (e2/2m.c)Im[O| Q| oK)+(oK|L|0)]. (13) 
Since L is referred to a fixed origin, we write L = >> k,,, preferring to reserve the 
symbol J,, for the angular momentum of the eontens of the mth residue, referred 
to an origin in that residue; we note that 

kom = (0|L| om), (14) 
like p,m, isa pure imaginary vector. Accordingly, 
Rox = (te? ame) u (om| cK) (cK | om’) 
X ifem - Ron + "(2e,) 1 2G Ben Van Ses Ret han | 


nZ=m’ nzm 


¢ 9 a r r | f 
= 2 z: (<, e [dX Jam - es €, J ~ + z. Grn - Ram’ €, J a it (15) 


T#o n=m’ n=zm 


The partial rotation, per residue, due to this transition is therefore 


heF,.\ Rex 
A eg eel SS in. Ty 


e?N (¥eK" ye v*) 


where v,x = E,«/hc is expressed in wave numbers and we have set 


_ A8Noe? (" + *) (17 
he? 3 ’ ‘) 


No being Avogadro’s number and n, being the refractive index of the medium at 
frequence v 

If the frequency of the anger light is not too close to v, = e€,/hce, so that, for 
all K, |v, — vo} > | vex — ve, then the contribution to the residue rotation from 
transitions to all states K in the exciton band Yo — y, is proportional to 

y? » Rox 2y2 Ye 

ae Ro x(EoxK — €g)/ Eq. 
io. 


The summation in the first term reduces to 


‘ . 9 _ y 
(2e?, ‘2m.c) 7? 1WoK* Rom — 2 me :® (e,” — €,) ; hue. 


m rHo nXm 
Xx (€Qem’ Rn + €,9 rm’ Beat 3 


and, dropping second-order terms, the summation in the second becomes 


(1e? '2m,c)e St ® > Gom'R om’ Pp» (om| 5) q KR — éo| ok) (cK | am’). 
K 


m m’ 


Since & is locally diagonal, the last line is just 


| ss oe | - , f a 
dX Dd (om!) cK) ‘ | 0K’) (aK"| om’) = (om| 5 — e,| om’) = Vian" 
K K’ 


which vanishes ifm = m’. Accordingly, the summation in the second term becomes 


pe | lé . == 
(se*/Qmec) DD, cc Vii hn Rik 


m nem 


As a last transformation, we now put 
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Rom = lem + Im X Pom = lem + (Meée/th)Im X Tem (18) 


where 1,,, is the appropriate vector matrix element of the angular momentum re- 
ferred to a point r,,, fixed in the mth residue. 
The partial rotation [m,] may therefore be written as follows: 


[m, | = [m,i| + [m,;] ai [mon] + [mon | “fr [m,;"], 


where 


thF y? 
[mi] Dx ( Dor gtd y?) pa SO om) 


2m,N “ 
F, €,€ 5 
[m,;| “7 “Gy re Zz, :% 9 g J Wiggs a Ta) * (Gem Xx G-n); 
(v," m n#¥m r¥a (€,, — €,°) 
F, yr oo 
[mon] Bo G ool 2D J mn Ting oa Tn)°\Gom x Gen); (20) 
m nxm 


3 
AF 
[mon] na - (= ¥ 2 ie a ‘so p> €a V mn ‘“< ae 


m,N = m n#¥m 


hF, 
[m,;"] ey --(§ ey > p > €,”) Feat 


m,.N m n#m tr¥oe 


x (SMouiitia + €Aca'low) 


[m,;] represents the intrinsic partial rotations of the individual residues, whose 
electrons are subject to slightly asymmetric Hartree fields, arising from any asym- 
metric centers they may possess or from their polymeric environment. If the 
transition ~p —> yw, is associated with a large change in magnetic moment, and there- 
fore a small dipole moment, '® it is said to have a high dissymmetry factor;'! for such 
bands, [m,,;| dominates the partial rotation, and the remaining terms in [m,] may 
safely be dropped. [m,;], which shows the same dispersive behavior as [m,;], arises 
from the first-order corrections, {| ) — |)}, due to the correlative interaction be- 
tween states in different exciton bands. For strong electronic bands, this term can 
be large; in general, it is of the same order as [m,,;|, though, in coupled-oscillator 
theory, one is inclined to hope that it outweighs [m,;|._ The third contribution, 
namely, [m,,]|, shows different dispersion and is a unique feature of the degenerate 
case; it arises from the correlative interaction between different residues, within 
the same exciton band. It was at this term that Moffitt attempted to direct 
attention. The last two terms [m,,] and [m,,-], are small with respect to their 
imprimed counterparts, namely, [m,,] and [m,;], by a factor d,/|1, — T,|, where 
d, is the distance of r,, from the partial optic axis of the transition Yo, > Wem." 
Following Kirkwood,’ we therefore drop both of these. 

It may be noticed that, on summing over all exciton bands, the residue rotation 
may now be written 


[m] = [m,] = uo [m,i] 
y Pe pp rs as Tm)*(Gem X Gen), (21) 


m n 


where, by combining the [m,;| and [m,,]| contributions, we have removed the re- 
striction r # o from the summation over 7. Since vy, = vom, ¥; = Yen forall m,n, 





728 CHEMISTRY: MOFFITT ET AL. Proc. N. A. 8. 


it will be recognized that this formula is identical with that obtained by Kirkwood.’ 
But, whereas the validity of the original derivation was confined to the nondegener- 
ate case, we see that this restriction has been removed by the present treatment. 
If V»»’" is replaced by an appropriate dipole-dipole approximation, the more par- 
ticular results of the polarizability theory follow immediately. 

3. Application to Helices.—Let us now specialize our discussion and suppose that 
the successive residues are arranged to form a right-handed helix. The monomers 
may then be located by the equivalent points r,, in each, namely, by 


Tm = pCOS Gn 1+ p SiN am Jj + mZk, 


where P is the number of residues per turn (and not necessarily integral), PZ is 
the pitch of the helix, and a,, = 2rm/P. Similarly, we may put q,,, = | Ge! Gems 
where é,,,, is a unit vector whose radial, tangential, and vertical components, relative 
to the circumscribed cylinder, are ¢,,, €,;, and é€,,, respectively: 


Com = (Cer COS Am — Co; SIN Aq) + (€,, SIN Ay, + Cyz COB Ay)J + Coyh. (23) 


We suppose the helix is sufficiently long that we may neglect end effects. To 
simplify the discussion, we confine our attention to two transitions, Yo > y, and 
Yo > yo, say, supposing that contributions from these dominate the residue rotation. 
The optical rotatory dispersion is then described by 
j v?yyv2(Ay2 + Aoi) i] v?yy"Ay 
at Lan ao gues y acpi? maar yee 
\ (v4° — p*) (po* ar V;") (14° — p*)° 
v?ven(Ao + Ais) y 
+ 4 [my] + Pn ett 
: ( )) (y 


vo” — v®) (vy? — vy” 


or, equivalently, by 


[m] = [m,,;] + [mao;] 


eet ol i (v2 — v?) (v2? — v?) 


\ eA 4 ee Aa 


where 


Agni F 2a: Ven’ Qn 7 tal nh Cn). 
n(>m) 
In making a theoretical analysis of the “observed” dispersion curves for such an 
idealized system, one may adopt several different viewpoints which it is as well to 
enumerate. 

It should be recognized, first of all, that the calculation of terms involving A,, 
presents a very different problem from that of computing the [m,,]’s. The so- 
called one-electron theory'! is primarily concerned with the latter task, whereas 
the coupled-oscillator theory is best suited to the former. For the present, we 
cannot attempt to assess the intrinsic rotatory parameters and are forced to regard 
these as outside the scope of our work. We are therefore faced with two alterna- 
tives. 

On the one hand, we may hope that the [m,,|’s are appreciably smaller than the 
coupled-oscillator terms. In this case, we attribute the main part of the optical 


rotation to those terms of equation (25) which have been inclosed in braces. We 
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then base our calculations on the polarizability theory or make some other use of 
the optical properties of the individual residues—in particular, of their absorption 
spectra. This approach has been adopted by Fitts and Kirkwood in their work on 
polypeptides." ? 

Alternatively, we may observe that [m,,;] shows the same frequency dependence, 
namely, v?/(»,2 — v?), as the term accompanying it in braces in equation (24); and 
similarly for [mo;]. The other two terms, which involve Ay and Ag, show an un- 
usual dispersive behavior, however. We may therefore hope, on the basis of the 
observed frequency dependence, to distinguish between these two different sets of 
dispersive terms. Adopting this approach, we attempt empirically to separate 
contributions to the rotation due to the terms in braces from those due to the re- 
mainder of equation (24). The latter are amenable to calculation, and we are 
enabled to avoid the [m,,| difficulty without having to neglect it. This avenue 
has been explored by Moffitt and Yang’ in the present context. 

It should perhaps be added that neither method is very satisfactory. Thus 
it may be questioned whether the existing data actually succeed in separating the 
abnormally dispersing terms from their normal counterparts. And, again, it is 
by no means clear that the [m,;|’s make a negligible contribution to the observed 
rotation. However, while making these reservations, we continue to pursue such 
work, since it is unlikely that progress is to be made in any other way. 

Returning to the helical problem, we evaluate A,, by using equations (22) and 
(23). Itis found that 


5. = Ae + a* 


1 + | ‘ , _— , 
) = | el Peeler Ly Vinn*” sin? [x(n — m)/P}, 
n(>m) 


au F,| qe| °Z (Cor? + Cet?) DY, (nm — M)Vnn”” sin [24(n — m)/P]. (27) 
n(>m) 

This expression for A,, should be compared with that derived by Moffitt,’ under the 
same assumptions. (It will be recalled that this term arises from correlative inter- 
actions within the same exciton band and is therefore directly related to the correct 
zeroth-order functions | oK).) He found only the term A,,‘”, so that his analysis 
is in error. Since his mistake is rather interesting, we shall conclude by showing 
how it arose. 

Moffitt considered the helix as a one-dimensional crystal and employed exciton 
theory. By this method, the coefficients (om/oK) of equation (3), which deter- 
mine the correct zeroth-order functions, are evaluated by adopting the Born- 
von Karman cyclic condition, which leads to 
*exp (2rimK/N), (28) 


(om| oK) = N7~' 
where K runs from 1 to N, the total number of residues. The procedure is standard 
in the treatment of the optical properties of molecular crystals, as well as in almost 
all of crystal physics. It is not, however, quite correct, when the crystal is not 
in fact re-entrant. 

It is easy to see that the error incurred in the calculation of (0| Q| eK) by adopt- 
ing this boundary condition is of order N~'’*, although the main part of this term is 
of order N’”* and correctly given. In computing oscillator strengths per residue, 
which are proportional to N-'| (0| Q| cK)| 2, the error is therefore of order N-', 
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which is indeed negligible. However, in the present context, we wish to compute 
rotational strengths per residue, which are proportional to N~'(0| Q| ¢K)| (oK L| 0). 
Now (cK|L|0) involves the moments of the electronic displacements about some 
fixed (but arbitrary) point in the molecule; it therefore contains terms of order 
N’”* as well as of order N'“. (The former were innocuously symbolized by N’'”a,, 
N’’B, in Moffitt’s equations [48],* although, since the difficulty was not foreseen, 
they were not explicitly evaluated.) The scalar product of the N’”’ terms of 
(cK | L|0) taken with the N’”* terms of (0| Q| cK) vanishes identically—as, clearly, 
it must, because the residue rotation is not proportional to the number of residues. 
The predominant N’”* terms of (0| Q| cK) and the N’” terms of (cK L\0) yield 
a contribution to the residue rotation of order unity, namely, the term involving 
A,,‘?. But the scalar product of the N-'” errors of the former taken with the 
N*”* terms of the latter does not vanish and is therefore of order N; this again leads 
to a contribution of order unity to the residue rotation, namely, the neglected term 
A,,.”. It may be noticed that the N’” terms of (oK| L|0) can involve only per- 
pendicularly polarized components of g,, whose cylindrical direction cosines are 
e,, and é,,. These large moments are also perpendicularly directed, so that the 
neglected term cannot involve the vertical component e,,. Looking at equations 
(27), it is therefore clear why A,,‘” was correctly assessed but A,,‘” was missing. 

We may conclude that whereas the application of the Born—von Karman bound- 
ary conditions leads to a correct picture of the spectra of molecular crystals, they 
may not be used with impunity. In particular, they fail to give an accurate 
account of the optical rotatory properties of long helices. 
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Since the early work of Morgan and his students! it has been known that the 
mutant rudimentary (r, 1, 54.5) is sterile in the female sex. Rudimentary females 
are therefore produced by outcrossing and backcrossing with r o@ which are com- 
pletely fertile. Morgan and Lynch proposed an adequate explanation: the 
homozygous rr locus, when acting during oégenesis, produces something in the egg 
cytoplasm which affects adversely the predetermination processes and prevents 
the eggs from developing. In the heterozygote + /r the normal locus counteracts 
this ‘‘poisoning.”” The morphology of this lethal effect was recently studied by 
Counce.” She found in the eggs of some r alleles that development proceeded 
normally up to primary organ formation but stopped completely between 13-16 
hours. Later it might start again, which explains the few rare offspring from rr 
2 already obtained by Morgan. 

Fifteen years ago, while studying a new r mutant called r’””, I obtained strange 
results. While most rr 9 were sterile, some produced offspring. From the mating 
rr X r many normal daughters were obtained, a smaller number of normal sons, 
and a similar small number of r 9 andro. Rudimentary females mated to + ¢@ 
produced many normal daughters and a few normal as well as a few rudimentary 
sons. Normal daughters of rr 2 mated to r o& bred like heterozygotes r/+, and 
normal sons bred as + co’. The first idea was, of course, “contamination.” But 
many repetitions made with all precautions gave the same results, including the 
extreme sex ratio in favor of females. Many experiments were made to find en 


prbl 


hancers or suppressors which might explain the facts, all without result, and the 
work was dropped for the time being. As a matter of fact, Morgan had already 
had a comparable experience. Breeding r 9 with r o’, he had one positive result, 
namely, an offspring consisting of one rudimentary daughter and one rudimentary 
son, but also two long-winged daughters. He assumed the latter to be due to the 
mother’s not having been virgin. He mated the two unexpected females to rudi- 
mentary males and obtained both normal and rudimentary females and males, 
which he considered proof of prefertilization (1915, p. 241).!_ A few pages later (p. 
245) he reports again on offspring of seventy-four rudimentary females. In one 
bottle a few long-winged females and males appeared, “‘which must be due to error 
or contamination.” We shall see that the two results were not necessarily due to 
error or contamination but were typical for some rudimentary lines. 

In Drosophila the cytology of the meiotic divisions in the egg is topographically 
similar to that in Bombyx, where we proved? the presence of a recessive locus which 
induces development involving a polar body. The two second maturation spindles 
are arranged in tandem, and in telophase four nuclei are found vertically to the sur- 
face (see Fig. 1, in which only the X chromosome is drawn). The innermost is the 
egg nucleus and the others the polar nuclei. The second and third nucleus come 
to lie near each other later and are called the “polar copulation nucleus,’’ (PCN) 
which remains for a long time in the egg but is finally absorbed. 


731 
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Now it is possible that in the r eggs of the heterozygous female r/+ fertilized by 
an r sperm, one of the effects of the presence of two r chromosomes is to cause the 
PCN to act as an autogamic zygote nucleus and, further, to move into the polar 
plasm of the developing egg together with the descendants of the proper rr fertili- 
zation nucleus. Thus the developing rr 2 contains in its ovary rr cells, and also 
r/+ cells derived from the PCN of the F; egg. If all rr eggs of these rr 9 were 
sterile, only those eggs could develop which contained the r/+- PCN nuclei, (Fig. 1). 

If this nucleus behaved normally in 
the meiotic division, the offspring of 
such an r/+ egg from an rr mother 
would be females and males, one- 
half rudimentary, one-half normal. 
But there is also the possibility 
either that the PCN from F, re- 
mains together without meiosis or 
that, with correct meiosis, the new 
PCN takes over via autogamy (pseu- 
doparthenogenesis), no fertilization 
occurring. All such eggs will develop 
into r/+ 9. This sounds rather 
unusual. But recently Seiler* found 
in the moth Solenobia triquetrella that 
the parthenogenetic race develops 
from eggs with the PCN furnishing 
the embryonic nuclei. Thus it is to 
be expected that the offspring of ex- 
ceptionally fertile r 9 would consist 
of only r/+ daughters. But if most 

Fic. 1.—The two polar regions of a Droso- of the 7/+ nuclei in the egg of an rr 
pla ca. Abo: nuclews and three polar 9 behave thus, but a few undergo 
somes are drawn instead of nuclei. Left: fer- meiosis and normal fertilization, a 
tization of egg nucleus with by sperm with. few rr @ andr @ and + o? could be 
tion nucleus, PCN, which is r/+. At the op- present. 
at pole lies the polar plasma, determined to The results of this set of experi- 

vecome the germ track. Left arrow: movement 
of descendant of normal fertilization nucleus mentsare tabulatedin Table 1. Only 
nto the ge track gl arm: abnormal 5 of the pair matings rr Xr, but 18 
mass matings (20-30 rr 9), were 
successful. Some such were made with + o (Nos. 1, 2, 4), others with r @ (No.3). 
In No. 4 the rr mothers were not F R, 2 but rr 2 from more complicated crosses, 
which, however, must have had r/+ grandmothers—i.e., were identical in origin 
with the F, rr 2. No.3 shows the typical result of rr X r parents, the three others 
those of rr X + parents. Together there are 1,793 9, 52 o, a sex ratio of 34.5 
2:1. We further see that exactly half the males are + (r). In No. 3 there are 
also a few rr 9, about the same number as of males in one of the male classes. 
This shows conclusively that over 1,700 are the products of pseudoparthenogenetic 
development of the grandparental PCN and that, in less than 3 per cent of the eggs, 
normal meiosis and fertilization had taken place. Only in No. 3 can the daughters 
from fertilization (rr) be seen; in Nos. 1, 2, and 4 they are also r/+. 
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In Nos. 5-8, tests for the + daughters of rr 2 are given. Forty-two such 
daughters were mated to different kinds of males. All results show that the + 
daughters had been r/+. The numbers in the segregating classes point to normal 
meiosis and fertilization and not much pseudoparthenogenesis, if any. 


TABLE 1 
OFFSPRING OF FERTILE r 2 FROM Backcross (N X r) X r AND (FourtH Group) r 
EXTRACTED FROM LATER GENERATIONS BRED FROM THE BACKCROSS; FURTHER 
Tests (5-8) or NorMAL DAUGHTERS FROM rr MOTHERS 


Pair Mass FEMALES 
PARENTS MaTINGs MATINGs + 


eK 5 27 : 

Pi SR Get 7 910 10 
ox se = 
9 


—Mates— 
«fh. 


y r 


6 349 g ‘ 

(later gen.) X o* 5 506 10 9 
+ fromrr Xr, dr 792 495 707 469 
>) + from rr X +, o& brother r é 262 143 191 122 
i 
+ 


en 


oo 


) 


~ 


from rr X r stock j 831 286 597 522 
f ees 


from rr X brothers + f 365 207 69 


om 


At the time these experiments were made, only one other r stock, r®, was tested 
on a small seale, with the same results. Neither of these stocks still exists. But 
it appeared advisable to repeat the experiments with dominant markers in the 
X chromosomes of the non-r grandparents, which should reappear in F; because 
they are derived from the grandparental (F;) PCN, though not visibly present in 
the rr F, 2. Unfortunately, it turned out that one line, r*, had a considerable 
amount of fertility of rr 2, while the other, r**, produced almost completely sterile 
rr 2 without any development of a PCN. But in a very large set of experiments 
r® gave a single case of PCN development which did not admit of any other ex- 
planation. (The almost complete sterility and the presence of the rr parents 
exclude any contaminations, apart from the types themselves.) The figures for 
this test are as follows: 127 sterile bottles with +20 rr 2; 198 fertile bottles 
with +20 rr 9: 3,111 rr 9, 1,641 r o& in all; sex ratio approximately 2 2:1 7; 
average flies in bottles: (a) from r°/r* Q , 21; (b) from r*/r*® 9 , 8; (e) from r**/r39 
2,0.2. The only case of PCN development was6 92 B+/+r,107B+. 

Although this is thus far the only successful experiment using a dominant marker 
in the X chromosome, the experiments with r* and r** yielded an additional number 
of instances of PCN development, with the result that rudimentary females pro- 
duced nonrudimentary offspring. The rarity of the phenomenon for r* and r*® 
crosses, as opposed to r?”", becomes visible in the total result of half-a-dozen in- 
stances while testing almost 8,000 rr 2. Actually, the following cases could not 
have been explained if the foregoing facts had not been known; thus they afford 
further proof for the analysis. Backcrosses (+ Oregon  B f r**) X r® were made, 
and the rudimentary females (= B f r**/r*; also heterozygous for B f) were mated 
in mass (ca. 20 9) tor® o&. Of 6 fertile bottles, 5 gave the expected result, i.e., 
rudimentary females and males, the females half with the B/+ phenotype, the 
males half B fr. There was an unusual extreme preponderance of females, which 
remains unexplained. One bottle contained completely aberrant offspring, namely 
(phenotypes): 95 9 rr, 113 9 B/+,3 9 ++4+,1 9 B/+rr,1 9 B/B, 88 2 
Bf, S88 %r, 14 A Br, 2% +. It is obvious, first, that rr 2 mated with r 7 
produced largely nonrudimentary offspring, second, that the big noncrossover 
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classes prove that the rr mothers had produced gametes half B f, half r; and, third, 
that unusually large numbers of offspring had come from a fertile rr 9. The non- 
rudimentary classes show, in view of the former discussion, that they must have 
originated from B f/r nuclei in the gonad of an rr 2, which could get there only by 
activation of the PCN nucleus of the grandmother. This nucleus should have 
been +/B fr; but, by crossing over in the F; 9 , it could have been in an individual 
egg Bf orr. Fertilized by r in the backcross, a rare R F, 2 could be B f/r—i.e., 
not rudimentary—and therefore could not have been the mother of the F;. If, 
however, in the F; 2 the crossover PCN had been activated, as described, an R F; 
rudimentary female would have nuclei in the ovary of the constitution B f/r, and 
this is borne out by the experiment. The crossover classes + and B/+ rr (pheno- 
types) would be the result of ordinary crossing over in the B f/r germ-track nuclei 
of the R F, 2, and also the male classes B f r and + with the surprisingly large 
Bfrclass. The one Bar female might be a purely phenotypical variant of B/+ 
or a double crossover. 

Exactly the same thing happened once more in another combination, namely, rr 
° from the backcross (+ (Canton) X Bf r)r® mated to Bf ro. Eleven times 
(against 22 sterile bottles) the expected results were obtained: females Bfr/Bfr 
and B f r/r, malesr and B f r and a few crossovers, all rather small numbers (all 
together, 64-80-26-25, 13 crossovers), as expected in the presence of r**. One cross 
was aberrant and clearly derived from two fertile mothers, one with small but nor- 
mal offspring, the other with larger offspring from a B f/r mother, the origin of 
which must have been the same as in the last case. The data are as follows: No. 
4645: females Bf r/Bfr39, Bf r/r 27, Bf r/Bf61, Bfr/+ 12, Bfr/B1; males 
r59, Bf r5, Bf 68. A mother B f +/+ +r mated to B fr & would produce 
the second and third classes of females and the first and third classes of males. 
A normal mother would produce the first group of females and part of the second 
group, and part of the first group of males and the second group. But this latter 
may also be a crossover class from the first females with PCN nuclei, as are the 
females of the fourth group. Again one double crossover female is found. 

By chance, an analogue of the last case also was found. In the last example, a 
crossover egg from B f r/+ — Bf +/+ + r was the beginning of the unexpected 
features. In the present case the backcross, from which rr 2 were selected for 
mating with r o’, was (B/B X r*) X r*, which, in the absence of r**, produces many 
fertile rr daughters. Out of 18 mass matings with +20 rr 92, 17 gave the expected 
result: 404 rr 9,127ro. But one bottle segregated into 65 rr 2, 31 ro, 83 + 
2,39+ o&. The first idea is contamination; but this is excluded, since the only 
error imaginable, contamination with a B/r sister, could not produce + offspring. 
The relatively high fertility and the sex ratio indicate that we have the same situ- 
ation as in the last case: the F, mother came from a crossover egg (B +/+ r—> 
Br/+ +) with the + chromosome. Again the PCN r/+ was activated, and the 
F, rudimentary mother had r/+ nuclei in the germ track. 

Finally, a few more complicated cases were found, proving that other variations 
exist in regard to the abnormal behavior of the PCN. The following case is re- 
markable, though the numbers are very small, due to the presence of r**. The 
mating is B/+ rr Q (phenotype) from the backcross (Canton X r*) & Bf r*9 
with r* o’, Expected are femalesrr and B/+ rr,malesr and B fr. Only 8 flies 
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of the expected types hatched; in addition, one female of the phenotype B/+ and 
one 2 +. Both of these unexpected females were successfully mated to r? o. 
The B/+ 9° produced daughters (phenotypes) half B/+ rr, half + andafewrr 9, 
sons half +, half B fr. This female, then, had been B f r/+. As the mother 
should have been B f r/r, a + chromosome can be obtained only if there was un- 
equal crossing over of r, which is improbable as far as we know, or if the + polar 
body of F; was fertilized by a B f r** sperm and then moved into the germ track 
as did the PCN in the other cases. The other unexpected female, the + 2, was 
mated to r° and produced females and males half +, half r. This female had then 
been r/+-, which is possible if it was derived from the F; PCN in the typical way. 
The interpretation seems unavoidable and shows that in the presence of r** three 
different kinds of nuclei went into the germ track of the F; 2: normal ones, a 
PCN nucleus, and a polar-body nucleus fertilized by sperm. In the same series a 
second identical case occurred; from B f r/r 2 mated to r® <, 4 expected flies orig- 
inated and 2 females + and 1 female BB or B/+ (intermediate in phenotype). 
The unexpected females were mated to r o& and again turned out to be of the con- 
stitution +/r and +/B f r, respectively. The occurrence of 2 such unexpected 
cytological abnormalities in 7 fertile rr X r matings (out of 17) is remarkable. 

Discussion.—The phenomenon here described consists in the activation of the 
normally perishing polar copulation nucleus (second and third polar bodies) under 
the influence of the mutant rudimentary, and its entering the germ track of the F; 2, 
thus making it possible that a homozygous rudimentary female produces r and + 
gametes. In the mutant r?”', females rr become fertile only if this process saves 
them from sterility. In the mutants r° and r*® the same phenomenon can occur 
but only rarely, while fertile rr 2 behave normally as a rule, though only a few 
offspring survive. Actually about 8,000 rr 2, namely, r°r°, r°r**, and r*r**, were 
tested in batches of about 20 females (all single-pair matings made resulted in a 
single son). More than half of these were sterile, and it can be estimated that not 
more than 400 individual females produced offspring at all. Among these, only 
5 cases of PCN development occurred, i.e., 1-2 per cent of reproducing females 
and only a few hundredths of 1 per cent in all females tested. As a whole, the 
phenomenon belongs to the group in which recessive mutants also affect the cy- 
tological behavior of the meiotic divisions in the egg, of which the best-analyzed 
case, both genetically and cytologically, is the mosaicism in the silkworm (Gold- 
schmidt and Katsuki*). But, while here the action requires the homozygous re- 
cessive mutant (for oily skin), in the present case the mutant, r, acts in the hetero- 
zygote. It is the specific features of the insect development, namely, germ-track 
determination by polar cytoplasm into which cleavage nuclei immigrate and fur- 
ther the formation of a PCN, which make possible the surprising details of the 
abnormal behavior. It is remarkable that different r alleles act differently, and 
it can hardly be doubted that other, not yet tested alleles may present other 
variants of the phenomenon. 

Summary.—The mostly sterile females homozygous for the rudimentary alleles 
rez and r®, if fertile, produce large numbers of wild-type daughters, a few such 
sons, and a few rudimentary daughters and sons. The + daughters are genetically 
+/r. The explanation is that only those eggs of rr 9 develop that are derived 
from r/+ nuclei in the rr ovary. These originate from the polar copulation nu- 
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cleus of the grandmother r/+ which exceptionally entered the germ track. If 
these r/+ eggs develop by pseudoparthenogenesis (autogamy)—a possibility known 
from Solenobia—the offspring is female and wild type. If occasionally a few such 
nuclei undergo normal meiosis and fertilization, the possible classes of daughters 
as well as sons appear in small and equal numbers. A repetition with the alleles 
r® and r*® gave only a few such cases in about 8,000 females tested; r** females are 
almost sterile, r? females relatively fertile. All together, about 1 female in 20 was 
fertile. In many of these, the grandmaternal X chromosome was marked with 
Bar, but only once were Bar offspring obtained from non-Bar rr 2 by the activation 
of the grandmaternal polar copulation nucleus. A few more such cases occurred 
with complications by previous crossing over, all together in less than 2 per cent 
of the fertile females. 

1C, J. Lynch, Genetics, 4, 501, 1919; T. H. Morgan, Z. indukt. Vererbungs., 7, 323, 1912; Am. 
Naturalist, 49, 240, 1915. 
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3 R. Goldschmidt and K. Katsuki, Biol. Zentr., 47, 45, 1927; 48, 39, 685, 1928; 51, 78, 1931. 
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Studies on man have contributed relatively little to the flowering of genetics 
which has characterized twentieth-century biological science. The reasons for 
this are well known: the length of the generation time in man, the impracticability 
of test matings, the difficulties in controlling environmental variables, the role of 
cultural overlay in molding and obscuring genetic predispositions. There are, in 
my opinion, numerous evidences today that despite these drawbacks, the study of 
human genetics has finally come of age, to the point where studies on man will 
yield results on a level of significance with those emerging from studies of other 
organisms well known to the geneticist. These results will emerge, not from looking 
to man for the qualities which have made Drosophila and Neurospora so useful to 
the geneticist, but, in part at least, from taking advantage of certain special 
attributes of man and the populations in which he gathers, attributes not shared by 
other organisms. 

In considering the factors which have contributed to the coming of age of human 
genetics, it is convenient to distinguish between developments which are related to 
the “natural” evolution of the subject and other developments which are currently 
contributing to a somewhat “forced’’ evolution. As regards the ‘‘natural’’ evolu- 
tion of human genetics, three contributory factors are outstanding. The first is 
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the recent rapid expansion of medical and biological knowledge, with the result 
that in the evaluation of a problem involving the genetics of man there can be 
brought to bear a set of techniques which far exceed those available for any other 
organism. The second is the development, during the last twenty years, of statis- 
tical techniques which enable a maximum of genetic information to be extracted 
from two-generation data, such as one must often be content with in human 
genetics. The third factor of importance in the recent ‘“natural’’ evolution of 
human genetics is the changing content of medical practice, with the increasing 
control of contagious and parasitic diseases and the growing relative importance of 
congenital, constitutional, and degenerative diseases, to all of which genetic factors 
are so important. 

As regards the “forced” evolution of human genetics, it is apparent that the very 
lively discussions of the past several years concerning the genetic problems created 
by increasing exposure to ionizing radiation have gone far to focus attention on the 
absolute necessity of acquiring as soon as possible a far better understanding of 
human genetics than we now possess. These same discussions have also stimulated 
many able geneticists to direct a greater portion of their creative energies toward 
human genetics than would otherwise have been the case, and to give serious 
consideration to the special problems of this field, a development from which the 
field cannot but benefit. The activities of the National Academy of Sciences 
Committee on Radiation Genetics are a notable example of this. 

My assignment today is to review the special problems inherent in the study of 
human genetics. An exhaustive treatment of this subject would obviously be 
extremely time-consuming. Under the circumstances, I have chosen to restrict 
my remarks to the special problems which arise in connection with the subject of 
today’s symposium, the evaluation of the genetic risks of increasing exposure to 
ionizing radiation. 

In the evaluation of this problem there are certain types of information which are 
indispensable, in the sense of supplying the figures fed into mathematical treatments 
of the subject. These are (1) the spontaneous rate of mutation in man, be it on a 
per locus or a per gamete basis; (2) the induced mutation rate per locus or per 
gamete per unit of ionizing radiation; (3) the ‘accumulation factor,” i.e., the ratio 
of deleterious mutations already present in the population to those arising spon- 
taneously each generation through mutation; and (4) the manner of action of 
natural selection on past and contemporary human populations. Let us now 
consider briefly some of the special problems inherent in furthering our information 
concerning man in each of these four areas. 

The Spontaneous Rate of Mutation of Human Genes. The spectrum of phenotypic 
effects resulting from gene mutation ranges from complete lethality to extremely 
minute effects, apparent only under special circumstances. Two of the standard 
techniques used in the study of mutation rates in sexually reproducing forms con- 
sist (1) in the use of a special strain of animals in which one or more chromosomes 


carry complex chromosomal rearrangements associated with so-called ‘marker” 
genes, of particular use in the detection of mutations assoviated with lethal effects, 
and (2) in the use of special strains homozygous for a number of recessive mutations 
involving known loci, of particular use in the detection of mutations at these specific 
loci. Neither of these techniques is available for man. One is therefore forced to 
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rely on the very simple device of counting the rate of appearance of a variety of 
mutant phenotypes in human populations and then computing the mutation rate, 
either on the basis of the gene replacement rate necessary to offset the reproductive 
disadvantage accruing to the phenotype or, in the case of dominantly inherited 
traits, on the basis of the proportion of sporadic cases. The former approach 
requires the assumption that the population is in genetic equilibrium or nearly so. 
In the case of phenotypes determined by recessive genes, the former approach also 
necessitates certain critical assumptions concerning the effect of the gene when 
heterozygous. 

In order to avoid making such assumptions, which may introduce indeterminate 
errors, the tendency in recent studies on human mutation rates has been to concen- 
trate on phenotypes determined by dominant genes and to base the estimate solely 
on the proportion of sporadic cases of the disease occurring in a given population 
over some arbitrary time interval. Now the situation is simpler, but still, before 
one can argue from the frequency of the trait in the population to the mutation 
rate, one must meet two important qualifications. 

First, one must establish how frequently the phenotype in question is not the 
result of a germinal mutation. This may be difficult, but to fail todo so may intro- 
duce an error. A single specific example will suffice. Retinoblastoma is a type of 
cancer of the eye, developing at an early age, which in many family studies has 
been shown to behave as if due to a single dominant gene with a high degree of 
penetrance. Some years ago our group carried out a state-wide survey in Michigan, 
with reference to the frequency of occurrence of “sporadic” cases—i.e., individuals 
with no family history... We then calculated a mutation rate on the assumption 
that each of these individuals represented a mutational event. We specified at the 
time that this estimate would be incorrect if some of these individuals owed their 
disease to nongenetic causes or to somatic mutations, but we had no evidence that 
this was the case. Subsequently, several studies have suggested that when the off- 
spring of affected individuals who have reproduced are studied, there are fewer 
affected children than would be expected if all affected individuals owed their 
disease to germinal mutation.” * Our own estimate was consequently too high. 
This potential source of error is not serious, in the sense that the means of dealing 
with it are both clear and feasible, although temporarily, as in our case, the necessary 
data may not be at hand. 

Second, for an estimate to be entirely valid, one must establish the fact that the 
mutant phenotype in question always results from mutation at the same locus. 
This is not possible in man at the present time, so that all studies to date may actu- 
ally be measuring the sum of mutation at several loci. However, developments are 
in sight which will in some measure meet this problem. These are of two types. 
Studies on the linkage relationships of the gene(s) responsible for a given phenotype 
will reveal how many different loci may give rise to mutations resulting in that 
phenotype. For instance, elliptocytosis is a condition in which the erythrocytes 
are oval rather than round. Many extensive pedigrees are known in which the 
trait behaves as if due to a dominant gene. It has now been shown that in some 
of the pedigrees this gene behaves as if closely linked with the locus or loci responsi- 
ble for the Rh serological reactions, while in other pedigrees there is no apparent 
linkage (see references in Morton‘). The logical conclusion is that there are at 
least two genetic loci at which mutation may give rise to the ovalocytic trait. The 
second approach to a determination of the number of different loci at which there may 
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occur genes resulting in a given phenotype involves careful clinical studies of the 
trait. For instance, it has now been shown that there are really two types of hemo- 
philia, both sex-linked in their inheritance, with the serum from an individual of 
one type able to correct the bleeding defect in the other type (see references in 
Graham, Macfarlane,® Schulman and Currimbhoy’). One type is about four times 
as frequent as the other. Studies were carried out some ten to twenty years ago 
on the rate of mutation to the hemophilia gene.* * If the mutations resulting in 
these two types of hemophilia are ultimately shown to involve different loci, then 
the mutation-rate estimates resulting from these early studies must now be revised 
downward, to take into account the fact that two different loci are involved (cf. 
Graham). As another example, diabetes insipidus, a condition characterized by the 
excretion of large quantities of very dilute urine, frequently behaves as if due to a 
sex-linked recessive gene. It has now been shown that in some pedigrees typical 
of sex-linked inheritance, the condition can be remedied by the administration of 
pitressin, an extract of the pituitary gland. In these pedigrees the primary defect 
would seem to be in the pituitary gland. In other pedigrees the condition does not 
respond to pitressin. Here the primary defect apparently resides in the renal 
tubule.’ Clinical tools again allow us to differentiate two distinct types of the 
disease, a fact to be reckoned with if mutation-rate studies on this condition were 
contemplated. 

Thus far, studies have been carried out on the rate of appearance through muta- 
tion of about a dozen different dominantly inherited phenotypes. The results are 
summarized in Table 1. It is an interesting fact that despite all the qualifications 


TABLE I 
EsTIMATES OF THE RATE OF APPEARANCE THROUGH MuTATION OF CERTAIN GENES IN MAN* 
Mutant 
GENES 
PER 
100,000 
CHARACTER POPULATION GAMETES REFERENCE 
Autosomal Dominants 
Aniridia Denmark 0.5 Mgllenbach, 1947 
Penrose, 1957 
fDenmark 4.3-4.9 Merch, 1941 
Chondrodystrophy Sweden 7 Boék, 1952 
Japan 9-14 Neel and Schull, unpub- 
lished 
Dystrophia myotonica Northern Ireland 0.8 Lynas, in press 
Kpiloia England 8-1.2 Penrose, 1936 _ 
Marfan’s syndrome Northern Ireland 0.5 Lynas, unpublished 
Microphthalmos without men- 
tal defect Sweden f Sjogren and Larsson, 1949 
Multiple polyposis of the colon Michigan “< Reed and Neel, 1955 __ 
Neurofibromatosis Michigan Crowe, Schull, and Neel, 
1955 
Pelger’s nuclear anomaly Germany Nachtsheim, 1953 
England . Philip and Sorsby, 1944 
Michigan ‘ Neel and Falls, 1951 
Germany Vogel, 1954 
Waardenburg’s syndrome Holland Waardenburg, 1951 
Sex-Linked Recessives 
He hilis England 2 Haldane, 1935, 1948 
— | Denmark 3.2 Andreassen, 1943 
(Utah 10 Stephens and Tyler, 1951 
) Northern Ireland 4.5-6.5 Stevenson, 1953 


Retinoblastoma ) 


Pseudohypertrophie muscular 
dystrophy 


(Germany 4 Becker and Lenz, 1955 
England 4.3 Walton, 1955 

* After J. V. Neel and W. J. Schull, Human Heredity (Chicago: University of Chicago Press, 1954); L. S. Pen- 
rose, in Proceedings of the First International Congress of Human Genetics (Acta Genet. et Statist. Med. [Suppl.], 6, 
169-182, 1956); A. C. Stevenson, World Health Organization Publ. (in press). 
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I have erected, the average of these estimates, about 3 X 10, is not very different 
from the average commonly quoted for Drosophila, about 2 K 10~°. However, 
because of these very qualifications, as well as certain other methodological difficulties 
involving Drosophila as well as man, it would perhaps be wise to regard this as an 
interesting coincidence rather than as a demonstration of similarity, until such time 


as much more extensive data are available. 

In certain considerations it is important to have an estimate of the “‘total’” muta- 
tion rate, i.e., the average number of mutations per gamete. In order to convert 
figures on the rate of mutation at specific loci to a ‘total’? mutation rate, one must 
make certain assumptions concerning the number of mutable loci present in man. 
This is an extremely thorny problem for Drosophila or Neurospora, and even worse 
forman. Nor does it seem that here we can anticipate rapid advances with respect 
to human material. For the present, human genetics, in problems involving gene 
number, will have to be guided by the findings in other animals. Fortunately, 
there is an approach to the question of ‘“‘total’’ mutation rates which does not involve 
gene number per se. Recently, Morton, Crow, and Muller! have described a 
mathematical formulation which enables one to estimate the total mutation rate per 
gamete per generation on the basis of a comparison of certain attributes of the off- 
spring born to consanguineous as contrasted to nonconsanguineous marriages. 
This is a technique which should be quite useful. Unfortunately, the data now 
available which could be fed into these equations are of limited validity. This is a 
situation which time will certainly rectify, since there is currently a great deal of 
interest all over the world in the outcome of consanguineous marriage. 

There is one final point I should like to make concerning studies of spontaneous 
mutation rates. It is a well-established fact that mutations with very minor pheno- 
typic effects far exceed in number the mutations with major effects, with which it is 
usually convenient to work in experimental material. Better techniques for evaluat- 
ing both the frequency and the effects of these ‘‘small mutations” are urgently 
needed. A recent development in human genetics may point the way to approaching 
these small mutations. It now appears that the distinguishing feature between 
sickle cell and normal hemoglobin may be a simple amino acid substitution in a 
polypeptide chain.'? This difference is determined by a single gene substitution. 
One may speculate that the rather marked effect on the properties of hemoglobin 
induced by this gene substitution in the case of sickle cell hemoglobin are due to the 
fact that this) particular amino acid substitution is especially critical. The tech- 
niques are now available for analyzing the hemoglobin of different individuals, 
literally amino acid by amino acid. Such techniques promise to give much needed 
insight into the amount of concealed genetic variability in protein structure and, 
by inference, into the frequency of undetected mutation. Such mutation may be 
of importance either through a cumulative effect or as supplying a significant source 
of variation when environmental conditions change. 

The Induced Mutation Rate per Locus per Unit of Ionizing Radiation. The 
difficulties in reaching conclusions concerning the sensitivity of human genes to 
ionizing radiation are readily apparent. There are a number of circumstances 
under which substantial groups of human beings have received relatively large 
amounts of radiation. However, in no instance is the amount of radiation suffi- 
ciently large that one can survey the population involved with respect to the rate of 
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appearance of specific new mutant phenotypes, with any reasonable expectation of 
detecting a difference between the progeny of irradiated and controls. Rather, 
the approach must be much less specific, involving changes in the vital statistics of 
the progeny, i.e., such end points as sex ratio, stillbirth and neonatal death rates, 
frequency of occurrence of malformation, age-corrected death rates, etc. The most 
ambitious attempt to date to come to grips with this problem has been the Academy- 
sponsored study on the children of the survivors of Hiroshima and Nagasaki.'* 
This study was inconclusive, in the sense that while no clear-cut genetic effect of the 
atomic bombings could be demonstrated, neither, because of the magnitude of the 
“observational error,”’ could effects of an extent comparable to those produced by 
the irradiation of experimental material be excluded. There is, however, a more 
positive and sometimes overlooked aspect of these studies. They permit us to set 
upper limits to the sensitivity of human genes to irradiation. Thus, one of the 
important developments in radiation genetics during the past decade has been the 
demonstration by Russell'* that the mutation yield per unit of radiation was some 
ten times greater for the mouse—the first mammal to be so studied—than for the 
fruit fly. It was conceivable that man was in turn conspicuously more sensitive to 
the mutagenic effects of radiation than the mouse. We believe that the findings in 
Japan tend to exclude this possibility. Other positive aspects of the Japanese 
studies are that they have been collected and recorded in such a fashion that the 
findings may for analytic purposes be combined with such other sets of pertinent 
data as may some day become available and that they yield normative data con- 
cerning such questions as the outcome of consanguineous marriage, which contribute 
to a better understanding of the population genetics of man. 

Recent developments in tissue culture’ hold real promise of a somewhat more 
direct approach to this question of the sensitivity of human genes to irradiation, 
albeit in somatic rather than in germinal cells. It is now possible to establish 
“clones” of tissue culture cells, all the cells in a given culture having been derived 
from a single ancestor cell. It thus becomes possible to irradiate suspensions of 
cells and then “plate them out,” as it were, much as can be done with bacterial 
suspensions. Mutant clones can be detected by morphological or antigenic criteria, 
and mutation rates thus established. 

The Accumulation Factor. The term “accumulation factor’ has been used to 
designate the ratio of the number of deleterious genes which, on the average, are 
already present in the individual’s gene complex to the number arising each genera- 
tion in consequence of mutation. This figure serves as an important yardstick, 
both in estimating the strength of selective factors in the past and in providing a 
background against which increases in mutation rate must be viewed. One way 
of getting at this problem is the study of both the rate of death in, and the physical 
and mental characteristics of, the children resulting from consanguineous mar- 
riage. B6édk" has recently presented the first unbiased, unselected, and carefully 
controlled data on this point; his limited material suggests that there is an average 
of a minimum of three detrimental autosomal recessive genes in heterozygous form 
per individua!. On theoretical grounds, the true figure has been estimated to 
be considerably higher."*-2° As noted earlier, there is much current interest in 
this topic, and important data will almost certainly be forthcoming in the near 


future. 
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A second approach to this problem consists in the methodical collection of data 
concerning the frequency of inherited diseases. An estimate obtained from such 
data will of course be a minimum estimate. At the present time, our information 
concerning the frequency of most inherited diseases is so poor as to make frequency 
estimates of this type of relatively little value. This is a situation which can 
certainly be rectified, although for many years in the future this approach will 
yield very minimum estimates. A third means of estimating the accumulation 
factor consists of the development of clinical techniques which will make it possible 
to detect gene effects which have in the past gone unnoted. There is increasing 
evidence that many so-called “recessive” genes have some manifestation when 
heterozygous. Individuals who possess undesirable genes which find no or very 
slight expression are sometimes designated as ‘‘genetic carriers.” We are now 
able to detect, with greater or lesser degrees of accuracy, several dozen different 
carrier states (review in Falls and Neel?'). Some of the most pertinent examples in 
this respect are drawn from the field of hematology, where the relative ease with which 
hematopoietic tissue and its products may be sampled has led to relatively rapid 
advances. While it is certainly not imminent that by subjecting an individual to 
a battery of tests we can come up with a precise statement concerning the number 
of “undesirable” genes he carries, it should at least be possible to evaluate the situ- 
ation as regards quite a number of specific possibilities in the near future. 

The Manner of Action of Natural Selection in Past and Contemporary Human 
Populations. Of all the gaps in the background knowledge necessary to a precise 
evaluation of the genetic risks of ionizing radiation, that gap which will some day 
be filled by information concerning the manner of action of natural selection looms 
in my mind the largest. Many current treatments of the problem assume that the 
combined efforts of the physician and the sanitary engineer, coupled with agricultural 
and transportation developments that render famine unlikely, have brought the 
processes of natural selection which have shaped the human race to a virtual 
standstill. It is difficult to discern a really sound basis for this rather sweeping 
opinion, since we know so little about the workings of natural selection in either 
primitive or contemporary populations. There can be no question of some im- 
portant recent shifts as regards selective pressures. But whether these shifts 
are of a magnitude or type which actually mean that the human species is no longer 
capable of evolving to meet rapidly changing conditions, as some would maintain, 
is tome a very moot point indeed. There appear to be many opportunities to study 
natural selection at work today.??. These range from studies on causes of individual 
fertility differences to the relationship between diet and health. Recent develop- 
ments with respect to the abnormal human hemoglobins have provided a striking 
example of how the frequency of an apparently undesirable trait can be maintained 
through the mechanisms of balanced polymorphism. It is of the utmost importance 
to determine the degree to which the human species may be “buffered” against 
adverse genetic developments by such genetic systems. 

Mathematical Models. UHaving discussed the principal types of information we 
need in order to evaluate the genetic risks of increasing exposure to ionizing radia- 
tion, I do not want to leave you with the impression that, if this information were 
available, we would have only to feed it into certain equations and the answers 
would emerge like sausages from the other end. Despite the brilliant contributions 
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of Wright, Fisher, Haldane, and, more recently, Kimura, to the mathematical 
theory of evolution, this subject is still in its early stages. Most statistical models 
in this field have dealt with events at a single locus represented in a population by 
two alleles, with one of the alleles bearing a fixed selective relationship to the other, 
or with a relatively simple departure from this situation. In truth, of course, 
genetic changes in populations in consequence of ionizing radiation are the result of 
events at many thousands of loci, each represented by multiple alleles, with the 
selective value of a particular gene a function of the genetic makeup of the organism 
at many other loci. As the fifth general development necessary to an evaluation of 
this problem, then, I would mention the development of a mathematical approach 
which goes far beyond present techniques in recognizing the complexity of the prob- 
lem and which is evolved at ail times to take advantage of present-day electronic 
computers. 

Concluding Remarks. In concluding this very brief survey of some of the prob- 
lems of human genetics, I would like to face squarely one practical issue. The most 
precious asset of mankind today is the germ plasm which has enabled him to develop 
the culture we all enjoy. If current developments—of which increasing exposure 
to ionizing radiation is only one—threaten the integrity of this germ plasm and are 
paving the way for the decline of our species, then here is a problem equal in magni- 
tude to any with which we are currently confronted. I think it important to em- 
phasize that the kind of data necessary to an intelligent evaluation of the problem 
will not issue from the investigations of a few scattered geneticists but will require 
a greatly expanded effort on both human and nonhuman material, involving many 
teams of highly trained personnel. The annual financial support of research in 
human genetics in this country today is probably well below the cost of a single 
large jet bomber. While I by no means believe that genetic disaster is just around 
the corner, there is no doubt that the problems with which we are confronted 
call for an order of magnitude of research effort far in excess of what in the past, 
in the field of biology, has seemed adequate. 

* Contribution to the Symposium on Genetics and Radiation Hazards of the National Academy 
of Sciences, under the chairmanship of G. W. Beadle. 
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THE SCOPE OF GENETICS* 
By Curt STERN 
DEPARTMENT OF ZOOLOGY, UNIVERSITY OF CALIFORNIA, BERKELEY, CALIFORNIA 


“The experiments which will here be discussed took their origin from artificial 
fertilisations which have been carried out on ornamental plants in order to obtain 
new color variants.”” This opening sentence of Gregor Mendel’s paper points to 
one of the humble origins of genetics: horticulture. Significantly, as late as 1901, 
it was the Royal Horticultural Society of London that sponsored the English trans- 
lation of the then newly rediscovered classic. But the practical desires of the 
lovers of flowers were not the only incentives to the birth of genetics. All through 
the eighteenth and nineteenth centuries the species problem occupied the minds 
of biologists. On the one hand, Linnaeus had defined species in a way which 
lacked operational usefulness: There are as many species as there were diverse 
forms created at the beginning by the Infinite Being. On the other hand, Linnaeus 
and other plant hybridizers attempted to probe into the nature of species by experi- 
mental means. Were hybridizations between different species possible? Could 
they be induced artificially? Did they ever occur in nature? What kind of organ- 
ism was a species hybrid? Was it sterile like the mule? If fertile, what kind of 
offspring did it produce? 

In 1859, Darwin’s Origin of Species proved to be the successful break through 
the confines of a static conception of species within which biology had accumulated 
so much of its unordered bulk of empirical facts. Then plant hybridization be- 
came a potential tool for progress in a new range. Mendel himself, who began his 
studies before the appearance of Darwin’s book but wrote his report six years after 
that event, by then had recognized clearly the role of genetics in this respect. 
The detailed analysis of successive generations of hybrids ‘‘seems to be the only 
right way,” he wrote, “by which we can finally reach the solution of a question the 
importance of which should not be underestimated in connection with the history 
of the evolution of organic forms.’”’! 

Genetics and evolution were indeed to become intimately and fruitfully asso- 
ciated with one another; but the scope of genetics was destined to embrace many 
other regions. From variations in color and form, and from the species problem and 
hybridization, the emphasis soon shifted to the permanent units which underlie 
the externally visible phenomena. Mendel’s essential discovery had been the 
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recognition that the hereditary contributions of two parents were not each unde- 
fined wholes which in the offspring mingled into a new whole. On the contrary, 
he showed that the hereditary makeup consists of many separable elements which 
recombine in many ways in the sequences of generations and always reappear in 
unchanged nature. Bateson called these units factors, “‘makers.’’ Johannsen, in 
1909, coined the term “gene.”’ For him, ‘‘the word gene was . . . completely free 
from any hypothesis."" It was meant “solely to express the fact that characters of 
the organism are conditioned by special, at least partially separable and thus in a 
sense independent ‘states,’ ‘factors,’ ‘units’ or ‘elements’ in the constitution of the 
germcells.”’ 

Johannsen’s reluctance to go beyond the statements that “we can only operate 
with the knowledge that such constitutional elements exist’’ and that “no concept 
concerning the nature of genes is at present sufficiently well founded” was already 
somewhat out of date when he formulated these statements. While only the past 
few years have led to specific insights into the nature of the genetic material, the 
stage was set in 1900 to give the Mendelian gene its material counterpart as part of 
visible elements of cells, namely, the chromosomes. Roux and Weismann? had 
prophetically speculated on the nature of the chromosomes as bearers of hereditary 
material. Sutton and Boveri showed how the specific findings of the Mendelians 
fitted in with the specific observations of the cytologists. Assuming that the 
chromosomes carry the genes, the demonstrable distribution of the chromosomes 
in successive generations accounts for the Mendelian laws. 

The chromosomes are limited in numbers. Few organisms contain as many as a 
hundred in their nuclei, and most organisms possess many fewer chromosomes. 
The numbers of genes are larger than those of the chromosomes. Thus the sep- 
arateness of the chromosomes is not sufficient to account for the separateness of 
genes. This was foreseen early by Boveri but was established only later, by Morgan 
and Sturtevant. Each chromosome is a linear assemblage of numerous genes whose 
separability by recombination depends on physical exchanges of parts between 
chromosomes. 

The Mendelian discovery that the hereditary material of an organism is not a 
homogeneous mass but is subdivisible into genes has often been compared to the 
discovery of atoms. There is, however, a striking difference in the way genes and 
atoms can be handled by the experimenter. The chemist or physicist can actually 
isolate his elementary particles, singly or en masse, and study their properties by 
themselves. The geneticist can do this only very incompletely. Each gene is part 
of a complex interacting system. A gene becomes meaningful as a biological unit 
in terms of a co-operating, not an isolated, material. The geneticist corresponds to 
a chemist whose knowledge of the properties of hydrogen and oxygen atoms 
would be based solely on the role they play within protein molecules! 

Great progress has been made during the last decade in defining the genic material 
in chemical terms. The essential role of nucleic acid has been established. The 
transformation phenomena in the Pneumococcus micro-organisms in which extracts 
of deoxyribosenucleic acid (DNA) gained from one strain can transmit specific 
genic properties to another strain have given us genes in test tubes. The chemistry 
of DNA is complex, but it is chemistry. The test-tube genes can be changed 
by radiation and chemical agents: they can be made to mutate. Some day we 
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will be able to define in physicochemical terms the difference between DNA con- 
figurations which may cause redness of a flower or whiteness, well-functioning of a 
brain or idiocy. Clearly the gene is no longer Johannsen’s materially undefined 
element. But to the biologist the DNA in the test tube is of transitory interest 
only. Did he not know of its role in the cell—before the cell had released it or 
after the cell had again received it—it would be of little concern to him. 

“The central problem of biology is the physical nature of living substance. It 
is this that gives drive and zest to the study of the gene.’”’ Thus wrote L, Stadler 
in 1954, shortly before his death. Genetics has come close to the solution of this 
problem. Fortunately, the scope of genetics transcends by far the problem as 
defined by Stadler. He himself would have persisted forever in the search for 
still deeper insights. The central problem of biology, if there is one single one, is 
indeed inseparably bound up with the physical nature of living substance; but this 
problem is on a level different from that of the physical nature of the gene. The 
central problem of biology has often been hinted at: the organization of physical 
matter to form an organism, the emergence of biological systems, of co-ordinated 
cellular activities, of development, behavior, evolution. 

Organization and emergence are words which admittedly do not lead to analysis. 
They would not have given us the knowledge of the DNA gene. But they are 
concepts in their own right, with their own usefulness. They call us back to the 
organism, which is the object of our study as biologists. We rejoice in the success 
of our analyses of parts, but the facts found are biologically meaningful only within 
the realization of their role in the whole. When the physical nature of the material 
basis of a living system has become completely known and the dynamics of the 
biochemical interactions within it can be followed in all their individual steps, 
chains, and cycles, the integrated nature of these processes, the possibility of their 
existence, and the history of their assembly will still be a cause of wonder. 

The scope of genetics as a central biological discipline rests on the fact that genes 
are those parts of biological systems which give continuity to them. The primitive 
experiment of cutting an amoeba in two and finding that only the nucleated part 
survives and multiplies has been refined until we know that it is the genic material 
on which survival and multiplication of the organism depends. Most of the cellular 
material may be removed without permanent harm. Regeneration soon takes care 
of replacement. A gene, however, once removed from a cell, is not regenerated. 
Within a cell, a gene always comes from a gene. Again, this attribute of the gene 
is part of the organization of the whole. The material structure of the gene 
seems to carry within it the potentiality of its replication, but actual genic multi- 
plication depends on the presence of a whole assembly of structures and processes. 

The discovery of the cyclical nature of many biochemical sequences in cellular 
metabolism has been one of the significant events of our times. Genes are agents 
in these cycles not only because they determine the existence of enzymes, without 
which the reactions would not proceed sufficiently, but because some of these 
reactions themselves must lead to the making of more genes. We admire a jongleur 
who tosses and keeps rotating several rings moving in space. Figuratively speak- 
ing, a cell consists of thousands of such rings functionally bound together not by an 
outside mover but by themselves. The rings are not all of equal rank. Destroy 
some, and the rest will re-create them. Destroy a single agent in others, and the 
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whole system will disintegrate. These single agents are the genes. Through 
them, cellular physiology is brought within the scope of genetics. Rightly, more 
than thirty years ago, Muller called the gene ‘‘the basis of life.’’ 

Development is of the essence of most living forms. The emergence of a multi- 
cellular organism from a single cell is one of the miracles which daily occur around 
us. The embryologist has studied this miracle by observation and experiment. 
Strangely enough, for a long time embryology and genetics remained separate in 
methods and concepts. The embryologist was interested primarily in those basic 
developmental processes which occur alike in all individuals of a species, while the 
geneticist was concerned with differences among individuals. Moreover, in an 
opposite way, another alternative between similarity and difference found the two 
disciplines without common ground. Embryology is confronted with differenti- 
ation, the origin of unlikeness of tissues and organs, while genetics has postulated 
a genic identity of the cell nuclei in organs as different as skin, brain, and intestines. 

The gap between the two disciplines is being narrowed at present. It has always 
been obvious that genes participate in developmental determination. How else 
would they exert their influence, by means of which we discover genes for growth 
and form? More specifically, developmental geneticists have studied the embryo- 
logical pathways along which different genes produce different end results. Start- 
ing from another area, the embryologists are in the process of resolving the macro- 
scopic and microscopic events of development into their constituent biochemical 
reactions. Embryology becomes an aspect of cellular physiology—and the role 
of genes in cellular physiology thus includes their role in processes which underlie 
growth and development. 

Even the paradox of the postulated uniform genetic constitution of cells not- 
withstanding their differentiation is being resolved. Initial regional differences in a 
developing system, as in an egg, have long been known. They are consequences of 
such processes as diffusion of nutritive materials into the cell and of metabolic 
gradients. These general differences serve as evocators of specific gene actions, so 
that certain genes come to control one type of differentiation in one part of the 
embryo and other genes another type of differentiation somewhere else. Genes are 
part of the specificity of developmental phenomena. Genes determine whether a 
brain will develop here and fail to do so there, whether a wing will form or a leg, 
an eye or an antenna. The scope of genetics reaches to the problems of growth, 
form, and their origin. 

Evolution early was close to the mind of the geneticist. Its study is now un- 
thinkable without the framework of genetics. The differences between species are 
largely if not exclusively due to differences in their genes. Evolution is based on 
mutability, the origin of new varieties of genes, on the distribution and recombina- 
tions of these varieties within the members of populations, and on the effects which 
these diverse genotypes have on the fitness of individuals and of populations as 
whole entities. The New Systematics and the Modern Synthesis in evolution, to 
borrow phrases from Julian Huxley, are witness to the part which genetics has 
taken in the transformation of these subjects. 

Evolution is a historical field. Just as the history of the earth, or that of astro- 
nomical systems, is based on the properties of elementary physical particles, so the 
history of life on our planet is based on the properties of its genic material. The 
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scope of physics and chemistry extends to astronomy and geology without depriving 
these sciences of their own spheres of endeavor. To say—as it has been said— 
that a lion or a palm is only some genes’ way of making more genes expresses deep 
insights which link the lion and the palm to genetics—although little work on the 
genetics of lions or palms has been carried out. But the lions and palms remain 
valid topics for biological study as lions and palms. How evolution brought forth 
these specific species is not explained simply by referring to their genes. 

The analysis of animal behavior has rightly been re-recognized as another central 
discipline within biology. The inheritance of behavior patterns and the inherit- 
ance of plasticity in the range of possible behavioral phenomena are being redis- 
covered in general and analyzed in detail. Genes determine the degree of explor- 
atory behavior in rats, intelligence or dullness in running through a maze, fear 
reactions in birds to an object resembling in outline a specific predator, or the de- 
tails of a complex song. This is amazing and intelligible at the same time. If 
chemicals can change intimately one’s reactions, if under the influence of drugs 
spiders change the pattern according to which they weave their webs, why should 
not genes as biochemical cellular agents do the same? If genes can help to organize 
the differential development of a complex nervous system, why should they not be 
able to incorporate in it the subtleties of inborn reactions? The wonder is that 
there 7s behavior. Genetics provides one of the analytical tools for its study. 

The scope of genetics, with its enlightenment of all branches of biology, makes it 
an obvious partner in the pursuit of applied biology. From horticulture to plant 
and animal breeding in general, the concepts and methods of genetics—be it indi- 
vidual or strain genetics in terms of hybridization and segregation, or population 
genetics in terms of equilibria and selection—form indispensable parts of procedure. 
Problems of pest control involve the existence of innately sensitive strains and the 
origin of resistant ones. Host-parasite relations depend on a balance of genetically 
determined interactions between the two components of the combating teams. 

Human genetics is both a basic and an applied science. Knowledge of the 
genetic makeup of the species Homo sapiens, in its individuals and varieties, is the 
latest gain in the long history in which anatomy and physiology, in ever more 
penetrating modes, have provided insight into the biological nature of Man and 
anthropology into that of men. Eugenics, the planning toward rational improve- 
ment of the genetic makeup of human populations, has been regarded as a utopian 
dream, and a dangerous one in addition. Many of our concerns have been recog- 
nized as being excessive or even unnecessary, and our eagerness to act precipitously 
has rightly decreased. But, while the attempt at improvement beyond the present 
average could well be relegated to the future, even the most determined opponents 
of positive eugenics have now realized with the fervor of converts that negative 
eugenics, the attempt to inhibit the increase of harmful genotypes, is an immediate 
necessity. 

In the affairs of mankind, as a community of explorers and of transformers, 
genetics has indeed come to occupy a significant place. 

* Presented as the introduction to a Symposium on Genetics held under the chairmanship of 
G. W. Beadle at the ninety-fourth annual meeting of the National Academy of Sciences in Wash- 
ington, D.C., April 15, 1957. 
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' Tt is amusing to note that the English translation, made under Bateson’s influence, changed 
Mendel’s careful wording “nicht zu unterschitzen’”’ to the overstatement “cannot be overesti- 
mated.” 

* The name of August Weismann has been repeatedly misspelled in recent years. It should 
suggest a wise man, not a white man. 


DYNAMIC PROGRAMMING AND STATISTICAL 
COMMUNICATION THEORY 


By RicHarp BELLMAN AND RoBEerT KALABA 
RAND CORPORATION, SANTA MONICA, CALIFORNIA 


Communicated by P. A. Smith, May 28, 1957 


1. Introduction—The purpose of this paper is to present some applications of 
the functional equation technique of dynamic programming" * to the study of some 
multistage stochastic decision processes arising in statistical communication theory. 

Our starting point is a paper by Kelly* in which it is shown that the rate of 
transmission, as obtained by Shannon from considerations of coding of information,‘ 
can be obtained from a certain multistage process with a suitable criterion function. 
In this paper we shall complete a result of Kelly’s and considerably extend the 
scope of the investigation. Further results and proofs of the theorems stated be- 
low will be presented in a subsequent publication. 

2. An M-Signal Noisy Channel.—Consider a noisy channel which is called upon 
to transmit any of M different signals, which we name 1, 2, ... , M, in succession. 
Let 


pi; = the conditional probability that the 7 signal has actually 
been sent whenever the 7 signal is received ; 

qi; = the probability that the 7 signal is received at any particular 
time. 


A gambler, upon receiving a particular signal, is required to place bets on what 
he believes the transmitted signal to have been. He is allowed to bet a quantity 
2; that the 7 signal was transmitted, subject to the restrictions that }>z, < 2, his 


L 


initial capital, and z; > 0. If he bets correctly, he then receives r,z;,, otherwise 
nothing. This process continues for N stages, with a payoff at the end of each 
stage. Assuming that the transmitted signals are independent of each other and 
that the gambler wishes to maximize the expected value of a function ¢(w) of the 
final total at the end of the process, the problem is to determine an optimal wagering 
policy. 

3. Dynamic Programming Formulation.—Let us define the sequence of functions 


v(x) = the expected value of ¢(w) obtained using an optimal 
N-stage wagering policy, starting with a capital of x. (3.1) 


for N= 1, 2,.2., anda > 0: 
Then the principle of optimality yields the recurrence relations 
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If we specialize the function ¢(w), we obtain a noteworthy result: 
THEOREM 1. /n the case where ¢(w) = log w, we have 


f(x) = log x + Nk, 
where 


M 


M M 

k=: > a4 Max > py log (v2 +1-)> 2). (3.4) 

i=l Be, S1j=1 l=e /; 

The optimal policy is independent of the number of stages remaining, independent 

of the quantity of money available, and independent of the sequence {q;\. It is deter- 
mined by the maximization in equation (3.4) 

A particular case of the above in which it is required that }>z, = x yields the ex- 


t 
pression — >> p, log p, of Shannon. This furnishes an interesting link with infor- 
L 


mation theory. The foregoing result resolves a problem left open by Kelly. 

There is an analogous result if ¢(w) = w*, a> 0. 

4. Generalizations —The results of the preceding section may be generalized 
in many directions, in particular, to time-dependent processes and to the case where 
there is a continuum of types of signals. In both cases the functional equation 
technique is applicable, and the analogue of Theorem 1 holds. 

Another interesting type of process to consider is that in which the p;; are fixed, 
but unknown, constants. For an expository account of the problems encountered 
in this area we refer to Robbins; ef. also Bellman*® and Robbins.” 

As we shall show in a subsequent paper, a number of these problems may be 
treated by means of the foregoing techniques. 

5. Correlated Signals —Let us now consider the case where the signals are not 
independent. Although a large variety of questions of this type may be formu- 
lated, the following discussion of a simple process will illustrate the general method 
that may be employed. 

Assume that there are only two types of signals, say 0 and 1, that the probability 
of correct transmission at any stage depends upon whether or not the preceding 
signal was transmitted correctly, and that the gambler bets at each stage a certain 
quantity of his resources that the signal he receives was actually sent. Let 


p. = the probability of correct transmission of the kth signal if the 

(k — 1)st signal was transmitted correctly. (5.1) 
rs the probability of correct transmission of the kth signal if the 

(k — 1)st signal was transmitted incorrectly. 


Define the sequence of functions 


f(x) = the expected value of the logarithm of the final capital ob- 
tained from the remaining k stages of the original N-stage 





Vou. 43, 1957 MATHEMATICS: M. HALL 


process, when one has a capital of z and the information 
that the N — kth signal was transmitted correctly, and 
uses an optimal policy. 
the corresponding expected value in the case where the 
N — kth signal was transmitted incorrectly. 
Then 
f(z) = Max [py—eiife-i(a + y) + (1 — pw—esidge-i(2 — y)], 
O<y<z 
gx(z) = Max [ry—eiife-i(a + y) + (CL — ry—esadge-i(a — y)]. 
O<sy<r 
[t follows inductively that 
f(x) = log x + ay, g(x) = logx + b, 
where a, and b, are independent of x. The recurrence relations for the sequence 
\a,, b,} are readily obtained from equation (5.3). 
1 R. Bellman, Dynamic Programming (Prineeton, N.J.: Princeton University Press, 1957). 
2 R. Bellman, ‘““The Theory of Dynamic Programming,’ Bull. Am. Math. Soc., 60, 503-515, 


1954. 
3 J. Kelly, ‘A New Interpretation of Information Rate,’’ Bell System Tech. J., 35, 917-926, 1956. 
‘C. Shannon, ‘“‘A Mathematical Theory of Communication,” Bell System Tech. J., 27, 379-423, 


623-656, 1948. 
5 H. Robbins, ‘“‘Some Aspects of the Sequential Design of Experiments,’”’ Bull. Am. Math. Soe., 


58, 527-536, 1952. 
®R. Bellman, “A Problem in the Sequential Design of Experiments,’’ Sankhya, 16, 221-229, 


1956. 
7H. Robbins, “‘A Sequential Decision Problem with a Finite Memory,”’ these ProcEEpDINGs, 


42, 920-923, 1956. 


SOLUTION OF THE BURNSIDE PROBLEM FOR EXPONENT 6* 
By MarsHauu HALL, JR. 
DEPARTMENT OF MATHEMATICS, OHIO STATE UNIVERSITY 


Communicated by Saunders Mac Lane, June 5, 1957 


The restricted Burnside problem for exponent 6 was solved by Philip Hall and 
They even found the order of the largest finite group of expo- 


nent 6 generated by k elements, this order being 


Graham Higman.' 


(2) + () 


’ 


203° ¢ 

+ * @ an ‘ 

where a = 1 + (k — 1) 3! and b = 1 + (k — 1)2*. A proof is 

sketched here that a finitely generated group of exponent 6 is necessarily finite 

solving the Burnside problem for exponent 6. The proof will be published in detail 
elsewhere. 

It has been shown by Levi and van der Waerden? that there is a group B(3, k) 

generated by k elements and of exponent 3 whose order is 3‘, K = k + (3) + ()), 
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and that every other group of exponent 3 with k generators is a homomorphic image 
of B(3, k). The following relations hold on commutators in a group of exponent 3. 


(x, y~') = (a, y)™ 
(t, y, 2) = (y, 2, 2%) = 
(z, y, 2, w) = 1. 


It is well known that a group of exponent 2 is Abelian, since, from a? = 
(ba)? = 1, we find ba = (ba)! = a~'b~ = ab. 

THEOREM. A group G of exponent 6 generated by a finite number k of elements is 
finite. 

Lemma 1. The cubes of the elements of G generate a normal subgroup M of index 
of most 3*, K = k + (3) + (§). 

This follows directly from the result of Levi and van der Waerden. 

Lemma 2. M is generated by a finite number of elements of order 2. The derived 
group M’ of M is of index a power of 2in M, and M’ is generated by a finite number of 
elements of the form abab, where a? = 1, b? = 1. 

M, being of finite index in a finitely generated group, is itself finitely generated, say 
by a, ..., Q@m. M is also generated by elements of order 2, namely, the cubes of 
elements of G. Thus each a can be expressed in terms of a finite number of ele- 
ments of order 2, and the finite number of elements of order 2 needed to express 
a, ..., @m Will be a set of generators for M. If M is generated by 1, ..., x, with 
x =1,i=1,...,¢t, then M’ is of index at most 2‘in M. Also, M’ is generated 
by the commutators 2;~!x;~'x,2; and their conjugates, and so by a finite set of these. 
Hence M’ is generated by a finite number of elements of the form abab, where a? = 
1, b? = 1. 

The rest of the proof is devoted to showing that M’ is of exponent 3 and hence 
finite from the result of Levi and van der Waerden, and even from the less precise 
result of Burnside’s original paper.* Since M’ is of finite indes in G, the finiteness 
of G will follow from the finiteness of 17’, providing the theorem. 

Lemma 3. If H is generated by x, ..., X, and tf every subgroup of H generated by 
Jour of the x’s is of exponent 3, then H is of exponent 3. 

Here we may take H as generated by x, ..., x,, with defining relations z* = 1 
for every z in a subgroup of H generated by four 2’s. Then B(3, n) is a homo- 
morphic image of H. Philip Hall’s collecting process applied to words in x, ..., Xp, 
together with the validity of relations (1) whenever 2, y, z, w are taken from the set 
of x’s and their inverses, will show that the order of H is at most that of B(3, n). 
Hence H is indeed B(3, n), and H is of exponent 3. 

Lemma 4. if H = }2x,, x2, 3, xs} is of exponent 6 and if each of the subgroups gen- 
erated by three of the x’s is of exponent 3, then H is of exponent 3. 

This is proved by calculations using relations such as ((2, x2).2324)® = 1, together 
with the relations on three of the 2’s. 

Lemma 5. If H = }x, a, b} is of exponent 6 and if x? = 1, a® = 1,68 = 1, vax = 
a~!, xbx = b~!, then \a, b} is of exponent 3. 

This is the most difficult lemma of all and involves some complicated calculations 
with relations. 

Lemma 6. If H = }x, a, b} is of exponent 6 and if x2 = 1,a° = 1,68 = 1, xax = 


a~', eba = b, then \a, b} is of exponent 3. 
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Since x(bab—!)x = (bab—')— and x(b—ab)x = (b—'ab)—', it follows from Lemma 
5 that {a, bab} and }a, b—'ab} are of exponent 3. With these facts and further rela- 
tions the lemma follows. 

Lemma 7. If H = }{x, a, b, c} is of exponent 6 and if x? = 1, a’ 
vax = a~', xbx = b-', xex = c—, then a, b, c} is of exponent 3. 

From Lemma 5 it follows that | a, b{, }a, c}, and |b, c} are of exponent 3. Thus 
x(a, b)x = (a~', b-') = (a, b). Hence from Lemma 6 } (a, b), ce} is of exponent 3. 
Also, x(aba)x = (aba)~', whence, from Lemma 5, } aba, c} is of exponent 3. From 
these and similar relations the lemma follows. 

LemMa 8. If x? = 1,a;° = 1,1 =1,...,n, xa = a,~', and if }x, a, ..., ay} is 
of exponent 6, then { ay, ..., d,| is of exponent 3. 

This is a direct consequence of Lemmas 3, 4, and 7. 

LemMa 9. Jf H = a, b, c} is of exponent 6 and a® = b* = c? = 1, then H’ is of 
exponent 3. 

We easily find that H’ is generated by the five elements 


a, = abab, a, = abcach, 
a acac, a, = abcbca, 
a; = bebe, 


and then show that the elements a;, 7 = 1, 2, 4, agas~!, a3~!a5, a3ajas~', a2 tajas, 


t = I, 2, 4, are transformed into their inverses by the element a and that they are of 

order 3. By Lemma 8 these elements generate a group K of exponent 3. But K 

is normal of index 3 in H’. Hence H’ is of exponent 9 and so of exponent 3. 
Lemma 10. Ifa? = b? = c? = d? = 1 andif H = {a,b,c, d} is of exponent 6 and 


a = abab, 8 = ceded, then T = } a, 8} is of exponent 3. 

Write 8, = ceded, 8, = acdeda. Then, by Lemma 9, }8;, 62} is of exponent 3. 
We show that each of a, 81:82~', 8:~'82, BiaB2—', B,~!aBe is transformed into its in- 
verse by the element a and is of order 3. By Lemma 8 they generate a group U of 
exponent 3. U is normal of index 3 in V = ! a, Br, Bot , whence V is of exponent 9 
and so is of exponent 3. Hence 7’ as a subgroup of V is of exponent 3. 

Lema ll. Ifa? = b? = c? = d? = e? = f? = 1 andif H = ja, b,c, d, e, f} is of 
exponent 6 and a = abab, 8 = ceded, y = efer, then T = | a, B, y} is of exponent 3. 

Write 8 = 8; = cded, B, = acdeda, y = y; = efef, y2 = aefefa. Then, by Lemma 
9, |B;, Bo} and } 1, yo} are of exponent 3. We show that a, W(61, y:)aW (Be, y2)~', 
and W(@;, y:)W (62, y2)~!, where W(u, v) is any word in two generators u and v, are 
transformed into their inverses by the element a and are of order 3. By Lemma 8, 
they generate a group U of exponent 3. U is normal of index 27 in V = } a, 81, 11, 
Bs, v2}. V/U is of exponent 3, and so V is of exponent 9 and so 3. Thus T as a 
subgroup of V is of exponent 3. 

LemMMA 12. M7’ is of exponent 3. 

This comes from Lemmas 2, 3,4, and 11. This proves J/’ finite, and so the finite- 
ness of G, which is our theorem. 

* This research was supported by the Office of Ordnance Research. 
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PARABOLIC EQUATIONS 
By JoHn NASH 
THE INSTITUTE FOR ADVANCED STUDY, PRINCETON, NEW JERSEY 
Communicated by Marston Morse, June 10, 1957 


This note outlines the proofs of theorems on the continuity of solutions of linear 
parabolic and elliptic partial differential equations. These a priori continuity 
theorems are important for their applications to nonlinear equations. Our original 
motivation was the existence problems for fluid flow, especially the general com- 
pressible viscous flow, and the turbulence problem. We were led to the problem 
solved here by a process of reduction, isolating what seemed the essential difficulty 
facing an attack on the general nonlinear equations. 

We study parabolic equations of the form 

of. or] oT Bes is "a 
2 or, [ ij(X1, Xe, ei dndids 3 | = aL F V-(¢ VT] = T, (1) 


i, j 


where C is a symmetric real matrix varying in space and time and constrained only 
by universal upper and lower bounds ¢: > ¢, > 0 on its eigenvalues. Such equations 
characterize diffusion or heat conduction. 

Turorem. [f T(x, t) is a solution of equation (1) satisfying |T| < B whent = to, 
then for tz > t, > to 


| , pee 3 a a dee a/(2a+2) 
| T (2, t,) - T (22, te) | < BA | ( ti) + ( os ) | (2) 
V ty 1 to ty a to 


where a and A are constants that depend only on c, ¢, and the space dimension n. A 
similar result obtains for an elliptic boundary-value problem. 

TueoreM. I[f T(x) satisfies V-(C-VT) = O in a region and T < B on the 
boundary, then 


: jn — a] yet 

| T(x) — T(®)| < Bp| Tr | (3) 
min (d(x), d(#)) 

where a is the a of equation (2) and D depends only on ¢;, c. and n. 

We attack via fundamental solutions, which are conveniently unified in a single 
characterizing function S(xtit). S is defined and positive whenever ¢t > 7. For 
fixed Z, 7 S satisfies equation (1) and has unit weight: S dx = 1, where dx con- 
ventionally stands for the volume element in the Euclidean space. For fixed x, t 
S is similarly a fundamental solution of the adjoint equation 

2 : os 
Ve [C (&, t) -V S| =lc OL ° 

For any solution 7’ of equation (1) bounded at time 7, we have a basic representa- 
tion identity for t > 7, 

T(x, t) = JS S(atét) T(z, 2) dz, (4) 


and also the special case 
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S(xtrt) = Jf S(xtzt) S(alet) di. 
The inequality 
| T(x, t) — T(x’, t’)| < max| T(z, 1)| f/| S(atel) — S(x’t'zl)| dz, (6) 
& 
which follows from equation (4), shows that continuity of the solution 7’ is equiva- 
lent to L, functional continuity of S in the appropriate sense. 

Our argument has three main stages, in which different inequalities are estab- 
lished. The major result of the first stage is important in itself and is used in all 
later stages. It controls the rate of “spreading out” of a fundamental solution and 
seems the key to the whole problem. We begin by considering a fundamental 
solution 7’ = S(xt00) with singularity at the origin. We derive, first, 


<a (7) 


As in inequality (7), constants that depend only on ¢, ¢2, and n will be indicated by 
the letter k. Relation (7) can be derived from equation (5) and an L, inequality, 


ST? dx < kt~"”. (8) 
This, in turn, follows from the inequality 
S'\9T |? dz > kL fT? dz)"t?” (f T dz)-*, 
| 


which holds for any positive function 7’, and another inequality derived thus: 
d : 
" ST? dz 2S TT dx =2 Sf TV:-(C-VT] dz 


—2 f VT-C:VT dx 
l 
= ms S T* dz > 2, S |VT\? dz. 

’ 
We define a moment M as 

M=frldr= f Vait+ ... + 10° T dz, (9) 
and the estimate to be established, 

M < kv, (10) 

is the first main result. We also define an “entropy” Q as 


Q=—S T log T dz. (11) 


From relation (7) it follows that 


Q>k+ (5) log t, (12) 


and from simple extremal considerations using only T > 0 and { T dx = 1, we 
obtain 


M > ke®"*. (13) 
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Differentiating equation (11), we find 
O= — f (1 + log T)T dx 
—f (1+ log 7) V-[(C-VT] de 
JS Vilog 7T)-C-9T dx 
JS Vilog 7)-C-¥ (log T)(T dx), 


so that 
oo) > Sf \C-¥ log T| 2(T dx) 


[S| C-¥ log T| (T dx)}? 
[S |C-VT| dx}? 


4° y nm | px 
Since, analogously, M = — f Vr-C-VT dex, and | Vr = 1, we have a basie fact: 


eo > (M)*. (14) 


The three inequalities (12), (13), and (14) together enable us to bound both M 
and Q both above and below as functions of time; in particular, we obtain relation 
(10). 

In the second stage of our work we consider two fundamental solutions 7 and 7, 
with adjacent singularities (7,, 0) and (22,0). We show that there is a positive fune- 
tion e(r, | — 2s ) such that 

S min (7), T:) dx > er, | ay — x|) >0 
or (15) 
S \T — T.| dx = S [(T + T. — 2min (7), T2)] dx < 2 — 2e(r, | x1 _ x2| ) 
when ¢t > r. Since the family of equations of the form (1) with C bounded as we 
assume is invariant under the transformation « = Ax’, t = \*t’, any a priori bounds 
applying to solutions must be preserved under this transformation. This gives us 
a principle of dimensional analysis. In particular, the form of relation (15) must be 
© ‘yy ial T , 
S min (7, T:) dx > «€ re) (16) 


ley i ve i‘ 


where the function ¢ depends only on ¢;, ¢2, and n. 
Our proof of this depends on a study of the quantity 


F(é, t) = S *{log(T; + 6) + log(T. + 8)] dx, (17) 


where ¢ is aC” function, 0 < ¢ < 1, and ¢ = 0 outside a compact region containing 
the singularities of the two fundamental solutions, and 6 is a small positive constant. 
For small 6, the following bounds on F are valid: 


(a) Fis and (b) F > k log é. (18) 


These depend only on equation (17) and the facts T > Oand f T dx = 1. If we 
assume, for purposes of contradiction, that JS min (T;, T2) dx < «fort < 7, and 
that ¢ is sufficiently small in relation to 6, then we can show that for ¢ < 7 
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FP > kt"(log 6)? + kt”? (log 6). (19) 


The derivation of relation (19) is somewhat involved; we use the differential equa- 
tion and relations (7) and (10), but it requires no special insight. Now relations 
(18) and (19) are inconsistent if 6 is small in relation to r. Therefore, there is a 
lower limit on ¢ as a function of 7 (referring to the hypothesis above on ( min 
(T,, Ts) dx), and since £ min (7, T2) dz is an increasing quantity, relation (15) 
and therefore relation (16) are established. 
The final phase of our attack obtains the inequality 
t 


a ae |t1 — 22| \* 
JS |\T, — T2| dz < 2 — 2e ose oe. Sew —I], (20) 
Ty — 2° Vt 
where A and a@ depend only on ¢, ¢, and n. This more specific result is obtained 
by iterative use of the result (16). For this we use relation (16) with 7 in a specific 
ratio to | 2, — 22|?, so that e(r/| a, — 2s 2) has a specific value @. Then we show 
in the iterative argument that there is a sequence of times ;, fo, . . ., t», . . . such that 
the two fundamental solutions satisfy 


S T; — T.| dx < 3 (1 a | for ¢ > tp. (21) 


2 
As a part of the inductive argument, we estimate a moment quantity 
N,= S rlT, — T2| de at time f,. (22) 


using relations (21) and (22) for time ¢t,_; and (4) and (10) to estimate N,, for time 
t,. N, enables us to control the time interval t,_,; — ¢, necessary for fulfilment of 
relation (21) for time t,+:. Finally, we see that 


tn < ky(ke)” (23) 


? 


from which, with relation (21), we obtain directly the special case of relation (2) 
where ¢t; = ¢: (space continuity). Then a quite simple argument using relations (5) 
and (10) enables us to reduce time continuity to space continuity and obtain the 
general form of relation (2). 

Elliptic boundary-value problems can be modeled as parabolic boundary-value 
problems, and we can use relations (10) and (2) to derive (3). Consider a parabolic 
boundary-value problem with a fixed boundary B bounding a region R and T pre- 
scribed on B for all times, thus determining 7 within # for all times. For such a 
problem there is a measure (Green’s) function with properties analogous to the 
fundamental solutions in free space. For each point and time (2, f2) in the region 
and previous time ¢, there is a positive measure with total weight 1 defined. It is 
clearer to divide this into a part p(22, t2) applying to the boundary cylinder B,, ;, 
between times ¢; and ft and a part o(22, fs, t:) applying to R at time ¢, (R,,). In 
terms of this measure the values of 7 are given as follows in terms of prior values: 


T (x2, to) - Sa, te T(z, t)d px, ts) (2, t) + Sr, T(x, ty) dO (xs, te, ty) (21). (24) 


The measure satisfies an idéntity involving the characterizing function: 


S(tefatih) = Sb, ,, S(axtayty) dpe, ») (2, 0) + Sr, S(xtayty) does, tr, 0)(%1)- (25) 
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In equation (25) the second integral can also be written as do /dx,, the conversion of 
the measure into a function, and x; can be any point in space. Using equation (25) 
with relation (10), we control the portion of the weight of the total measure in the 
boundary part dp in terms of the size of ~/t. — t, compared with the distance 
d(x) from x2 to B. Then we use equation (25) again to establish the approximation 
of da to S(xetort;) dx;. This reduces the problem of estimating the continuity of 7' 
at time t. to the free-space problem where relation (2) applies. The elliptic bound- 
ary-value problem corresponds to the parabolic one where everything is independent 
of time and a steady state obtains. The best continuity control is obtained by 
selecting t2 — t; in proper relation to | x2 — 22’|, d(x2) and d(x’), x2 and x’ being 
two points in R. Proceeding crudely, we obtain relation (3). 

Previously, relation (3) was known only up to two dimensions and was obtained 
by somewhat special methods, This limited the study of nonlinear elliptic prob- 
lems to one and two dimensions. Relation (2) was, at most, known only in one 
space dimension. The key inequality (10) contains more strength than we have 
used here to derive relations (2) and (3). 

In all our discussion we have ignored smoothness questions, because, in obtaining 
the results we obtain, the essential thing is that no quantitative smoothness assump- 
tions are used, while any qualitative smoothness, such as C*, assumed for the coef- 
ficients of equation (1) is quite harmless. 

The author wants to thank the Alfred P. Sloan Foundation, the Institute for 
Advanced Study (Princeton), Massachusetts Institute of Technology, New York 
University, and the individuals who helped in various ways. 


RADIOACTIVE FALLOUT 
By Witiarp F. Lipsy 
MEMBER, UNITED STATES ATOMIC ENERGY COMMISSION 
Communicated May 8, 1957 
1. INTRODUCTION 
The radioactivity produced by the fission reaction changes its characteristics 
continuously and rapidly following the explosion of an atomic weapon, and the 
conditions of firing are of extreme importance in determining the rate at which the 
radioactivity descends to earth. As a result, there are in general three different 
kinds of radioactive fallout, the relative importance of which is determined by the 
nature of the weapon, principally its yield, and the conditions of firing. The first 
of these three types is the local fallout, which is insignificant unless the fireball 
touches or comes close to the ground but which, in case the fireball does touch the 
ground, can amount to a major fraction, in some instances as much as 80 per cent 
of the total debris. This type of fallout consists of radioactivity which is carried 
down by the larger particles. It consists largely of matter drawn up into the fire- 
ball from the surface, which is either totally or partially vaporized. Under these 
conditions so much matter is vaporized by virtue of the fireball’s touching the 
ground that the particle sizes formed in the freshly cooled vapor are large. 
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The second and third types of radioactive fallout are world-wide in nature and 
consist of finer material ; they are divided according to whether the material happens 
to lie in the lower part of the atmosphere, the troposphere, where rain and weather 
phenomena occur, or in the higher part of the atmosphere, the stratosphere, which 
is free of such precipitating mechanisms. The tropospheric fallout occurs in a mat- 
ter of two or three weeks or a month or so. It occurs largely as a result of rain 
and snow, and water precipitation in general, and falls in the general latitude of the 
test site. The stratospheric fallout, in contrast, takes years. We are not com. 
pletely certain, but it appears that an average time of something like ten years, on 
perhaps somewhat less, is a reasonable figure, and during this time the distribution 
becomes nearly world-wide. When the stratospheric fallout manages finally to 
pass into the troposphere, it is quickly removed by the same type of mechanism 
that brings down the world-wide tropospheric fallout, namely, rain and moisture. 

The precipitating mechanisms consist in general of the collision of the tiny par- 
ticles with moisture droplets in clouds, together with the interception of particles 
by falling raindrops. The first mechanism was recently suggested by Dr. Green- 
field in connection with “Sunshine” problems. In addition to the scavenging 
action of rains and fogs, there is definite evidence for a considerable probability of 
pickup on direct contact of air with surfaces such as the leaves of grass and trees. 
Frequently, grasses are found to have higher strontium 90 content than would cor- 
respond to the soils in which they grow, and this is due undoubtedly to direct pick- 


up. 
The dissemination of strontium 90 and all fallout is greatly dependent upon the 
firing conditions. There is every evidence that important factors include not only 


contact of the fireball with the surface but the nature of the surface—whether it 
is land or water, the type of soil or the composition of the water, whether fresh or 
sea water. Also, the height to which the fireball rises is important—in particular, 
the height relative to the tropopause, the dividing layer between the troposphere 
and the stratosphere. Yield is the main consideration here. A rough rule is that 
megaton weapons push through the tropopause into the stratosphere, and kiloton 
weapons stay below the tropopause in the troposphere. 

Thus we see immediately that kiloton weapons deposit their fission products 
much more quickly than do megaton weapons. Of course, this is of less importance 
in so far as the long-lived fission products, such as strontium 90 and cesium 137, are 
concerned, but it is of more importance for the shorter-lived fission products. As 
a general rule, an air-fired kiloton weapon will deposit its radioactive fallout in a 
period of between two weeks and one month, on the average, after the detonation, 
whereas an air-fired megaton weapon will deposit its radioactive fallout over many 
years—on the average, about ten years. Thus the effects which are due to the short- 
lived fission products are larger for a given amount of fission energy release in 
kiloton weapons than they are for air-fired megaton weapons. Considering the 
average age of the kiloton fission products to be one month, the external gamma- 
ray exposure from a megaton of fission fired as, say, 50 bombs of 20 kilotons each 
would be 30 times that for a single bomb giving one megaton of fission energy— 
if both were fired well up in the air. The fission products from the small bombs 
fired in Nevada would fall in the latitudes 10° N. to 60° N. in about one month, 
while the larger bomb would give fallout over essentially the whole earth in about 
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ten years. For strontium 90 effects there is relatively little difference per unit of 
fission yield, since even the residence time in the stratosphere is small compared 
with the 28-year half-life of radioactive strontium and the 27-year half-life of 
radioactive cesium, which is produced at slightly higher yield than strontium 90 and 
which appears to be disseminated in about the same way. 

The contents of radiostrontium and radiocesium in the stratosphere are shown 
hy direct measurement to be roughly the same, although the radiocesium is some- 
what higher, possibly because of the slightly higher fission yield. The content of 
radiocesium in rain water is comparable to that of strontium 90. Also, the content 
of radiocesium in the human body as measured by Marinelli at Argonne and by 
Anderson and Langham at Los Alamos agrees well with the fact that it has an aver- 
age residence time in the human body of about five months, as compared with 
many years for strontium 90. The radiocesium data are very interesting because 
of their bearing on the fallout dissemination mechanism and the confidence with 
which we can establish the probable future behavior of radioactive strontium. 
The data confirm previous suggestions as to the dissemination mechanism, that 
is, we find that radiocesium fallout, except of the local variety, is carried down 
very largely in the form of moisture droplets and that there is some direct pickup 
by leaves and grass on surfaces. It is captured and held tightly by the top two 
inches of most soils, so that the water which falls and runs off in the form of rivers 
is clean by the time it has drained a short distance through soil. All this is very 
similar to the radiostrontium behavior. 

The plants pick the strontium 90, and radiocesium to a lesser extent, out of the 
soil and also off their leaves and take it into their systems. There appears to be 
a discrimination mechanism which operates in most plants, so that the strontium 90 
content of the plant is considerably less, relative to its calcium content, than in the 
case of the soil. On the average, the discrimination factor between the topsoil 
and plants against strontium relative to calcium seems to be about 1.4. When 
the cows eat grass, they further discriminate by about a factor of 7 in making milk; 
hence there is an over-all protection factor for strontium 90 from the topsoil to 
milk of about 1.4 K 7 or 10. Also, there is a further discrimination factor against 
strontium relative to calcium in the human body. This factor is not known accu- 
rately but is known definitely to be at least as large as 2 and is thought possibly to be 
as high as 8. Research is now in progress to settle this. Therefore, there is a 
series of protective factors which makes the concentration of radiostrontium de- 
rived from milk relative to calcium in human bone not over one-twentieth and 
possibly as little as one-eightieth of that in the topsoil. Of course, it should be 
pointed out that a considerable part of the fallout is picked up directly on the leaves, 
and to this the factor of 1.4 does not apply; thus, for this fraction of the fallout, the 
protective factor may be reduced to 14. Since milk is the source of most of our 
calcium, this means that the actual ratio of radiostrontium concentration in new 
human bones relative to that in the topsoil should approach these numbers. 

It must be realized that, although only a small part of the calcium is derived 
from vegetables and meat, a similar calculation must be made for this portion and 
the total average ratio obtained. It seems that the meat-vegetable over-all dis- 
crimination factor is about 10; if 20 per cent of the calcium is derived from such 
sources on the average, the average over-all factor will be between !/13 and 1/30. 
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The experimental data on new human bone in children appear to give a smaller 
figure, '/g, as mentioned later. 

A matter of importance in connection with the amount of strontium 90 which 
one would expect to be deposited in human bone as a result of atomic weapon de- 
tonations is the calcium concentration in the topsoil. Since calcium is so similar 
to strontium, it seems very likely, and the evidence confirms this, that high avail- 
able calcium content of the soil will reduce the probability of strontium 90 being 
taken up into the plants. Of course, this probably does not have nearly as great an 
effect on the uptake of the material which is picked up directly on the leaves. We 
might expect, therefore, that soils which are particularly low in calcium might 
show higher strontium 90 contents for the grasses grown on them. This is, in 
fact, so, and sheep and goats and cattle feeding on such pasture displays a higher 
strontium 90 bone content. 

How such calcium deficiencies in the soil should affect the strontium 90 uptake 
by the human population is a most important question. One sees immediately 
that food distribution systems are such that the food supply is derived from large 
areas, and that there is consequently a sharp reduction in the sensitivity of the 
human population to calcium deficiencies in local soils. This is brought out par- 
ticularly well by the data on the radium contents of human bones and their obvious 
lack of strong dependence on the radium contents of local waters. But for people 
who consistently drink milk from cows grazing on such ground there should be 
definite effect on the amount of radiostrontium uptake, and the effect should be 
proportional to the radiostrontium content of the milk. So the question resolves 
itself largely into this: What are the strontium 90 contents of the foods that people 
in such regions actually consume? We find, on inspection of the food-eating habits 
and calculation of the strontium 90 intake relative to calcium, that the increase in 
average strontium 90 concentration of the food due to the low calcium content of 
the particular soils can hardly be more than fivefold for a soil caleium deficiency of 
fifty fold. That is, whereas normal soil carried about 20 gm. of available calcium 
in the top 2.5 inches, a region with soil of only 0.4 gm. per square foot would pro- 
duce a human-body burden equilibrium of about five times that which the normal 
soil would produce. 

In order to understand the hazard of radiostrontium, which is generally agreed 
to be the most hazardous of the long-lived fission products, we try to establish the 
maximum permissible concentration both for occupational workers and for the 
population in general. These numbers have been set at 1 microcurie and 0.1 
microcurie for the standard man, respectively. That is, an occupational worker 
may carry | microcurie of strontium 90 in his body, whereas the average standard 
adult in the general public should not have over 0.1 microcurie of strontium 90. 
This last figure corresponds to a concentration of 100 micromicrocuries per gram 
of body calcium, or what we call 100 “Sunshine Units’’; that is, | micromicrocurie 
of strontium 90 per gram of body calcium is defined as | Sunshine Unit. 

Now we must try to see in some other way how our normal experiences can be 
brought to bear on the question, ‘‘How dangerous is atomic weapons testing from 
the point of view of radioactive fallout?” At the present time we have in our 
bodies about 0.1 or 0.2 Sunshine Unit, and children have about one-half of a Sun- 
shine Unit, A little later I will speak about the question of the variation from 
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these average values; but, assuming at the moment that these are the values, what 
is the threat or the hazard from these quantities? Obviously, they are much 
smaller than the 100-Sunshine-Unit tolerance figure mentioned above. To obtain 
a comparison with normal experience, let us consider the fact that we know in a 
genetal way the magnitude of the radiation levels to which we are normally sub- 
jected by the cosmic rays, the potassium in our own bodies, and the uranium, tho- 
rium, and potassium in the ground and in our surroundings. We know that these 
quantities amount to something like 150 milliroentgens per year for an average 
person in this latitude. But we also know that there are considerable variations 
with conditions. 

For example, a person living in a brick house may very well get 25-50 milli- 
roentgens per year more than one living in a wooden house, because of the natural 
radioactivity of the bricks. It is also very well known that, whereas at sea level 
in this latitude the cosmic-ray dosage is 37 milliroentgens per year, at 5,000 feet 
altitude, as in Denver, Colorado, the dosage from cosmic rays is 60 milliroentgens 
per year—a difference of 23 milliroentgens per year. What is this in terms of stron- 
tium 90 body burden? 

First, we must consider what part of the natural radiation, if any, is similar to 
the radiation of strontium 90 in biological effect, so that we can say without doubt 
and hesitancy that the physiological effects, whatever they are, will be the same 
for the same energy absorbed. Fortunately, the cosmic rays seem to fit this 
bill. In other words, we are at liberty to compare the cosmic-ray radiation dosages 
with the dosages from radiostrontium in our bone structure. The reason this is 
permissible is that the ionization density along the tracks of the mu-mesons which 
are the principal cosmic-ray components at sea level and at altitudes of 5,000 feet 
are nearly the same as those of the yttrium 90 beta rays, the principal radiation 
which radiostrontium emits; that is, radiostrontium has a radioactive daughter, 
yttrium 90, which emits a very energetic beta ray, and the ionization density along 
the track of this radiation is very similar to that of the mu-mesons of the cosmic 
rays and their disintegration electrons; it is generally accepted by health physicists 
and radiobiologists that radiations of the same ionization density have very similar, 
if not identical, biological effects for the same energy absorbed. The high energy 
of yttrium 90 gives it an average distance of penetration in tissue of 2 millimeters; 
hence any effect of local nonuniformity of deposition of strontium 90 in the bone is 
removed. The cosmic-ray exposure is, of course, uniform throughout the bone 
structure. Therefore, we can equate cosmic-ray dosage with strontium 90 dosage, 
and it is thus possible for us to say that the difference between one altitude and 
another is equal in effect, other effects being equal, to a certain number of Sunshine 
Units in bone. Now, to follow this thought through, 1 Sunshine Unit is equal to 
3 milliroentgens per year. Therefore, the difference in annual cosmic-ray radi- 
ation dosage between Washington, D.C., or any place at sea level in this latitude, 
and Denver, Colorado, is equal to 8 Sunshine Units, that is, 16 times the present 
body burden of equilibrium bone or bone near equilibrium as we see it in young 
children who are growing now. 

Therefore, we must see whether anything in our experience indicates that these 
differences are significant in terms of the occurrence of the principal effects ex- 
pected of radiostrontium, namely, leukemia and bone cancer. Now, of course, 
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when one looks for such vital statistics, one finds that they are very hard to ac- 
quire. However, the National Institutes of Health of the Department of Health, 
Education, and Welfare have given us statistics for the occurrence of leukemia 
and bone cancer for the year 1947 for the three cities New Orleans, San Francisco, 
and Denver. They are shown in Table 1. It is clear from this table that there is 


TABLE 1 
OccURRENCE OF BonE CANCER AND LEUKEMIA 
(New Cases per Year per 100,000 Population) 
Bone Cancer Leukemia 
Denver 2.4 6.4 
New Orleans 2.8 6.9 
San Francisco 2.9 10.3 


no obvious effect of altitude, and it is also clear that there are other factors which 
are noticeably more important than cosmic-ray dosage. Of course, there may still 
be a considerable effect of altitude hidden in large fluctuations caused by other 
factors which presumably are largely unknown, and we cannot say that this proves 
anything. It does, however, give us some assurance from normal experience that 
the effect of 8 Sunshine Units will not cause a detectable increase in bone cancer or 
leukemia. 

This fits well with the laboratory data on animals and the limited experience on 
humans with radium. That is, | microcurie, being 1,000 Sunshine Units, is still 
considered to be pretty safe on the basis of the laboratory data. It is set as a toler- 
ance for occupational workers, and it is therefore reasonable that 8 Sunshine Units 
should give an effect so small as to be very, very difficult to detect. It is, I think, 
helpful for us, however, to realize that the present body burden of strontium 90 in 
new bone from the weapons tests that have occurred in the past is equal to the in- 
crease in cosmic-ray intensity that goes with an increase of some 400 feet in alti- 
tude, a very small fraction of the difference in cosmic radiation intensity between 
Denver and sea level. Therefore, at the same time that we consider the possible 
effects of strontium 90 from such concentrations, we may deduce from our everyday 
ordinary experience limits on the effects to be expected. None of the evidence on 
the occurrence of bone cancer or leukemia as a function of altitude has given us 
any reason to believe that the present tolerance limits are in any way in error. 
The present body burdens in new bones are small compared to these limits. 

Separate from the strontium 90 effects are the effects of general gamma radiation, 
the radiation that is received mainly from outside the human body and which 
comes mainly from the very young fission products in the local fallout area, but 
which can come in smallest part from radiocesium accumulating on the ground in 
the case of the stratospheric fallout or, more importantly, from the shorter-lived 
fission products deposited by the tropospheric fallout. Of course, weapons tests 
are so conducted as to avoid exposures to local fallout; hence our present discussion 
of the effects of weapons will be restricted to the much smaller gamma-ray doses 
from the off-site tropospheric and stratospheric types of fallout. In time of war, 
of course, it would be the local fallout which would be of more direct concern, 
next to blast and thermal effects, and it is to this aspect of fallout that Federal Civil 
Defense Administration addresses itself in the main. In regard to nuclear tests, we 
have to study the effects on human genetics and the possible effects of such doses of 





764 PHYSICS: W. F. LIBBY Proc. N. A. 8. 


radiation on health. Let us again apply the criterion of normal human experience to 
this. Measurements have shown that the general average intensity of fallout gamma 
rays from tests is from 1 to 5 milliroentgens per year. Now the general magnitude 
of the effects to be expected from this can be compared with the natural radiation 
intensity. We find, as mentioned earlier, that such things as living in a brick 
house instead of a wooden house can produce as much as 25-50 milliroentgens extra 
dosage per year, that there are certain areas in the world where the average dose 
in this country of 150 milliroentgens per year is exceeded by tenfold, that people 
living on granitic rock as compared with those living on sedimentary rock receive 
about 70 milliroentgens per year more dosage due to the higher content of uranium 
and thorium in these rocks, and that people living at higher altitudes have a higher 
natural cosmic-ray dosage. Also, of course, we know that medical uses of X-rays 
can be considerably larger than any of these fallout dosages. 

We do have experience and valid evidence that the somatic effects other than 
cancer and leukemia, that is, the effects of radiation on ordinary human health, 
require dosages which are very much larger—of the order of 25-50 roentgen units— 
in order to be observed as changes in the blood, and 100-200 roentgens for injury 
symptoms, whereas the dosages we are speaking of, from test fallout, are about one 
hundred thousand fold smaller. 

As for genetic effects, these are extremely difficult to evaluate, since so little 
is known about human genetics. But, judging from experience with plants, in- 
sects, animals, and lower organisms, there is every reason to expect some genetic 
effects of radiation. The question is how much radiation is required for a given 
level of effect. There are a certain number of mutations in every new human gen- 
eration. Are these largely induced by natural radiation, or are they mainly of 
chemical, or rather biochemical, origin, or both? From a chemical point of view, 
it seems likely that not all the spontaneous mutations in the human or any other 
species are caused by radiation effects, because it seems likely that radiation acts in 
inducing mutations mainly via molecules which it generates in the human cell, and 
that the mutations are caused by these chemicals and therefore in a sense are 
chemical in nature. Now, if this is so, and the radiation-induced mutations are 
nearly always caused by chemicals which are produced in the first instance by radi- 
ation, then chemicals which are not themselves produced by radiation, but have 
other origins, can cause mutations; thus it seems likely that a major part of the 
natural or spontaneous mutations in any species is not radiation-induced. This 
point js an important one to settle, for the reason that we have to compare the 
effects of fallout radiation with the fraction of the natural spontaneous mutations 
due to the radiation to which we are normally subjected. In other words, if the 
normal mutations are all due to radiation, then the effects of the additional radi- 
ation from general test fallout, or from other sources of radiation such as atomic 
power or the medical uses of isotopes and X-rays, will be larger. It seems likely, 
and many genetic authorities agree on genetic grounds with this conclusion, that a 
major portion of the spontaneous mutations of the human species is due not to 
radiation but to other causes. Therefore, a fraction of the spontaneous mutations 
in the human species is taken as being due to irradiation. Now, what this fraction 
is, it is difficult to say; but Professor H. J. Muller has estimated that this might be 
10 per cent. Therefore, one estimates that the 150 milliroentgens per year from 
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natural radiation now causes about 10 per cent of the spontaneous mutations and, 
therefore, that the test fallout, if continued indefinitely, will, at the present level 
of about 1-5 milliroentgens per year, cause an increase in the natural spontaneous 
mutation rate of something like one-fiftieth of 10 per cent, or 0.2 per cent of the 
spontaneous mutations. In the extreme, if it should prove that all the spontaneous 
mutation rate is radiation-induced despite the chemical arguments, the effect would 
be ten times as great, or 2 per cent. Dr. Gordon Dunning, of the Division of Biol- 
ogy and Medicine of the Atomic Energy Commission, estimated 1.4 per cent in 1955, 
based on the assumption that the annual fallout was equal to the highest amount ex- 
perienced heretofore in any one year, and that 10 per cent of spontaneous mutations 
were due to natural radiation (Scientific Monthly, 81, 265, December, 1955). This 
effect is one which is comparable to moving to a slightly different locality and is 
much less serious than changing from one house to another or doing any of a dozen 
things. The only important point is that genetic effects appear only if large num- 
bers of people are subjected to them. Therefore, we would expect that the effects 
of large populations changing their environment—for example, living at a higher 
altitude or living in a region of naturally higher radioactivity—should cause genetic 
effects if test fallout does so. An examination of vital records should be made 
to test for such effects, and the Atomic Energy Commission is doing so as best it 
can. The United Nations Scientific Committee on the Effects of Atomic Radi- 
ation has been comparing the data on natural background dosages, and it is hoped 
that this study will be continued and that the search will be made for observable 
effects of variations in the natural background dosage, for it is certain that any 
effects due to gamma rays from fallout must be already present in much larger 
measure due to the natural dosage. 


II. VARIATION IN INDIVIDUAL STRONTIUM 90 BURDENS 


What is the likelihood that even if the average strontium 90 content is well 
within tolerance limits, a few individuals should exceed tolerance limits? Let us 
consider first the case which will ultimately hold, the situation of complete equilib- 
rium with the environment in so far as the strontium 90 burden is concerned. 
The only way we can make judgments about the expected individual variations from 
the mean concentration is by direct. experiment on human-body composition, not 
only for strontium 90 but for other analogous constituents. Most of the recent 
data on the strontium 90 body burden are from odd bits of bone removed during 
surgical operations, but fortunately we have actual data for the strontium 90 
content of the entire bodies of some several dozen stillborn children’ in the city 
of Chicago in the year 1953. A strenuous effort is now being made in Project 
Sunshine to continue this series and also to check the human-bone data by anal- 
yses of complete skeletons. We present the distribution of the strontium 90 
data for the stillborn children in Figure 1. Data for the occurrence of ordinary non- 
radioactive strontium in human bones also have been published.* These obvi- 
ously refer to the full steady-state condition and are obviously at least as nearly 
in equilibrium with the environment as the fallout radioactive strontium ever will 
be. These data are presented in Figure 2. The occurrence of radium in the 
human body also has been used, since it is chemically similar to both calcium and 
strontium and therefore is a bone-seeker, and because it is obviously also in steady- 
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state equilibrium. The data used were those of Palmer and Queen* (Fig. 3) 
And, finally, we use the recent data on occurrence of normal potassium in human 
bodies as determined by Anderson and Langham‘ at the Los Alamos Scientific 
Laboratory, presented in Figure 4. All these data show a normal frequency dis- 
tribution as indicated by the theoretical curves. The respective widths of the 
curves (standard deviations) are 36 per cent for radiostrontium, 40 per cent for 
normal strontium, 40 per cent for radium, and 18 per cent for natural radiopotas- 
sium. It is completely clear from these data that the distributions agree with 
one another in general shape and that the magnitude of the distribution of the 
strontium 90 contents of the Chicago stillborn babies was not in any way anoma- 
lous. Therefore, we shall take the distribution curve for radiostrontium to be the 
same as for the normal strontium data. The occurrence of nonradioactive normal 
ordinary strontium in the bones should certainly tell us what the equilibrium dis- 
tribution will be for radioactive strontium, and from it we should be able to learn 
the points about distribution which we cannot yet learn in any detail from the 
radioactive strontium itself. Turekian and Kulp noted in their study of normal 
strontium in human bone that in a given region the deviation from the average was 
about 34 per cent of the average, that is, for human bone from the regions Colorado, 
Texas, Cologne, Bonn, Venezuela, Chile, Vancouver, China, and India. In each 
instance the ratio of the standard deviation from the mean itself was taken and the 
average calculated, to obtain 34 per cent. Therefore, we take 34 per cent as the ex- 
pected standard deviation from the mean for a given locality for the eventual stron- 
tium 90 equilibrium burden in human bones. 

With this result we can, assuming a normal error-curve shape of the distribution 
of probabilities, answer the immediate question: What is the probability of an 
individual exceeding the tolerance even though the mean does not? On the basis 
of this analysis, we find that at steady state and in equilibrium the variation from 
the mear will constitute an error curve with a shape corresponding to the standard 
deviation, being one-third of the mean. Therefore, at steady state among people 
living in a given locality, only one person in about 700 will have more than twice 
the average strontium 90 burden, and the chances of anyone having as much as 
three times the normal burden will be about one in twenty million. 

Now, what about the nonequilibrium distribution, when the strontium 90 is 
finding its way into the biological system? Obviously, the burden will be much 
lower here, but the deviation from the mean will probably be much higher per- 
centagewise, particularly in adults, where most of the bone has been deposited 
before strontium 90 was produced. The present strontium 90 content- of adults 
depends very much on the growth rate and the metabolic activity of the various 
bones in the given individual’s body which happens to be sampled. However, the 
specific concentration of the strontium 90 deposited will not exceed that in new 
bone developed entirely in the present biological environment, i.e., the local concen- 
tration in adult bone will not exceed that for the whole bone in young children, 
whose total bodies are composed of the mixture of strontium 90 and calcium which 
now is present in food. Since the present child-adult ratio is about 4 to 1 for aver- 
age total strontium 90 content, the factor of concentration in adults’ active bone 
regions may be as much as fourfold greater than the whole-body average. Thus 
the apparent spread for random bone samples taken from adults should be very 
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large compared with the true equilibrium spread for these reasons. As equilibrium 
is approached, however, the spread must decrease very, very markedly. 

The data on human bones indicate a very wide scatter, but it seems extremely 
clear that the variation is a reflection of the fact that the main skeleton of adult 
individuals is not in equilibrium with the present food supply and that the vari- 
ations reflect the different rates at which the various bones in the bodies of indi- 
viduals are coming into equilibrium with the food supply in the general biological 
environment. A study of whole skeletons taken from one given locality, now 
under way as a part of Project Sunshine, will clarify the point about the variations 
among individuals in their rate of coming into equilibrium with the general biologi- 
cal environment. This study is being carried on in Dr. Kulp’s laboratory. 

It should appear from these studies that the variation from the mean of adults 
will be larger than the factor of one-third which apparently is normal for the types 
of equilibrium distribution considered above. It is, of course, very important to 
establish the truth of this prediction clearly. However, the general agreement 
in shape of the distribution curves for such widely different materials as normal 
potassium in whole bodies, radium, and normal elementary strontium in fragmen- 
tary human bone, and actual fallout radioactive strontium in the whole bodies of 
stillborn children give us good reason to believe that there is nothing extraordinary in 
the distribution of radiostrontium in human bone. 


III. VARIATION OF THE STRONTIUM 90 BODY BURDEN WITH LOCALITY 


Most important of the causes of variation of the strontium 90 content of indi- 
viduals with locality is, of course, the amount of fallout in a given region. The 
general rules about the intensity of fallout have been described above. For air- 
fired megaton weapons our present indication is that the fallout is almost world- 
wide, and for reasons of simplicity and in the absence of better information at the 
present time, we work on the model that this is a uniform distribution, over the en- 
tire world, of material that falls from the stratosphere. Further evidence and 
data on this point are rapidly being collected which will undoubtedly settle the 
stratospheric horizontal mixing question. 

At the present time, the general latitudes in the Northern Hemisphere between 
10° and 60° N. have the highest strontium 90 content. In the United States, 
which, because of proximity to the Nevada test site, has unusually high fallout, there 
are at the present time about 25 millicuries per square mile of strontium 90. For 
average soil this means a concentration in the topsoil of about 50 Sunshine Units. 
With the factors of discrimination mentioned above, this means that an equilibrium 
body burden between 1.7 and 3.9 Sunshine Units is to be expected. Actually, the 
present body burden in young children indicates that the lower value is probably 
more realistic. The present body burden in children—about 0.5 Sunshine Units- 
probably was derived from an average strontium 90 content in the topsoil of some- 
thing like 15 millicuries per square mile, or about 25-30 Sunshine Units during the 
time the strontium 90 was being acquired. Thus we find that the experimental 
value for the ratio between the body burden of young children and the average con- 
centration in the topsoil is about 50 to 1—rather closer to the higher range of the 
laboratory results than to the lowest range. Table 2 contains the latest data for the 
total strontium 90 fallout as measured in United States soils, and Figure 5 displays 
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these data graphically. The northern part of the United States has about 20-30 
millicuries of strontium 90 per square mile, the southern states somewhat less. 
The low figure of 7 millicuries per square mile for Grand Junction, Colorado, is 
probably due to local climatic and sample-site conditions. 


TABLE 2* 


H®ALTH AND Sarety LaBoratory 1956 Survey or U.S. Sorts ror Strontium 90 SAMPLES 
TAKEN BETWEEN OCTOBER 8 AND 13, 1956 
Deprun DIsSINTEGRATIONS/ Mc/Mi? 
SamPLinG SITE (INcHEs) Min/Gm Som Mc/Mi? Av. 
Albuquerque, N.M. 0-2 0.078 + 0.001 .5+0.1 3 
0.075 + 0.001 .2+0.1 
2-10'/ 0.008 + 0.002 440.9 3.4 11 
0.005 + 0.002 .4+0.8 
Atlanta, Ga. 0-2 0.35 +0.007 +0.3 
0.42 +0.009 : 
2-6 0.018 + 0.004 
0.021 + 0.003 
Binghamton, N.Y. 0-2 0.32 +0.007 
0.35 +0.007 
0.019 + 0.003 
.024 + 0.005 
Boise, Idaho - .23 +0.006 
.26 +0.006 
.012 + 0.002 
.015 + 0.002 
Des Moines, Iowa .31 +0.907 
.31 +0.007 
.028 + 0.002 
.024 + 0.003 
Detroit, Mich. .26 +0.006 
.27 +0.006 
038 + 0.003 
044 + 0.008 
Grand Junction, Colo. 10 +0.001 
091 + 0.001 
11 +0.019 
070 + 0.013 
2-10'/2 <0.002 
<0.002 
Jacksonville, Fla. 0-2 11 +0.009 
2-6 013 + 0.004 
020 + 0.005 
Los Angeles, Calif. 0-2 12 +0.008 
14 +0.009 
2-7 .009 + 0.002 
.006 + 0.002 
Memphis, Tenn. .27 +0.006 
26 +0.006 
.028 + 0.003 
029 + 0.003 
New Orleans, La. .24 +0.006 
22 +0.006 
009 + 0.002 
006 + 0.002 
New York, N.Y. 21 +0.006 
29 +0 .007 
072 + 0.004 
068 + 0.004 
Philadelphia, Pa. 17 +0.005 
16 +0.005 
029 + 0.003 
026 + 0.003 
Rapid City 8.D. 29 + 0.006 
34 +0.006 
053 + 0.004 
.045 + 0.003 
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TABLE 2*—(continued) 

Depru DISINTEGRATIONS 
Sampuinea Siri (INCHES) Min/Gm Soin Mc/M?? 
Rochester, N.Y. 0-2 0.22 +0.006 16 
2-6 0.013 + 0.002 2.520 
013 + 0.002 2.6+0 
Salt Lake City, Utah 0-2 32 +0.007 22 +0 
33. + 0.007 23 +0 
31 +0.007 22 +0 
2-8 016 + 0.002 §.7+0 
016 + 0.002 5.9+0. 
Seattle, Wash. 0-2 46 +0.011 17 +0 
44 +0.010 16 +0 
2-6 .051 + 0.007 9.421. 
.052 + 0.004 9.620. 
* Strontium extracted with 6 N HClat room temperature. Replicates represent individual soil aliquots 


taken after sampling and air drying. Each error term represents one standard deviation due to counting 
error. 
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The differential rates at which the fallout has been occurring probably are best 
measured by the so-called ‘‘pot collection” method. A bucket with vertical walls 
of appreciable height is placed out in the open and allowed to collect the total 
fallout for a given period, including the rain, snow, dust, etc. The bucket is left 
out whether it has rained or not and covers the total fallout for a given period. 
Figures 6 and 7 give the curves so obtained for the New York and Pittsburgh 


Numbers are in mc/mi? at individual site. 


Fig. 5.—Strontium 90 in United States soil (HASL, October 8, 1956) (HCI extraction method). 


areas, together with the estimated errors of measurement. It is interesting to 
note the changes in slope and to correlate them with the occurrence of test activities 
and the relatively short-lived tropospheric fallout. The minimum slopes which 
appear during quiet periods when no one is testing are the stratospheric fallout of 
which we have spoken, and these slopes, when we have enough pots operating all 
over the world, will, when taken together with the results of the measurements 
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of the amounts of radiostrontium and radiocesium in the stratosphere, give an 
accurate value for the stratospheric residence time and settle the mixing question. 

In addition to the intensity of fallout, the question of the fraction of the radio- 
strontium, and, for tropospheric fallout, the radioiodine of 8-day half-life, that is in 
assimilable form is an important one. So far, most fallout strontium appears to 
be completely water-soluble and therefore most assimilable, though continued 
tests on this point should be made. Direct leaf pickup of course promotes assimila- 
tion of the strontium, because the plant differentiation against strontium when it 
assimilates it from soil thus is avoided. Another factor is, of course, the concen- 
tration of available calcium in the soil. By ‘‘available’”’ calcium we mean calcium 
which is available to plants and not the total calcium in the soil. It is known that 
soils which are high in available calcium produce plants of lower radioactive stron- 
tium content; that is, the radioactive-strontium-to-calcium ratio in the plant is 
lower as a direct consequence of the lower concentration of radiostrontium in the 
available soil calcium. In addition, as mentioned previously, plants tend to prefer 
calcium to strontium with a discrimination factor of about 1.4. Sheep which grow 
in certain areas of Wales have shown concentrations in their bones approaching 
150 Sunshine Units, while sheep and cattle growing in the United States have 
hardly ever exceeded one-fifth of this. The Welsh soil in certain areas is very low in 
calcium and, as a result, grows grasses of high radiostrontium content. Of course, 
it is clear that fertilization with calcium will immediately relieve this difficulty; 
but in the absence of such fertilization, the question is: How serious is the effect 
of calcium deficiency in promoting strontium 90 pickup through the food chain? 

As was remarked earlier, there is an averaging which occurs in food distribution 
systems, and calcium-deficient soils are naturally rather poor producers; as a con- 
sequence, the weight of the food so produced is less than for a good, well-fertilized, 
well-balanced soil. This factor reduces the flow into the general food system of 
material of exceptionally high strontium 90 content. It will therefore probably 
be sufficient to consider the radiostrontium of milk, since milk is the main source 
of calcium, in order to test for the radiostrontium content of the food in given areas. 
Direct measurements have shown that a factor of 5 encompasses the total variation 
due to all factors, including calcium deficiencies in acid soils. 

The general intake must depend on the food distribution pattern, and the rela- 
tively small fluctuation in milk contents must reflect this. The number of indi- 
viduals who rely totally on the food output of soil of very low calcium content is 
very small indeed, but it must be true that these individuals, if they grew up on 
such a provincial, isolated farm, would have as much as 10-50 times the normal 
average strontium 90 content. The normal calcium concentration in soils in the 
United States is about 20 gm. per square foot for the top 2.5 inches, and about the 
poorest soil known has about 0.4 gm. available calcium per square foot for the top 
2.5 inches—a deficiency factor of 50. 

It is clear from a detailed examination made by the author for people living in 
-aleium-deficient areas with normal food distribution patterns that a factor of 5 
is about as large an effect as can be expected from a fifty-fold deficiency of calcium 
in the soil. The food from outside the calcium-deficient area reduces by a factor 
of about 10 the increase in strontium 90 pickup rate which would be expected 
from the calcium deficiency in thé soil if people lived entirely off the soil for their 
whole growing period of 20 years or so. 
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The food of lowest strontium 90 content is fish flesh, because of the great dilu- 
tion the fallout receives by the hundred meters of sea water above the thermocline, 
which rapidly mix with the fallout within a few hours or days. This means that 
the specific concentration of radioactive strontiuny or any other fallout constituent 
in sea water, is relatively very much lower than it would be in soil. For example, 
100 meters of sea water has 370 gm. of dissolved calcium per square foot, as com- 
pared with the average of 20 gm. per square foot for the top 2.5 inches of soil 
which absorbs and holds the fallout radiostrontium. Therefore, in principle, sea 
food and fish are lowest among foods in content of radiostrontium fallout. 

IV. EFFECTS OF CONTINUED TESTING AND GENERAL CONCLUSIONS 

In summary, then, we see that the present body burden of strontium 90 from 
atomic weapons tests in the United States corresponds to the radiation dosage to 
the bones which would result from a few hundred feet increase in altitude, and the 
present vital statistics show no observable effect on the occurrence of bone cancer 
or leukemia of much larger changes in altitude. The tolerance figure of 100 Sun- 
shine Units, or 0.1 microcurie for an average individual, or 100 micromicrocuries 
per gram of body calcium, that is recommended now is about two hundred times 
the present level for new bone in the United States, and it will not be exceeded by 
fallout from weapons tests in any foreseeable circumstances. 

The distribution of strontium 90 burdens among individuals for a given locality 
will be a normal error curve with a standard deviation of about one-third of the 
average concentration. This means that about one individual in 300 will have 
more than twice the normal average value for a given locality, and that about one 
in several million will have three times the average value. 

The effect of locality is more important, however, particularly in the effect of 
calcium deficiency in the soils. Careful consideration of this question indicates 
that there will be a very few individuals who show a strontium 90 content which 
is strictly inversely proportional to the available calcium concentration of the soil 
in their region. This is due to the fact that food distribution systems automatically 
average over a wide area, and people assimilate their calcium slowly. Most people 
drink milk and eat cheese and other calcium-bearing foods from a rather wide area, 
and this effect reduces by an estimated factor of 10 the potential effect of calcium 
deficiency in the local soils. 

On the basis of laboratory experiments, the human-body concentration of stron- 
trum 90 at equilibrium will be between 13 and 30 times less than that in the topsoil. 
The present data indicate that the higher figure is closer to the truth, and so we will 
be conservative in taking the figure of 20 for this ratio. Therefore, the present 
burden of 50 Sunshine Units in the topsoil of the United States may eventually lead 
to as much as 2.5 Sunshine Units in human bones but more likely will lead to about 
Ld, 

Of course, as testing continues and more fallout occurs, the levels will rise. The 
strontium 90 that still resides in the stratosphere at the present wiil fall out, ac- 
cording to our expectations, at a rate which just about compensates for the decay 
of the material already deposited, so that no great additional increase from this 
source is to be expected from weapons fired in the past. If the testing should con- 
tinue at about the same rate it has averaged over fhe last five years, then we should, 
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at equilibrium, after an infinite time, approach a level of 8 times the present rate, 
since the average life of strontium 90 is 40 years. This assumes that future testing 
will be conducted so as to give in each future five-year period the same as the last 
five have. And so we would expect in the United States at that time an average 
human strontium 90 concentration of 20 Sunshine Units, with the conservative 
factor of 20 between the topsoil concentration and the concentration in human 
bone, or 5 Sunshine Units if the factor of 80 is used. In other words, in the United 
States something between 5 and 20 Sunshine Units would be the equilibrium con- 
centration of human bones if testing continued indefinitely at the average rate of 
the past five years. This level would be approached after only a few decades. 
After 28 years the level would be half this equilibrium value, and after another 28 
years—56 years total—from an arbitrary beginning which we have set as 1952, 
we would expect in the year 2008 three-fourths of the equilibrium figures. So 
somewhere between 4 and 15 Sunshine Units of strontium 90 in human bones in the 
United States might result from the present type of testing being continued for 
the next 50 years. 

In those particular areas in the world where the soil is low in calcium, this level 
might go fivefold higher. At the present rate of testing, we might indeed approach 
the figure of 100 Sunshine Units, the tolerance limit for large populations, at the 
beginning of the twenty-first century for these certain limited regions in the world. 
The observed conditions in these regions could be relieved, however, by fertiliza- 
tion of the soil with calcium, using either calcium nitrate or lime, as appropriate 
from other considerations. 

Project Sunshine continues to study, the problems of world-wide fallout—the 
stratospheric inventory of radiostrontium and fadiocesium, the occurrence of these 
isotopes in the soils and water and the biosphere all over the earth, the biological 
effects at certain levels of contamination with strontium 90, and to a lesser degree 
with cesium 137, and the possible genetic effects of the low gamma-ray dosages asso- 
ciated with world-wide fallout from atomic tests. All these are studied not only 
with the object in mind of devising methods of protection against atomic warfare, 
but also with the thought of possible application in the remote event of industrial 
accidents which may happen in connection with certain of the peaceful uses of 
atomic energy, particularly atomic power. Certainly an understanding of the basic 
principles of world-wide fallout is applicable to the control and safe handling of 
isotopes. All of this is done in collaboration with the United Nations Scientifie 
Committee on the Effects of Atomic Radiation, and it is to be hoped that, as the 
data appear, all the countries in the world will join together in this international 
effort to understand better the effects of the great new fact in life, the nuclear atom. 

1W. F. Libby, “Radioactive Strontium Fallout,’”’ these ProcEEpINGs, 42, No. 6, 365-390, 
1956; University of Chicago, Project Sunshine Bulletin No. 12, August 1, 1956. 

? K. K. Turekian and J. L. Kulp, Science, 124, 405, 1956. 


’ Hanford Report, HW-31242. 
4. C. Anderson, R. L. Schuch, W. R. Fisher, and W. Langham, “Potassium and Cesium 


Radioactivity in People and Foodstuffs” (in press). 





THE LOCALIZATION OF BASIC PROTEINS IN THE NUCLEI OF 
LARVAL DROSOPHILA SALIVARY GLANDS* 


By Epwarp ©. Horn anp Catvin L. Warp 
DEPARTMENT OF ZOOLOGY, DUKE UNIVERSITY, DURHAM, NORTH CAROLINA 
Communicated by J. T'. Patterson, June 14, 1957 


The identification and specific localization of chemical constituents within the 
cell nucleus have lagged considerably behind the qualitative and quantitative 
characterization of these classes of substances, with the possible exception of one 
component, the deoxyribonucleic acid (DNA). In addition to the latter, the 
nucleus contains ribonucleic acid ? (RNA), histone or basic protamine*® (BP), and 
at least one other type of protein, the chromosomin of Stedman and Stedman,‘ the 
tryptophane-containing protein (Tr.Pr.) of Mirsky and Pollister,’ or the non- 
histone protein fraction which Hamer*® has subjected to comparative amino acid 
analysis. The localization of these substances in the live nucleus is dependent 
almost wholly on inference from a variety of analyses and observations made on 
fixed preparations. The much disputed and maligned Feulgen reaction points 
to the DNA as a major constituent of the chromosomes. Studies with ultraviolet 
absorption spectroscopy,' enzymatic attack, and pyronin staining? indicate that 
the RNA of the nucleus is predominantly nucleolar. The distribution of the BP 
seems to be much less precisely delineated. For example, Caspersson’ has reported 
its presence in the chromosomes and nucleolus of Drosophila salivary gland nuclei 
from his ultraviolet absorption studies. Vincent,’ from his chemical studies of the 
isolated nucleoli of starfish oécytes and of maize, could find no basic protein. 
Bloch and Godman® have reported a proportional quantitative relationship be- 
tween DNA and histone synthesis per rat fibroblast or liver cell. They have also 
shown cytochemically that, in a general manner, DNA and histone have the same 
nuclear distribution. 

In the present study the Feulgen reaction for DNA localization, combined with 
the fast-green-staining method of Alfert and Geschwind" for basic protein, has been 
employed with Drosophila salivary gland material to determine more precisely the 
localizations of BP with respect to DNA. 

The method employed is essentially that of Bloch and Godman: salivary glands 
of D. virilis fixed with trichloroacetic acid (TCA) were partially crushed under a 
cover slip; these preparations provided some whole nuclei. After the slide was 
frozen on a block of “dry ice,” the cover slip was flicked off and the preparation 
stained directly by a modified Feulgen technique following 24 hours’ lipid extrac- 
tion in 95 per cent alcohol. Other salivary squashes were prepared similarly, using 
the conventional acetic acid fixative that permits rupture of the nuclei, yielding 
fully extended chromosomes. After a Feulgen treatment in which TCA was sub- 
stituted for HCl, photographs were taken of a number of chromosomes or of whole 
nuclei; the position of each was recorded on a calibrated mechanical stage. The 
preparations were then re-exposed, and, after TCA extraction (15 minutes at 90° 
C.) to remove the DNA, were stained by the fast-green technique of Alfert and 
Geschwind. The recorded fields were then rephotographed. 
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The results are illustrated in the accompanying figures. In both the ruptured 
(Figs. 1-3) and unruptured nuclei (Fig. 4) there was a clear band-for-band stain- 
ing correspondence between the Feulgen-treated chromosomes (Figs. 14, 2A, and 
3A) and their DNA-extracted fast-green-stained remnants (Figs. 1B, 2B, and 
3B). Diffuse heterochromatin of the chromocenter stained after both treatments. 
The nucleolus was completely unstained by the Feulgen reaction, although its 
presence was frequently indicated by the Feulgen-positive reaction of the nucleolus- 
organizing region (Fig. 1A, arrow). After the fast-green treatment, the nucleolus, 
which often appeared vacuolated, and its organizing region were clearly evident 
(Fig. 1B). Their position invariably was neither against the nuclear membrane 
nor in a central position but rather near the chromocenter. 

The whole nuclei (Fig. 4) clearly showed a heavier fast-green reaction in the 
region of the nuclear membrane; the exact location of this concentration could not 
be determined, primarily because the gland cytoplasm also presented a strongly 
positive fast-green reaction. In unbroken nuclei the materials of the interchromo- 
somal space gave no reaction with either procedure. The placement of the ends 
of chromosomes in whole nuclei showed no detectable pattern. 

A number of these observations merit further comment. The precise corre- 
spondence of histone and DNA in the banded portions of the salivary gland chro- 
mosomes confirms Caspersson’s’ conclusion drawn from ultraviolet absorption 
studies on Drosophila salivary glands. Furthermore, Caspersson reported finding 
histone in the nucleolus, a discovery also confirmed by our data. The failure of 
the chromosomal interband regions to stain by either cytochemical reaction does 
not, of course, prove the absence of nucleic acid or basic protein in these positions 
but certainly indicates a lower concentration. Vincent’s'' argument relating what 
appears to be an unusually high protein content (40-90 per cent) reported for the 
nucleolus with its staining properties does not seem to apply to these studies. He 
argues that with such concentrations the abundance of diamino acids per unit 
volume will be greatly increased and that therefore any reaction based on the de- 
tection of the amino groups (e.g., Sakaguchi) will always be a relatively strong one. 
Alfert and Geschwind” indicated that the intensity of their fast-green-staining 
reaction seemed less dependent on the absolute number and type of basic group 
involved than on the over-all balance of acid to basic groups. The positive reac- 
tion of the nucleolus, then, in these studies is interpreted to indicate the presence 
of BP; the concentration would be a foolhardy guess. 

Vacuolated nucleoli have been reported (see Vincent'! for review) under other 
circumstances. It should be noted here that the fast-green-stained material shows 
vacuolation of the nucleolus particularly clearly in preparations of the whole 
nuclei. The vacuoles vary in size and in distribution throughout the nucleolus. 

The similarity of staining reaction demonstrated by the diffuse heterochromatin 
of the chromocenter to that of the euchromatin lends support to the notion that 
these two types of chromatin are fundamentally alike chemically, differing perhaps 
only physically and temporally. 

Since the interband regions must be accounted for chemically and since the 
remaining ill-defined protein fraction (chromosomin, Tr.Pr., acid-insoluble pro- 
tein, residual protein, non-histone protein) remains unaccounted for in the nucleus, 
at least these portions of the chromosomes might by inference be thought to consist 





Fics. 1-4.—Photomicrographs of Drosophila salivary gland preparation stained first with the 


Fuelgen reaction (Figs. 1A, 2A, and 3A) and then, following DN extraction, with fast green at 
pH 8 (Figs. 1B, 2B, and 3B). ’ Figure 4 is of a whole nucleus prepared with fast green at pH 8 
after DNA extraction. The arrow in Fig. 1A indicates the nucleolus-organizing region. 
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of this fraction. Such a scheme is not meant to imply limitation of this protein 
fraction to the interband region but rather should lend credence to a chromosome 
structure consisting of a “backbone” constructed of this protein along which the 
nucleoprotein complex is distributed finitely with spatial and chemical irregularity. 

The large interchromosomal spaces seen in fixed intact nuclei would be predicted 
from the work of Buck and Boche™ on living salivary glands. They observed 
simultaneous chromosome shrinkage and nuclear swelling in live material as it 
was fixed with acetic acid. The present material cannot under the circumstances 
provide evidence for or against the existence of ‘‘nuclear sap.’”’ One should be 
reminded, however, that Buck and Boche characterized the living salivary chromo- 
some as “‘behaving chemically somewhat like a sponge, in that it is capable of taking 
in and giving out an astonishingly large proportion of its fluid content, without 
permanent alteration of its fundamental internal structure.” If at fixation, with 
its concurrent shrinkage of the chromosomes, materials are “given out,” then they 
must not be Feulgen-positive nor fast-green-positive; in other words, at fixation 
neither DNA nor histone is released to the interchromosomal space in quantities 
detectable by the cytochemical techniques employed. 

The fast-green-positive reaction of the cytoplasm of these glands is especially 
puzzling. The Alfert and Geschwind technique has been applied by us to a wide 
variety of tissues, both embryonic and adult, but never with the strongly positive 
cytoplasmic reaction seen for salivary glands. It is possible, of course, that some 
wholly unrelated substance elaborated by these larval structures produces a posi- 
tive reaction following the application of the fast-green technique (the results are 
identical whether or not the Feulgen reaction precedes the fast green). A modi- 


fied Sakaguchi reaction" for the estimation of arginine showed the cytoplasm to 
be strongly positive, suggesting by inference that the positive reaction of the 
cytoplasm to the basic fast green was due to the presence of histone outside the 


nucleus. 

Summary.—The successive treatment of Drosophila salivary gland nuclei with 
the Feulgen reaction for DNA and the fast-green technique of Alfert and Geschwind 
for basic protein has demonstrated an exact correspondence in position of these 
two chromosomal constituents. The nucleolus also appears to contain some basic 
protein. 

We are especially indebted to Miss Mary Elizabeth House for preparing the 
majority of the tissues and for running the eytochemical reactions. 


* This work was supported in part by research grant RG-4203 from the National Institutes of 
Health, Public Health Service. 
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A LABORATORY STUDY OF THE EFFECT OF REMOVAL OF NEWBORN 
ANIMALS FROM A POPULATION* 


By L. BAsth SLOBODKIN 
UNIVERSITY OF MICHIGAN, ANN ARBOR, MICHIGAN 


Communicated by G. E. Hutchinson, June 18, 1957 


If Daphnia are maintained in the laboratory under conditions of constant tem- 
perature and constant food ration, a characteristic pattern of population fluctu- 
ations is observed.!~* Population size at any temperature is determined exclusively 
by food increment.* Natural populations of Cladocera also seem limited by food 
supply. 

Removal of a constant percentage of the newborn animals from a Daphnia pulerx 
population reduces the size of the residual population as a nonlinear function of the 
removal rate. The number of animals removed (i.e., the ‘‘yield,” if we think of the 
Daphnia population as a fishery) is inversely proportional to population size. 
Total fluctuation, measured as the mean absolute difference between successive 
censuses of a single population, with the census made prior to the removal of the 
newborn, is independent of removal rate. In other words, if the number of animals 
in the population is counted every four days, the mean absolute difference between 
this number and the corresponding number determined four days earlier or later 
from the same population is independent of the number of newborn animals which 
are removed from this population on the completion of each census. 

Fluctuations in residual population size are reduced in proportion to the removal 
rate of the newborn. That is, the mean absolute difference between the number of 
animals remaining in any particular population after removal of the newborn at 
any census and the number remaining in the same population after removal of the 
newborn at the subsequent or previous census is smaller, the higher the percentage 
removal of newborn. 

Let Vy be the total variability in numbers, computed as the mean absolute 
difference between successive censuses of a single population, prior to the removal 
of newborn; V p, the residual population variability, computed as the mean absolute 
difference between successive censuses of a single population, after the removal of 
the newborn; Vy, the yield variability; Py, the mean residual population size when 
F per cent of the newborn are removed at each census; Po, the mean size of a popu- 
lation from which no young are removed; Yr, the mean number of animals removed 
from a population in which the removal rate is F'; F, the percentage of the newborn 
removed at each census. We can assign® the difference between V7 and Vp to the 
fluctuations due to the newborn animals which have been removed. That is, 


Vr ee Vp = Vy. (1) 
The data of Slobodkin and Richman® indicate that 


i Pr 
Vr F+2(1—- FP) 


for all values of F. 
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Assuming the contribution of each animal to population fluctuation to be con- 
stant, equation (2) can be translated into the verbal statement that fluctuations in 
yield are proportional to removal rate of newborn, while total fluctuations are the 
sum of yield fluctuations and fluctuations in the residual populations. Residual 
population fluctuation is proportional to twice the percentage of newborn per- 
mitted to remain in the population. 

The fact that population size is bounded would imply that any change in popu- 
lation size produced by an animal that will remain in the population must even- 
tually be compensated for by an equal change of opposite sign. The animals which 
are removed, on the other hand, contribute only once to population fluctuation. 

More recent data (which will be presented in detail elsewhere) indicate that both 
yield and population size are linearly and directly proportional to food supply over 
a wide range of food increment and removal rate. Some of these data will be dis- 
cussed here in terms of yield and population size per unit food. 

Since Vr is independent of removal rate, we can assume that 


Vp _ Pr 


Vr Po 
Solving equation (2) for Vy, rewriting Vp in equation (3) as V7 — Vy, and sub- 


stituting our solution from equation (2) for Vy, we can write 


= o Vr — VrF id + Ai - ud (4) 
I 0 | T 


whence 


F - 
>, = P pk e 
i (: F + 2(1 — “3 ss 


From equation (5) we can predict population size at any removal rate, given 
the size of the populations at 0 per cent removal. The agreement between calcu- 
lated and observed values of Pr can be seen in Figure 1. 

Note that the number of newborn animals is computed at each census as the 
difference between the sum of living and dead animals found at the current census 
and the number of live animals returned to the container at the previous census. 
There is reason to believe that some of the dead animals decompose completely in 
the interval between censuses, producing an overestimate of a priori removal rate. 
The populations shown in Figure 1 were maintained at a priori removal rates of 
50 or 90 per cent. Utilizing the fact that occasionally, when no newborn animals 
are being produced, a certain number of animals cannot be accounted for, it was 
possible to estimate the percentage of the dead animals that decompose completely 
between censuses for each population, and adjust the estimate of actual removal 
rate accordingly. This has been done in Figure 1. Adjustment was not possible 
at 90 per cent removal, since newborn were always present. 

The number of animals removed from the populations can be thought of as 
the “yield” from the populations. On the simplest assumption, yield from any 
self-regulating fishery should be equal to a conversion factor multiplied by the 
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difference in maintenance energy requirements between a fished and an unfished 
population, that is, 


Yr = K(Po — Pr), (6) 


where K is a conversion factor, in this case equal to 0.36. 

The agreement between observed and calculated values of yield is shown in 
Figure 2. The aberrant points (ft) are from two populations in which mean yield 
had not yet been stabilized. One of these was maintained at such a high food 
level that toxic effects of the algae may have been significant. 


Yield/unit food 














Fig. 2. 


If the removal process is altered so that the number of animals removed remains 
a constant proportion of the newborn, but adult animals are removed instead of young, 
neither yield nor population size is directly proportional to food supply, the pre- 
dictive equations are no longer valid, and it. can be demonstrated that an appre- 
ciable portion of the food given to the populations is never actually consumed. 
The above equations therefore seem relevant only if the food is limiting. 

Summary.—If a constant proportion of the newborn are removed from self- 
regulating populations of Daphnia, population size and yield are predictable from 
simple equations requiring only knowledge of the size of a normal population and 
one conversion constant. When adult animals are removed, destroying the simple 
food dependence of the population, these equations are no longer applicable. 


* These experiments were supported by funds from the Rockefeller Foundation (GA BMR 
5661) and the Phoenix Memorial Project of the University of Michigan (Project 993). Miss 
Mimi Ho did most of the very laborious population counts. 
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